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MEMORANDUM FOR DTIC (Acquisition) 

(Attn; Pat Mauby) 

SUBJECT: Distribution of USAF (AFOSR Summer Research Program (Air Force 
Laboratories) and Universal Energy Systems, Inc., and the Research Initiation Program 

FROM: AFOSR/XPT 
Joan M. Boggs 

110 Duncan Avenue, Suite B115 
Bolling AFB DC 20332-0001 

1. All of the books forwarded to DTIC on the subjects above should be considered 
Approved for Public Release, distribution is unlimited (Distribution Statement A). 

/ 

2. Thank ypti for processing the attached information. 



Chief, Technical Information Division 
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Master Index For Graduate Students 


AhoAd, ZBAd 
BS 

Blactxical Enginscring 
Portland Stata Onivorsity 
Portland, OR 97207-0000 

AnaoliBO, Andrew 
MS 

Mechanical Engineering 
Coluasbia University 
Mew York, NY 10027-0000 

Bapty, Theodore 
MS 

Box 1649 Station B 
Vanderbilt Oniversity 
Nashville, TM 37235-0000 


Field: Electrical Engineering 

I,aboratory: PL/WS 

Vol-Page Mo: 8-19 

Field: Mechanical Engineering 

Laboratory: RL/ER 

Vol-Page Mo: 9-3 


Field: Electrical Engineering 

Laboratory: AEDC/ 

Vol-Page Mo: 11- 1 


Barber, Brian 
BS 

Electrical Engineering 
Washington University 
St. Ijouis, MO 63130-0000 


Field: Electrical Engineering 

Laboratory: RL/OC 

Vol-Page Mo: 9-10 


Bard, David 
MS 

Mechanical Engineering 
Columbia University 
Mew York, NY 10027-0000 

Barnett, William 
BS 

CLOUD & Aerosol Science L 
University of Missouri 
Rolls, MO 65401-0000 

Bartell, Craig 
BS 

Chemistry 

Wright State University 
Fairborn, OB 45435-0000 

Bartelt, Timothy 
MS 

Physics 

Clarkson University 
Potsdam, NY 13676-0000 


Field: Applied Mechanics 

Laboratory: PL/SX 

Vol-Page No: 8-14 


Field: Chesdstry 

Laboratory: PL/LI 

Vol-Page No: 8-6 

Field: Chemistry 

laboratory: WL/PO 

Vol-Page No: 10-26 

Field: Physics 

Laboratory: RL/ER 

Vol-Page No: 9-4 
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GSRP Participant Data 


Brocklehurat, William 

Field: 

BS 

Aeroapace Engineering 

Laboratory: 

Oniveraity of Cincinnati 
Cincinnati, OB 45221-0000 

Vol-Page Mo: 

Buma, Paul 

Field: 

BS 

Electrical Engineering 

laboratory: 

Auburn Oniveraity 

Auburn, AL 36849-0000 

Vol-Page Mo 

Carlen, Edwin 

Field: 

BS 

Electrical Engineering 

Laboratory: 

Oakland Oniveraity 

Rocheater, MI 48309-4401 

Vol-Page Mo: 

Carter, Charity 

Field: 

MS 

Electrical Engineering 

Laboratory: 

Stevena Xnatitute of Tech. 
Hoboken. KJ 7030-0000 

Vol-Page Mo: 

Caalin, Barry 

Field: 

BS 

Coa^uter Science/Engineer 

Laboratory: 

Wright State Oniveraity 

Dayton, OH 45435-0000 

Vol-Page Mo 

Charley, David 

Field: 

BS 

Electrical Engineering 

Laboratory: 

Oniveraity of Cincinnati 
Cincinnati, OB 45221-0030 

Vol-Page Mo 

Cianeroa, John 

Field: 

BS 

Paychology ' 

Laboratory: 

Califoraia State Oniveraity 

Loa Angelea; CA 90032-0000 

Vol-Page No 

Craig, Ralph 

Field: 

BS 

Mathematica: 

Laboratory: 

Worth Carolina State Oniv. 
Raleigh, WC 27695-8205 

Vol-Page No 


n 


A«jro«pao* XnginMring 
PZ./aK 

8-10 

Bl*etxical Engineering 
HL/MH 

10-21 

Electrical Engineering 
HL/AA 

10 - 1 

Xnterdiaeiplinary Engr. 
RX./IR 

9 - 8 

Coaputer Science 
WL/m. 

10-14 

Coa^utcr Science 
WL/AA 

: 10 - 2 

Payehology 

RL/XP 

7- 8 

g 

Hatheaatica 

Pt/VT 

8- 15 



GSRP Participant Data 


Cttxtxa, Josaph 
MS 

Zoology 

Onlvaroxty of CallfozBxa 
Davis, CX 95S16-0000 

Osrnoll, Julio 
BS 

Mochsnicsl Enginooring 
Univoxsity of Tonnossoo 
TullshoBS, TN 37388-0000 

Davis, Wyatt 
BS 

Mschanical Enginooring 
Washington Stato Dnivorsity 
Pullman, WA 99164-2920 

DoVilbiss, Alan 
BS 

Physics 

Dnivorsity of Colorado 
Colorado Springs, CO 80901-7150 

Dotwilor, Duano 
BS 

Aoronautical Engineering 
Ohio Stato Dnivorsity 
Columbus, OB 43235-0000 

DiPietro, Jr., Anthony 
BS 

Aerospace Engineering 
Virginia Polytechnic Institute 
Blacksburg, VA 24060-0000 

Dobransky, Mary 
MS 

SystexBS Science 
Binghamton Dnivorsity 
Binghamton, MY 13902-6000 

Doss, Ellen 
BS 

Biophysics 

Dnivorsity of Scranton 
Scranton, PA 18505-0898 


Field: Biomedicine 

Laboratory: AL/OE 

Vol-Pago Ho: 7-25 

Field: Mechanical Engineering 

Laboratory: AEDC/ 

Vol-Page Ho: 11- 2 

Field: Mechanical Engineering 

Laboratory: FJSBL/ 

Vol-Page Ho: 11-10 

Field: Physics 

Laboratory: FJSRL/ 

Vol-Page Ho: 11-11 

Field: Aerospace Engineering 

Laboratory: WL/FI 

Vol-Page Ho: 10- 7 

Field: Aerospace Engineering 

Laboratory: AEDC/ 

Vol-Page Ho: 11- 3 

Field: Advanced Technology 

Laboratory: RL/XP 

Vol-Page Ho: 9-15 

Field: Biophysics 

Laboratory: AL/AO 

Vol-Page Mo: 7-1 
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GSRP Participant Data 


OuBois, Travis 
BS 

Hsohanieal Enginaaring 
Old Ooaioion Onivarsity 
Morfolk, VA 23529-0247 

Farquhar, John 
MS 

Instructional Taohnology 
Onivarsity of Gaorgia 
Athens, GA 30602-0000 

Fischar, Varlyn 
BS 

Materials Seianea & Engr. 
Onivarsity of Florida 
Gainesville, FL 32611-0000 

Foster, Robert 
BS 

Mechanical Engineering 
Onivarsity of Iowa 
Iowa City, lA 52240-0000 

Frank, Jonathan 
MS 

Mechanical Engineering 

Tala Onivarsity 

Maw Haven, CT 6520-0000 

Fuller, Joan 
BS 

Chasiistry 

Onivarsity of Alabaaia 
Tuscaloosa, AL ** 0 

Gavora, Mark 
MS 

Instructional Systems 
Florida State Onivarsity 
Tallahassee, FI, 32304-3551 

Gaierman, Robert 
BS 

Mechanical Engineering 
Oniversity of Tennessee 
Tullahoma, TM 37388-0000 


Field: Mechanical Engineering 

XiSboratory: FJSBX,/ 

Vol-Page Mo: 11-12 

Field: Instructional Technology 

Laboratory: AL/BR 

Vol-Page Mo: 7-17 


Field: Physics 

Laboratory: WL/ML 

Vol-Page No: 10-15 

Field: Mechanical Engineering 

Laboratory: ML/PO 

Vol-Page Mo: 10-27 


Field: Mhehanieal Engineering 

Laboratory: WL/PO 

Vol-Page Mo: 10-28 

Field: Chemistry 

Laboratory: FJSRL/ 

Vol-Page No: 11-13 


Field: Instructional Systems 

Laboratory: AL/HR 

Vol-Page No: 7-18 


Field: Mechanical Engineering 

Laboratory: AEDC/ 

Vol-Page No: 11- 4 
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GSRP Participant Data 


Oluck, K«vin 
BS 

Psychology 

Trinity University 

San Antonio, TX 78212- lOOO 

Oottlob, Lawrence 
MS 

Psychology 

Arizona State University 
Tsagpe, AZ 85287-0000 

Griffin, Sti^n 
BS 

Engineering 

University of Texas 

San Antonio, TX 78249-0665 

Grosskopf, Kevin 
BS 

FAC 101 

University of Florida 
Gainesville, FL 32611-0000 

Bam, Patricia 
BS 

Technical Coszounication 
University of Washington 
Seattle, WA 98195-0000 

Hickman, Mary 
BS 

Physics and Astronomy 
University of Wyoming 
Laramie, WT 82071-0000 

Hovey III, Leland 
BS 

Computer Science 

SUNT Institute of Technology 

Utica, MY 13504-3050 

Jain, Sanjay 
BS 

Northeastern Ohio University 
Rootstown, OB 44272-0000 


Field: Psychology 

Laboratory: AL/BR 

Vol-Page Mo; 7-19 

Field: Psychology 

Laboratory: AL/BR 

Vol-Page Mo: 7-20 

Field; Mechanical Engineering 

Laboratory: PL/VT 

Vol-Page Mo: 8-16 

Field: School of Bldg. Constmet 

Laboratory: WL/FI 

Vol-Page Mo: 10- 8 

Field: Comminications 

Laboratory: AL/HR 

Vol-Page Mo: 7-21 

Field: Astrophysics 

Laboratory: PL/LI 

Vol-Page Mo: 8-7 

Field: Mathematics 

Laboratory: RL/XP 

Vol-Page Mo: 9-16 

Field: Matural Sciences 

Laboratory: WBMC/ 

Vol-Page Mo: 11-16 


V 




GSRP Participant Data 


Jmftmxmon. Ellmn 

ns 

Biology 

Trinity Univorsity 

San Antonio, TX 78212-0000 

Jonny, Jaaon 
BS 

Matoriala Scioneo 
Camogic-Mollon Onivorsity 
Pittsburgh, PA 15213-0000 

Jolly, Mohanjit 
BS 

Aoronautics/Astronautiea 

MIT 

Caipbridgo, MA 2139-0000 

Jonos, Claud 
BS 

Blcctrical £ Coaputer Eng 
Onivorsity of SW Louisiana 
luifayotto, LA 70504-3890 

Kolly, John 
BS 

Electrical Engineering 
Onivorsity o£ Florida 
Gainesville, FL 32601-0000 


Field: Biology 

Laboratory: AL/OB 

Vol-Page Mo: 7-26 


Field: Materials Science 

Laboratory: ML/ML 

Vol-Page Mo: 10-16 


Field: Aeronautics/Astronauties 

Laboratory: PL/BK 

Vol-Page Mo: 8-11 


Field: Petroleum Engineering 

Laboratory: RL/C3 

Vol-Page Mo: 9-1 

Field: Electrical Engineering 

Laboratory: PL/LI 

Vol-Page Mo: 8-8 


Kisdaall, Lucia 
MS 

Mathematics 

Worcester Polytechnic Inst. 
Worcester, MA 1609-0000 

Kundich, Robert 
MS 

899 Madison Ave., Ste 801 
Univ. Tennessee-Memphis 
Menphis, m 38163-0000 

Kvasnak, William 
MS 

Mechanical Engineering 
Clarkson Onivorsity 
Potsdam, MY 13699-5725 


Field: Applied Mathematics 

Laboratory: PL/GP 

Vol-Page Mo: 8-1 


Field: Biomedical Engineering 

Laboratory: AL/CF 

Vol-Page Mo: 7-9 


Field: Mechanical Engineering 

Ijdsoratory: AEDC/ 

Vol-Page Mo: 11- 5 
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GSRP Participant Data 


L«ix , John 
BS 

Civil Enginooxing 
Onivoxsity of Mow Oxloana 
Mow Oxloana, lA 70148-0000 


riold: 

l^aboxatoxy: 

Vol-Pago No 


Lao, Daniol 
BS 

Bloetxieal Enginooxing 
Southoxn Xllinoia Onivoxaity 
Caxbondalo, XL 62901-0000 


riold: 

Itaboxatoxy: 

Vol-Pago Mo 


Loiwoko, Roboxt 
BS 

Aoxonautical Enginooxing 
Ohio Stato Onivoxaity 
Columbua, OH 43210-0000 


riold: 

luiboxatoxy: 

Vol-Pago Mo 


Litvin, Koxxy 
MS 

Eloctxical Enginooxing 
Coxnoll Onivoxaity 
Xthaca, MT 14853-0000 


riold: 

I,aboxatoxy: 

Vol-Pago Mo 


Lukox, Stophon 
NS 

Biological Seionco 
Onivoxaity of Alabama 
Tuacalooaa, AL 35487-0344 


riold: 

Xddsoxatoxy: 

Vol-Pago No 


Maatexaon, Gina 
BS 

Induatxial Enginooxing 
Aubuxn Onivoxaity 
Aubuxn, AL - 0 


riold: 

Laboxatoxy: 

Vol-Pago Ho 


McMuxtxy, J. 

MS 

Mechanical Enginooxing 
Louiaiana Tech Onivoxaity 
Ruaton, LA 71272-0046 


riold: 

lAboxatoxy: 

Vol-Pago No 


Menendoz-Baxxoto, Molani 
MS 

Phyaica 

Onivoxaity of Puoxto Rico 
Mayaguoz, PR 681-0000 


riold: 

Laboxatoxy: 

Vol-Pago No 


Xnduatxial Ooaign 
ML/rX 

10- 9 

Eloctxical Enginooxing 
RL/ER 

9- 5 

Aoxoapaeo Enginooxing 
PL/WS 

8-20 

Eloctxical Enginooxing 
RL/OC 

9-11 

Maxino Seionco 
PL/OP 

8 - 2 

Xnduatxial Enginooxing 
AL/OE 

7- 27 

Mechanical Enginooxing 
ML/MM 

10-22 

Phyaica 

PL/6P 

8- 3 
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GSRP Participant Data 


IteaMSsehaitt, Jam 
MS 

El«otxie«l Xnqxnmmxxng 
Polytaehnie Onlwrsity 
rAzningdalAf MT 11735~0000 

Mills, Jeffrey 
MS 

Cheietry 

Xndiena Oniversity 
Bloomington, XM 47405-0000 

Monks, Gery 
BS 

Chemiatry 

Oniversity of Sorsnton 
Scranton, PB 18S10-C000 


field: 

laboratory: Pl/VS 
Vol-Page Mo: S-21 

field: rhemistry 

Xiaboratory: PX,/BS 

Vol-Page Mo: 8-12 

field: rhemistry 

laboratory: AL/OS 

Vol-Page Mo: 7-28 


Moore, Michael 
MS 

Electrical Engineering 
Vanderbilt Oniversity 
Mashville, TM 37235-0000 

Moriarty, Daniel 
MS 

Nuclear Engineering 
MIT 

Cambridge, MA 2139-0000 

Magaraja, Chandra 
BS 

Mechanical Engineering 
Oniversity of Dayton 
Dayton, OB 45469-0000 

Maghski, David 
MS 

Electrical Engineering 
Oniversity of Cincinnati 
Cincinnati, OB 45221-0000 


field: Electrical Engineering 

Laboratory: AEDC/ 

Vol-Page Mo: 11- 6 


field: Mnelear Engineering 

Imboratory: PL/OP 

Vol-Page Mo: 8-4 


field: Mechanical Engineering 

Laboratory: NL/PO 

Vol-Page No: 10-29 


field: Electrical Engineering 

laboratory: NL/ML 

Vol-Page No: 10-17 


Meaffer, Ronald 
BS 

Polymer Science 
Oniversity of Akron 
Akron, OB 44325-3909 


field: Miysioo 

Laboratory: WL/ML 


Vol-Page Mo: 10-18 




GSRP Participant Data 


Waw ll, Tim 
MS 

Physios 

Onivorsity of North Tsxss 
Osnton, TX 7S203-0000 

Nguysn, Sonny 

as 

Blootricsl Knginooring 
Wright Onivorsity 
Osyton, OH 45401-0000 

Hutoson, T63d 

as 

Blootricsl Enginooring 
Wright State Onivorsity 
Dayton, OH 45435-0000 

Olson, Michaol 

as 

Bloctrical Enginooring 
riorida State Onivorsity 
Tallahassee, FL 32315-0000 

Panok, Heather 
MS 

Biology 

University of Scranton 
Scranton, PA 18510-0000 

Petroziello, Joseph 

as 

Biology 

Onivorsity of Scranton 
Scranton, PA 18510-0000 

Pobst, Jeffrey 

as 

Aerospace Engineering 
University of Southern Calif. 
Los Angeles, CA 90089-1191 

Povich, Clinton 

as 

Aerospace Engineering 
University of Kansas 
Lawrence, RS 65045-0000 


Field: Physics 

Laboratory: PL/LI 

Vol-Page No: 8- 9 

Field: Electrical Engineering 

Laboratory: WL/PO 

Vol-Page No: 10-30 

Field: Electrical Engineering 

Laboratory: WL/EL 

Vol-Page Me: 10- 4 

Field: Electrical Engineering 

Laboratory: WL/MN 

Vol-Page No: 10-23 

Field: Biochesdstry 

Laboratory: AL/AO 

Vol-Page No: 7-3 

Field: Biology 

Laboratory: AL/AO 

Vol-Page Ho: 7- 4 

Field: Aerospace Engineering 

Laboratory: PL/BK 

Vol-Page Mo: 8-13 

Field: Aerospace Engineering 

Laboratory: FJSRL/ 

Vol-Page No: 11-14 
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GSRP Participant Data 


RadoBaky, Laon 
MS 

Chanieal Enginaaring 
Colunbia Onivaristy 
Haw York, HY 10027-0000 

Raad, Jason 
BS 

Elactrieal Enginaaring 
Cornall Onivarsity 
Ithaea, MY 14853-0000 

Rnichnayar, Francis 
MS 

Elaetlrieal Enginaaring 
Syrscusa Onivarsity 
Syracusa, HY 13244-1240 

Railay, Danial 
MS 

Physics 

Onivarsity o£ Ala b a m a 
Huntsvilla, AL 35899-0000 

Raistar, Craig 
BS 

Elactrieal Enginaaring 
North Dakota State Onivarsity 
Fargo, HD 58105-0000 

Reustar, Danial 
MS 

Elactrieal Enginaaring 
Onivarsity of Dayton 
Oeyton, OH 45469-0000 

Riasingar, Joyce 
BS 

Eeology/Evolutionaxy Biol 
Onivarsity of Hast Florida 
Pensacola, FL 32514-0000 

Ryan, Arthur 
BS 

Psychology 

Wright State Onivarsity 
Dayton, OH 45435-0001 


Field: Materials Science 

Laboratory: WL/EL 

Vol-Paga Ho: 10- 5 

Field: Physics 

Laboratory: RL/OC 

Vol-Page Ho: 9-12 

Field: Electrical Engineering 

Laboratory: RL/C3 

Vol-Paga Mo: 9-2 

Field: Physics 

I,sborstory: WL/MH 

Vol-Paga Mo: 10-24 

Field: Electrical Engineering 

Laboratory: AL/AO 

Vol-Paga Ho: 7-5 

Field: Electrical Enginaaring 

Laboratory: WL/AA 

Vol-Paga Ho: 10- 3 

Field: Marina Biology 

Laboratory: AL/EQ 

Vol-Page Ho: 7-14 

Field: Psychology 

Laboratory: AL/CF 

Vol-Paga No: 7-10 
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GSRP Participant Data 


Saauid, Mohaamad 
BS 

Machanical Enginaaring 
Onlvaraity of Maw Orlaana 
Haw Orlaana, IJt 70148-0000 

Sehroadar, Mark 
BS 

Elaotrieal Enginaaring 
North Dakota Stata Onivaraity 
Fargo, MD 58105-0000 

Schwing, Ronald 
BS 

Machanical Enginaaring 
Dniveraity of Iowa 
Iowa City, lA 52242-0000 

Saifart, Jaaon 
BS 

Paychology 

Dnivaraity of Dayton 
Dayton, OB 45410-0000 

Saydal XIZ, Robart 
BS 

Aaroapaea Enginaaring 
Dnivaraity of Miaaouri 
Rolla, MO 65401-0000 

Shahidi, Anooah 
MS 

Intelligant Syataaa Studi 
Dnivaraity of Pittaburgh 
Pittaburgh, PA 15260-0000 

Shahrokhi, Kiaball 
BS 

Machanical Enginaaring 
Vandarbilt Dnivaraity 
Maahvilla, TN 37232-0000 

Slatar III, Robart 
MS 

Machanical Enginaaring 
Dnivaraity of Cincinnati 
Cincinnati, OB 45221-0000 


Fiald: Machanical Enginaaring 

laboratory: HL/FZ 

Vol-Paga Mo: 10-10 

Fiald: Elaotrieal Enginaaring 

Laboratory: AL/AO 

Vol-Paga Ho: 7- 6 

Fiald: Machanical Enginaaring 

Laboratory: HL/PO 

Vol-Paga Mo: 10-31 

Fiald: Paychology 

Laboratory: AL/BR 

Vol-Paga Ho: 7-22 

Fiald: Aaroapaea Enginaaring 

laboratory: ML/FZ 

Vol-Paga No: 10-11 

Fiald: Zntalligant Syatana 

Laboratory: AL/BR 

Vol-Paga Mo: 7-23 

Field: Mechanical Engineering 

Laboratory: AEDC/ 

Vol-Paga No: 11- 7 

Fiald: Enginaaring 

Laboratory: HL/FZ 

Vol-Paga Ho: 10-12 
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GSRP Participant Data 


Soeei, Xd 
MS 

Materials Soianea £ Engs 
Onivarsity of Virginia 
Charlottasvilla, VX 22903-2442 

Sodoy, Mauraan 
BS 

Biology 

Trinity Univaraity 

San Antonio, TX 78212-0000 

Soransan, Bryant 
MS 

Elaetrical Enginaaring 
University of Tennaasaa 
Tullahoma, TO 37388-0000 

Spencer, Robert 
BS 

Electrical Engineering 
Cornell University 
Ithaca, MI 14853-0000 

Stadlar. Bathania 
BS 

Materials Seianca 
Massachusetts Inst, of Tech. 
Cnabridge, MA 2139-0000 

Stansbary, Donald 
BS 

Aarospaca Enginaaring 
Univaraity of Missouri 
Rolls, MO 65401-0000 

Starchville, Jr., Thomas 
MS 

Aarospaca Enginaaring 
Pennsylvania State University 
University Park, PA 16802-0000 

Starks, Michael 
BS 

Electrical Engineering 
Boston University 
Boston, MA 2215-0000 


Field: 

Laboratory: ML/ML 
Vol-Paga Mo: 10-19 


Field: Biology 

Laboratory: AL/OE 

Vol-Paga Mo: 7-29 


Field: Mathematics 

I,aboratory: AXDC/ 

Vol-Paga Mo: 11- 8 


Field: Electrical Enginaaring 

Laboratory: RL/OC 

Vol-Page Mo: 9-13 


Field: Materials Science 

Laboratory: RL/ER 

Vol-Paga Mo: 9-6 


Field: Aerospace 

Laboratory: ML/MN 

Vol-Page Mo: 10-25 


Field: Aerospace Engineering 

Laboratory: ML/FZ 

Vol-Page Mo: 10-13 


Field: Elaetrical Engineering 

Laboratory: PL/OP 


Vol-Page Mo: 


8- 5 





GSRP Participant Data 


Stauffar, Josaph 
MS 

Managenant 

Univaraity of Iowa 

Iowa City, lA 52242-1323 


Fiald: 

Laboratory: 

Vol-Paga Mo 


Stangar, Vincant 

as 

Elactrieal Enginaaring 
Onivarsity of Cincinnati 
Cincinnati, OH 45221-0030 


Fiald: 

Laboratory: 

Vol-Paga Ho 


Stohs, Jonathan 

as 

Phyaics 

Univaraity of Haw Mexico 
Albuquerque, KM 87131-0000 


Field: 

Laboratory: 

Vol-Paga Ho 


Stroman, Ronald 

as 

Pharatacy 

Onivaraity of South Carolina 
Coluabia, SC 29208-0000 


Fiald: 

Laboratory: 

Vol-Paga Ho 


Sullina, Taahia 

as 

Cheauatry 

Univaraity of Georgia 
Athena, GA 30605-0000 


Fiald: 

Laboratory: 

Vol-Paga Ho 


Sutcliffe, Ronald 

as 

Cheaiiatry 

Southwaat Texaa State 
San Marcoa, TX 78666-0000 


Fiald: 

Laboratory: 

Vol-Paga Ho 


Swaaney, Haltar 
MS 

Electrical Engineering 
Onivaraity of Maine 
Orono, ME 4469-0000 


Field: 

Laboratory: 

Vol-Paga Mo 


Swindal, J. Field: 

MS Laboratory: 

Applied Phyaica 
Yale Univaraity 
Haw Haven, CT 6520-0000 


Management 

AL/HR 


7-24 


Electrical Enginaaring 
HL/EL 


10 - 6 


Math 

PL/VT 

8-17 


Pharmacy 

AL/AO 

7- 7 


Chamiatry 

AL/EQ 

7-15 


Chamiatry 

AL/OE 

7-30 


Electrical Engineering 
RL/IR 


9- 9 


Applied Phyaica 
HL/PO 


Vol-Paga Ho: 10-32 



GSRP Participant Data 


Taaeillo, Mark 
MS 

Klaetrical EnginMxing 
Binghamton Onivorsity 
Binghamton/ NX 13902-0000 

Thomas , Stuart 
B8 

Engineering 

North Carolina State Univ. 
Raleigh, NC 27695-0000 

Tipton, Kevin 
MS 

Nutrition and Food Sciene 
Auburn University 
Auburn, AL 36849-5605 

Tomow, Carina 
BS 

Quantitive Psychology 
Ohio State University 
Columbus, OB 43210-0000 

Trick, Kimberly 

MS 

Chemical Engineering 
University of Dayton 
Dayton, OH 45401-0240 

Heaver, Michael 
MS 

Aerospace Engineering 
Georgia Institute of Technolog 
Atlanta, GA 30332-0150 

Heitzman, Peter 
MS 

Thayer School of Engineer 
Dartmouth College 
Hanover, NH 3755-0000 

Hhitmore, Christine 
BS 

Mechanical Engineering 
Washington University 
St. l/ouis, MO - 0 


Field: Electrical Engineering 

Laboratory: RL/OC 

Vol-Page Mo: 9-14 

Field: Chemical Engineering 

Laboratory: AL/EQ 

Vol-Page No: 7-16 

Field: Eoology 

Laboratory: AL/CF 

Vol-Page Mo: 7-11 

Field: Psychology 

Laboratory: AL/CF 

Vol-Page No: 7-12 

Field: ChesuLcal Engineering 

Laboratory: HL/ML 

Vol-Page No: 10-20 

Field: Aerospace Engineering 

I,aboratory: AEDC/ 

Vol-Page No: 11- 9 

Field: Electrical Engineering 

Laboratory: FJSRL/ 

Vol-Page No: 11-15 

Field: Mechanical Engineering 

Laboratory: AL/CF 

Vol-Page No: 7-13 


XIV 




GSRP Participant Data 


Wld*MX, Chazlas 
MS 

Electrical Engineering 
Syrecuae Oniversity 
Syracuse, NY 13244-0000 
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PARALLEL TURBINE ENGINE INSTRUMENTATION SYSTEM 


Ted A. Bapty 
Research Faculty 

Department of Electrical Engineering 

Abstract 

The complexity of modern turbine engines requires significant amounts of testing to ensure 
correct, cost effective development. Ground testing enables highly accurate control over the 
external conditions along with the ability to acquire a large number of measurements of internal 
structures. 

.Analysis of turbine engine testing data is a computationally intensive, dynamic process. The 
combination of high bandwidths, large numbers of channels, and constantly changing processing 
requirements place many demands on an instrumentation and data analysis system. 

The Computer Assisted Dynamic Data Acquisition and Monitoring System(CADDMAS), de¬ 
veloped in conjunction with Arnold Engineering Development Center, has been in operation for 
several months in support of simulated altitude turbine engine stress testing. The system uses 
parallel processing to sustain over 400 million floating point operations per second(MFLOPS) 
producing on-line data analysis plots, such as Campbell Diagrams, Phase Campbell Diagrams, 
and spectral energy plots. Until recently, these plots took several days to weeks to obtain. The 
plots are displayed on multiple high speed color displays using standard engineering graphics. 
Hardcopy plots are available on demand. 

I describe hardware and software architecture of a greater than 100 processor parallel instru¬ 
mentation system for turbine engine aeromechanical stress analysis. The hardware consists of 
a mixture of hetereogeneous processing nodes arranged in a configuration closely matching the 
analysis algorithms. Specialized processors are used for numerical, input/output, graphics, user 
interface, and storage. 

The software uses model-based program synthesis techniques to dynammically reconfigure 
signal analysis algorithms to meet changing requirements. 
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PARALLEL TURBINE ENGINE INSTRUMENTATION SYSTEM 


Ted A. Bapty 


Introduction 


During the development and operational lifespan of a turbine engine design numerous tests 
are performed to ensure design integrity and performance over time. Ground testing forms a basis 
for this testing. The engine is installed in a controlled environment test cell to simulate a wide 
variety of conditions. The conditions are adjusted to simulate the entire operating envelope of 
the engine, to predict the operation of the engine in the field. 

Turbine engine stress testing is a critical phase in the testing process. In this test mode, 
the engine is instrumented with hundreds of sensors measuring the stress on various internal 
components. The stress sensors are mounted on blades and stators, directly measuring dynamic 
vibrations. The signals are telemetered out of the engine into the data acquisistion system. In 
addition to these signals, a wide variety of external parameters such as airflow, temperature, 
pressure, etc, are measured. 

Turbine Engine testing generates massive volumes of data at very high rates. A single engine 
may have several hundred sensors, with bandwidths up to tens of kilohertz. Test periods last for 
several hours. The entire bulk of data generated is in the terabytes for a single test. 

Complete analysis of all acquired data has been technically impossible in the past, due to the 
speed and capacity of standard computers. Historically, only a small fraction of the data has 
been fully processed by periodic sampling of the data stream. This analysis was also performed 
off-line, with delays of days to weeks before results were available. This latency prohibits the use 
of results to guide the testing process. 

On-line processing of these signals is critical to improving the testing process. The results 
from the on-line computations should be available to enable mid-test decisions, interactive test 
planning. Furthermore, on-line examination of information can assist in catastrophe avoidance, 
by enabling the detection of pathological engine conditions. 


1-3 




Performing these computations and displaying the results in an interactive format requires an 
extremely high processing rate. Typi. 1 algorithms include frequency trzmsformation, peak detec¬ 
tion and validation and correlation. For the analysis of 24 channels of stress data, approximately 
400 MFLOPS of sustained computation is required. To process the hundreds of channels available 
in a typical test, an aggregate sustained computation rate in the GFLOPS range is necessary to 
meet requirements. Note that these are sustained computational rates, not to be confused with 
Peak MFLOPS advertised by computer manufacturers. The situation is further exacerbated by 
the requirement to acquire the volumes of data into the computer and to plot the results quickly. 
The aggregate data input is in the tens of MWords per second. Compounding these requirements, 
is the real-time requirement that the system must keep up with the data stream and display results 
interactively. 

Typical supercomputers advertise processing speeds in these ranges. Some support data acqui¬ 
sition systems capable of meeting the input bandwidth requirements. The problem with applying 
standard supercomputers to this type of task is that both of these facilities are typically not 
available concurrently. The computer must simultaineously accept and process the data. Super¬ 
computers also tend to be expensive to acquire and maintain, especially with high-bandwidth 
acquisition systems. The real-time nature is complicated by the batch-oriented operating systems 
typically found on supercomputers. 

The majority of algorithms required for turbine engine testing are applied independently on 
each of the channels. The high computational requirements and natural parallelism make parallel 
processing a cost-effective hardware approach. The data streams are physically separate channels, 
so a parallel I/O system is also applicable. 

To achieve the multi-GFLOPS goal, a large number of processors must be applied. I designed 
and built a parallel processing system based on the Inmos Transputer to process turbine engine 
stress data on-line. The bare Transputer however, lacks the computational capabilities required 
for this task. At 2 MFLOPS (Peak) execution speed, over 1000 processors would be required to 
achieve the 2 GFLOPS needed. The I/O capabilities for analog input and graphics output are 
not inherent in the chip, so these features are added, producing a heterogeneous architecture. In 
the section 2, I will discuss the hardware configuration. 

While the hardware is available (most of the components are available commercially) and the 
technology allows easy configurations of arbitrary topologies, the software technology is not as 
simple. The software is complicated by several factors: 
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1. Distributed Control: Instead of having one thread of control executing in the system, 
I have as many threads as there are processors. The problem lies in getting all of these 
threads to cooperate and accomplish a single task. 

2. Communication: As the computational data dependencies are divided across processors, 
there must be a method for transmitting intermediate results from one place to another. 

3. Real-Time: The data is coming into the system at a constant rate. The system must 
respond to all of the data at this rate, or data will be lost. 

4. Dynamically Changing Requirements: The processing requested from the system is 
constantly changing. Each test point has a different mix of required output processing. 

Modularity: Systems of various sizes and processing capacity are required. The compo¬ 
nents (hardware and software) must be designed to allow modules to be plugged together 
to quickly and easily form arbitrary systems. 

These problems must be solved in a way that will not make the system implementation cost- 
prohibitive to implement. In the Software Environment section I discuss the software methodology 
used to implement the parallel instrumentation systems such as the CADDMAS 24 channel pro¬ 
totype system. 


CADDMAS Hardwaie Architecture 


The CADDMAS system is a heterogeneous, massively parallel system designed for parallel 
instrumentation applications. The current prototype relies on the Inmos Transputer for the basic 
processing nodes of the architecture. 

A Transputer is used in conjunction with each of the processing nodes. It provides the general 
message passing and control fabric of the system. The basic nodes are augmented as necessary 
to get the necessary computational characteristics. The basic nodes are as follows: 

1. The Transputer Processing Element uses a T800 and dynamic RAM (1 or 4 MBytes) for 
general purpose data manipulation tasks, such as data accumulation and plotter display list 
generation. 
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Figure 1: Real-Time Data Visualization System 

2. The Numeric Processing Element uses a T800 and a Zoran floating point Vector Signal 
Processor to supply the bulk of the floating point computational horsepower. 

3. The Graphics Processing Element consists of a IBM/386 compatible personal computer with 
an accelerated graphics card. The PC also provides interfaces for mouse and keyboard input 
and limited digital storage. 

4. Front End Processor samples the analog data, and detects peaks, filters and downsamples 
the digital data. 

The CADDMAS 24 channel prototype application of these processors to on-line test data 
processing requires a mix of approximately 100 processors. Figure 1 shows a typical full-scale 
system. 

The next two subsections describe in further detail the custom-designed components. For 
detailed descriptions of the commercially-available components, please refer to the product liter¬ 
ature. 
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Numeric Processing Element 


The Numeric Processing Element uses a 20 MHz T800 for communications and control. The 
Transputer has access to 1024K Bytes of static memory. The Zoran Vector Signal (VSP) Processor 
can also be a bus master to access this memory. 

The VSP is modeled after an array processor architecture. The VSP has high speed internal 
memory to allow concurrent processing with the host T800. Internal organization of the VSP 
allows concurrent operation within the subsystem, to further increase processing speed. 

The VSP implements IEEE standard 32 bit floating point. Its instruction set includes high- 
level routines that support common signal processing tasks. Vector addition, multiplication, 
polynomial e.\pansion, min/max. accumulate, int to FP, and FP to int are single instructions. 
There are also instructions that implement FIR and HR filtering. Frequency transformation 
primitives (DFT and FFT butterfly) are also single instructions. 

The VSP can perform real 1024 point FFT’s in 2.25 ms, complex in 2.35 ms, and can perform 
logaritms and complex magnitudes in 2.1 and 1.9 usee. By comparison, the T800 takes 100 msec 
for the complex FFT. 

The features and performance of the VSP combine to produce a node architecture that is 
very efficient for performing standard signal processing operations. By comparison, the T800 is 
more than an order of magnitude slower. This high speed node aJieviates the need to subdivide 
algorithms like the FFT, thus reducing the required number of processors and decreasing system 
complexity. 


Front End Processor 


As the transputer is an entirely digital processor, it has no method for acquiring data from the 
external world. This data acquisition is the primary function of the Front End Processor. Data 
is first sampled by a high speed Analog/Digital converter. Two Motorola 56001 Digital Signal 
Processors then pack the data into data ensembles and scan the data for maximum and minimum 
stresses. The data is filtered and downsampled to the Nyquist rate for the system. The data is 
transmitted out an Inmos-compatible serial link for processing in the CADDMAS system. 

The functions of the FEP are divided into two separate boards to improve the isolation of the 
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analog signal from the noisy digital system. The boards are separated using opto-isolators at all 
interface points. 


Software Environment 


While hardware can and ha5 been constructed to meet these requirements,the software engi¬ 
neering can be very difficult and expensive. The problem is compounded by the requirements for 
real-time performance, parallel computation, flexibility, and in situ reconfiguration. 

If standard software engineering techniques were used, the problem would be broken into 
manageable segments and partitioned out over the hardware architecture. Interfaces would be 
defined between each of these modules and communication protocols would be constructed. At 
that point, the coding of components and communications routines could commence. The major 
problems would occur when system integration begins. At that point, any overlooked factors 
would become apparent. Interfaces would change, possibly causing changes to ripple throughout 
the entire system. These changes would continue to ripple throughout the system when the actual 
performance was measured, and when requirements varied slightly. This would be an extremely 
painful method to attempt the development of such a system. It is doubtful that the development 
effort would end in success. 

Vanderbilt researchers have been working with AEDC for the past five years to develop tools 
and techniques to manage the software complexity of this large, parallel, real-time instrumentation 
system. In order to address the problems associated with the development of CADDMAS, I have 
used model-based techniques coupled with a real-time macro-dataflow execution environment. 
The principal technique is to automatically generate a macro-dataflow computation structure 
from a declarative model describing the desired system. This new technique, introduced by the 
Measurement and Computing Systems Group of Vanderbilt, is well suited to a large class of signal 
processing, instrumentation, and control problems. The primary tools and techniques used here 
are part of the Multigraph Programming Environment. Using the model-based approach provides 
several attractive features; 

• Program Synthesis techniques allow the system complexity to be managed in the model- 
based format. A simple example of this is a signal processing system. In this case, a 
model is constructed pictorially which looks like a signal flow diagram. When the system is 
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ready for execution, that model is analyzed and the corresponding real-time dataflow graph 
is constructed and started on the execution platform. The algorithm is specified at the 
highest level. Implementation specific details are handled by the builder. Changes in the 
underlying hardware architecture or operating systems have no effect on the models. The 
models are portable across many platforms. 

• Builder/Execution Decoupling The builder is only loosely coupled to the execution 
system. None of the signal-flow algorithms interact with the builder, and therefore the 
builder is not a bottleneck during execution. 

• Scheduling, Synchronization.and Communication are implicit in the modeling paradigm. 
The system is synthesized from a library of executable subroutines written in sequential 
languages. 

• Dynamic Reconfiguration during execution of a dataflow graph is supported. One scenario 
is that conditions in the currently executing system reach a point where they require a 
structural adaptation of the system. These conditions re-trigger the program synthesis. The 
builder reconstructs the system to better match the new conditions. Since the executing 
graph is real-time and continuous, it cannot be stopped so these reconfigurations occur 
on-the-fly. 

The Model-Based Approach to Instrumentation 

Development of real-time instrumentation systems requires the crossover of two distinct fields 
of engineering: signal processing and computer engineering. The signal processing issues are 
mainly high-level, algorithm dependent issues while the computational complexity of the real-time 
events and parallel processing synchronization causes numerous low level, computer engineering 
issues to emerge. 

I have combined numeric and symbolic processing to ease the crossover of the two disci¬ 
plines [l]. The Multigraph Architecture provides a general framework aiding the construction of 
model-based engineering systems for real-time, parallel computing environments [2]. The main 
components of the architecture are shown in Figure 2 and described below; 
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Figure 2; The Muitigraph Architecture. 

1. The Symbolic Model Declarations contain modeling information encoded in a problem spe¬ 
cific declarative language. For signal processing systems, these declarations are typically 
of a form that describes the signal flow graph topology as well as specific parameters of 
the processing blocks in the diagram. On multiprocessor systems, models of the available 
hardware resources are also included. 

2. The Model Builder interprets the symbolic model declarations. It creates, initializes, and 
connects the computation nodes of a macro-dataflow graph. This executable dataflow graph 
is built as specified by the Symbolic Model. In multiprocessor systems, an arc of a graph 
that crosses processor boundaries automatically sets up inter-processor communication. 

3. The Library of Operations contains basic sequential programs to be assigned as the compu¬ 
tation scripts of the specific nodes of the dataflow graph. For signal processing applications, 
these are typically numeric routines implementing basic computations (e.g. a FFT or filter). 
These components are usually implemented in standard numeric languages: C, or Fortran. 

4. The Execution Environment (also called the Multigraph kernel) provides facilities allowing 
the model builder to dynamically configure the macro-dataflow graph. As well, the kernel 
provides a scheduler that controls the graph's execution. The kernels present themselves to 
the model builder as a single virtual machine even when a set of heterogeneous distributed 
processing elements are used. 

5. The Hardware and Communication Systems provide resources to the execution environment 
needed to run the application as well as a communication path between processors. The 
Multigraph kernel insulates the user from this level. 
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Extensions to the Dataflow Model 


The macro-dataflow environment supported by the Multigraph kernel has several unique ex¬ 
tensions; 

1. Actornodes are the computational operators of the dataflow graph. An actornode has a state 
associated with the kernel scheduler. It may be inactive, active, ready, or running. An 
inactive actor will never be chosen for execution by the scheduler. An active actor may 
be chosen but does not have the proper input control tokens to warrant execution at the 
preset time. A ready actor is waiting in a queue for the processing resources to become 
available for it to be executed. A running actor is currently being executed. Functions to 
create, destroy, connect input and output ports to datanodes, disconnect from datanodes, 
activate, deactivated, stop, check the status of. replace, set the context of, set the script of 
actornodes are all provided for the symbolic interface. An actornode itself is made up from 
several components: 

(a) The script is a piece of code which performs the required operation on data items prop¬ 
agated to the node. It may be coded in almost any language, symbolic or numeric. 
The script is written in a reentrant fashion so that it can be attached to more than 
one actor simultaneously. Special Multigraph kernel calls are available to a script pro¬ 
grammer for receiving the input control tokens and propagating the resulting token(s). 
As well, special memory management routines are provided by the kernel for increased 
performance. Further kernel calls are available for error handling within the kernel. 

(b) The context is a static local memory section for a particular actor. It provides two basic 
services for the actornode. First, it allows the script to have a place to save its state. 
Thus the actornodes functionality may be changed due to previous values received as 
input control tokens. Second, it provides an interface to the symbolic builder and 
control functions such that parameters may be set during the build operation that will 
affect the operation of the script during execution. For example, a particular signal 
processing control graph may have an actornode that receives an input vector of time 
domain data samples, performs a fast fourier transform, and then propagates the result. 
One implementation of such an actor might be to provide the radix and input data 
packet size information at build time to the actor through its context. In this case, the 


1-11 







context provides a simple interface for the model builder layer to modify the actor’s 
behavior. 

(c) The number of input and output ports are assigned to an actor when it is created. The 
ports themselves are the location passed to the kernel routine when an input is to be 
received or an output is to be propagated. For example, a cross correlation actor may 
receive two blocks of time domain samples, one from ea^h of the data streams to be 
correlated. Stream one is connected to one input port and stream two to the other. 
The result is propagated out the output port. Thus, this actor has two input ports 
and one output port. An actor with no input ports is always ready for execution. 

(d) The control principle determines what criteria will be used by the kernel scheduler to 
decide when an actor’s script should be executed. There are currently two choices for 
this, ifall and ifany. Under ifall criteria, the actor will only be scheduled for execution 
after all input ports have data available. The ifany criteria causes the actor’s script to 
be scheduled any time one or more data items are available for input. 

2. Real-time actornodes are a special variation of actornode. They have two scripts. The first 
script is executed under the same criteria as a normal actor. Before completing, this first 
script should make a special kernel call to inform the kernel that the actor is waiting for a 
particular event to occur (such as an I/O channel communication to complete). The actor 
will consume no more CPU time until the event occurs. The actor will then be restarted 
after the event specified, at which time the second script will be executed. 

.3. Datanodes provide a queuing and connection function between actornodes. In order to allow 
certain portions of the graph to have a higher priority than others, the large-grain dataflow 
system implemented by the Multigraph kernel is asynchronous. The asynchronous nature 
of the graph forces the kernel to provide dynamic scheduling functions. As well, a way to 
queue input tokens between actornodes must be provided. This queuing function is provided 
by passive kernel components called datanodes. Any number of actor output ports may be 
connected to a datanode. A datanode may be connected to any number of actor input ports. 
Datanodes may be enabled or disabled. Unless enabled, a datanode will not allow data to 
be propagated through it to an actor’s input j- ^rt. Datanodes have been purposely left as 
a separate component from actors rather than providing the queuing function by way of an 
actor’s input ports. Datanodes provide a simple interface for the symbolic layer to control. 
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build, and monitor an executing graph. Functions to create, destroy, connect to actor input 
and output ports, disconnect from actors, clear, read, write, enable, disable, check length, 
set length, and check connections of datanodes are all provided for the symbolic interface. 

4. Environments are used to protect system resources and provide a priority mechanism for 
sections of the dataflow graph. All actors are assigned to an environment. Each envi¬ 
ronment is given a numeric priority. The kernel scheduler always prefers actors connected 
to environments of a higher priority. Environments of the same priority are serviced in a 
round-robin fashion. Only one actor per environment will be executed at any one time. 
Thus, if two actors need to share a global resource, placing them in the same environment 
ensures mutual exclusion. Functions to create, destroy, attach to tasks, and set priority of 
environments are provided for the symbolic interface. 

Tasks provide a generic interface to the basic computational resources of the underlying 
machine. In a multitasking environment they are simply the different processing threads 
available to the Multigraph kernel. In a multiprocessor system they are the individual 
processors themselves. Environments are attached to tasks. The task and environment 
concepts allow the actual implementation of the underlying system to be hidden from the 
user. Thus, the exact same model-based graph building techniques may be used regardless 
of the underlying hardware architecture. Functions to create and destroy tasks are provided 
for the symbolic interface. 

In summary, the Multigraph kernel provides an interface to build, modify, monitor, and control 
an asynchronous macro-JtAaflow graph. The graph may be partitioned by the concept of environ¬ 
ments and tasks. The ac:ual dataflow graph looks slightly different from that of a classic dataflow 
graph due to the introduction of the passive nodes providing queuing functions, datanodes. The 
computation nodes are called actornodes. They have input ports, output ports, a context, a 
script, and a status. The use of distinct actornodes and datanodes provides the capability for the 
symbolic system to build, control, and monitor the dataflow graph. 
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Signal Flow Graph Modeling 


In order to assist the creation and management of a particular dataflow graph, various tools 
and languages ranging from basic textual input to a graphical editing environment have been 
developed [3], [4], [5]. 

The Multigraph tools HDL and GMB are frequently used to model signal processing systems. 
HDL stands for Hierarchic ;! Design Language and GMB stands for Graphical Model Builder. The 
purpose of the language is to represent signal processing systems in a non-algorithmic manner. 
Rather, the system is d* ribed in terms of its structure. As such. HDL is classifled as a declarative 
language. Block diagrams are often used to represent the signal flow among various modules of a 
complex signal processing system. 

The HDL declarations allow a simple means of describing these blocks and their intercon¬ 
nections. .A.S well, HDL adds an extension to this approach through the concept of hierarchical 
decomposition. .4ny one block of the modeled signal flow diagram may itself be represented by 
a set of blocks (a sub-component, a smaller signal-flow diagram). These hierarchies may proceed 
to arbitrary levels of abstraction. Complex blocks, made from one or more other blocks, are 
called compounds. Non-decomposable blocks are called primitives. The actions to be carried out 
by primitives are implemented with Multigraph actomodes. Each primitive has local parameters 
that have been inherited/derived through the hierarchy. After design of the top level signal flow 
graph of the hierarchy is complete, an interpreter (the Model Builder) converts it first to the 
primitive components of the hierarchy and then maps the result into the parallel execution envi¬ 
ronment provided by the Multigraph kernel. Each of the nodes of the resulting dataflow graph has 
a local state derived from the parameters inherited through the hierarchy (the actor’s context). 
Node neighbors are known only through the input and output connection arcs. The resulting 
execution ready program is inherently parallel as a result of the declared signal flow graph being 
parallel. As well, it has been automatically generated from the high level declarations of the 
signal processing system and may now be scheduled by the underlying environment. Thus, the 
computer engineering issues around synchronization have been removed. 

GMB is a generic iconic editor program. It runs under X windows and on PC platforms. 
GMB has been tailored in correspondence with the HDL language. Through this interface, the 
user may build the declarations of his desired signal flow block diagram in an entirely graphical 
manner. Figures 3 and 4 show screen dumps from the GMB environment. 
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Figure 4: The CADDMAS multiplexor signal flow graph. 
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Hardware Architecture Modeling 


In order to assist the creation and management of a parallel architecture suitable for a partic¬ 
ular application system, the Automatic Parallel Network Analyzer (APNA) tool has been devel¬ 
oped. The basic computational blocks of parallel instrumentation systems often allow processor 
interconnections that are not fixed. For example, transputer architectures with DSP coprocessors 
fit signal processing applications especially well since the hardware connection architecture can 
be configured to match the structure of the signal flow graph. 

As computational requirements increase, so does the need for more processors. When an ar¬ 
chitecture contains more than a couple of processors, managing hardware complexity becomes a 
serious problem. Testing to see if the network of processors has been interconnected correctly, 
searching for hardware errors, generating network information for an application loader, or creat¬ 
ing map files for a message passing system by hand is time-consuming for tens of processors and 
impossible for hundreds. 

Modeling the hardware architecture helps manage this complexity. Our approach again uses 
the generic graphical model builder GMB. Models of parallel hardware architectures can easily be 
expressed in graphical terms (i.e. by icons and their interconnections). Hierarchy based modeling 
is used to manage complexity. The lowest hierarchy levels model the basic processors while the 
highest level constitutes the complete system. Intermediate levels model boards and subsystems. 
A number of attributes are attached to the processors (e.g. speed or memory size) in order to 
capture additional aspects of the system. 

The declarative description of the hardware architecture is used to automatically produce load 
files and maps required by different software components of the system. The APNA tool loads 
the interpreted models and draws the network in a graphical window. It is capable of loading files 
of a number of different formats, comparing networks and displaying their differences, generating 
files for hardware diagnostic programs and a network loader, and displaying the attributes of the 
processors. It can selectively generate message passing maps, handling multiple paths between 
processors while optimizing the maps for minimum memory usage. Figures 5 and 6 show screen 
dumps from the GMB used to model hardware architecture and the APNA tool. 
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Figure 5: The top level CADDMAS hardware model set up in GMB 
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I'he CADDMAS 

The CADDMAS system was developed using the hardware and software methodology pre¬ 
viously described. The execution system was synthesized using the models pictured in 3. Each 
horizontal block represents the processing associated with a single Front End Procesor. The 
hierarchical decomposition of these blocks is shown in 4. Any changes to the processing per¬ 
formed on all signals is accomplished by modifying the structure of this or lower-level blocks. The 
hierarchical decomposition allows a top-down approach to these modifications. 

To construct a system of arbitrary size, these blocks are replicated in the editor to match the 
desired number of channels. The blocks are connected to the graphics block to enable display and 
printing of the processed results. 

.\t the same time, the hardware is modeled in the hardware modeling paradigm. The pro¬ 
cessor complexes, i.e. physical boards, are hierarically modeled as small blocks. Given these 
blocks, the overall architecture is modeled. The model interpreter uses these models to construct 
communications maps, processor load configurations, and test programs for the new architecture. 

Another program, the User Interface(UI), performs two critical functions in the CADDMAS 
system. The graphical user interface allows the user to configure various visualization screens 
interactively. The user can select the number of visible windows on a screen, the contents of 
each plot window, and the parameters of each plot, such as titles, labels, axis ranges, and plot 
type, and display window update rate. Stored configurations automate the operation of the user 
interface. The user can also print any window or all windows on a screen. 

The UI also manages the online reconfiguration of the executing system. As the user request 
plots, the computation structures must be created from the system models, and unused structures 
must be destroyed to keep from saturating system memory and processing resources. 


Conclusions 


The model-based approach has proven to be a very useful tool to help manage the complexity 
of this large, parallel system. I have now built three different sized CADDMAS systems using 
the procedure outlined above. The first system was a four channel (single FEP), 8 Campbell, 2 
graphics displays system. Actual configuration time for that first system could not be measured 
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because some of the tools were being finished along the way. The second system was a 24 channel, 
48 Campbell, 3 graphics displays, 3 laser printers system. It took one week to physically wire, 
and get the software running. That time included finding a previously missed bug in the message 
passing system and maJcing several new patch cables. The third system was a four channel (single 
FEP), 4 Campbell, 1 graphics display system made from spare parts. It only took four hours to 
physically wire and get the software running. This proves that the approach is indeed flexible and 
efficient. 

Future plans for the tools and techniques include; better methods to map the signal flow graph 
to the hardware, provide a dynamic load balance capability, automatically create and destroy 
pieces of the signal flow graph as required by user interface changes (this addition will greatly 
help the need for cross correlation operations since pre-building all possible cross correlations 
would require (N * (N - 1)) connections). 

Future plans for CADDMAS applications include a 200 channel 50 KHz system using Texas 
Instruments TI 320C40 processors in conjunction with some of the current transputer technology. 
The Multigraph, .APNA, and HDL tools are sufficiently portable that the added processor type 
should not be a major problem. 
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Abstadt 

Arnold Engineering Development Center (AEDC) wants to initiate ground test 
capabiUty of hydrogen-fueled aeropropulsion systems in the Aeropropulsion Systems Test 
Facility (ASTF) using a freeiet test configuration. Such testing vriil potentially produce 
combustible mixtures of hydrogen and air entering the exhaust ducting, which represents 
serious safety considerations. Reduction of the level of hydrogen in the exhaust ducting to 
below flammability limits requires a Hydrogen Disposition System (HDS). AEDC is 
investigating the technical feasibiUty of an HDS for ASTF that involves controlled burning of the 
exhaust gas using pUoted, vee-gutter flameholders. Due to defidendes in the literature and the 
hydrogen-air flammabUity database, AEDC has initiated a technology development program to 
investigate the flammability limits of premixed, turbulatt hydrogen-air flows over the proposed 
HDS hardware for the range of flow ccmditions in ASTF where flammability is uncertain. 

To predict the flame stability performance of the proposed HDS hardware, the author 
modeled the recirculation zone of the vee-gutter flameholder as a well-stirred reactor. The well- 
stirred reactor is a one-dimensional model for steady, homogeneous combustion. Application ot 
the conservation equations aiKl the chemical kinetics of hydrogen-air combustion resulted in a 
set of non-linear, algebraic equations that were solved numerically using the Lewis General 
Chemical Kinetics and Seivsitivity Analysis Code, LSENS. The flame stability limits predicted by 
LSENS will be compared to DeZubay's experimental database of hydrogen flammability limits 
and to the HDS technology development database when testing is complete. Results of toe 
comparison of toe analytical predictiorw and the experimental data will be reported when toe 

analysis is concluded. 
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Using a Well-Stirred Reactor M(^el 

Julie A. Darnell 


1.0 Introduction 

Arnold Engineering Development Center (AEEXI) wants to irutiate ground test 
capability of hydrogen-fueled aeropropulsion systems in the Aeropropulsion Systems Test 
Facility (ASTF) using a free-jet test configuration. Such testing will produce combustible 
nuxtures of hydrogen and air in the facility exhaust ducting. Reduction of the level of hydrogen 
in the exhaust ducting to below flammability limits requires a Hydrogen Disposition System 
(HDS). The HDS will increase the engine operation times over the test envelope before 
combustible levels of hydrogen are produced in the exhaust ducting. A coiKeptual design of an 
HDS for ASTF involves controlled burning of the exhaust gas using an array of piloted, vee- 
gutter flameholders, sunilar to those in a gas tuihine afterburner. Figure 1 is a simple schematic 
of the conceptual design for real-time afterburning. 



H2-Fuelecl 

Engine 



Atm 


Figure 1; Schematic of the HDS Conceptual Design for Real-Time Afterburning 1^1 

The facility and capital requirements for this design are feasible; however, technical feasibility of 
the design over the range of conditions expected in ASTF is uncertain.l^l 

Flammability limits of hydrogen-air mixtures are well documented at pressures near 
atmospheric;l^l however, the literature is deficient concerning hydrogen-air flammability limits 
on the order of 1.0 psia, the minimum operating pressure of ASTF. Due to these deficiencies. 
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AEDC is conducting a technology development program that involves a hiU-scale, section test of 
the proposed HDS hardware. This investigation will provide a technological database for the 
actual design of an HDS for ASTF. The database will detail the flame stability limits of 
ptemixed, turbulent hydrogen-air flows over the piloted, vee-gutter flameholder for the range of 
flow conditior\s in ASTF where flammability is uncxrtain.t^l 

The author investigated the viability of predicting the flame stability perfomumce of the 
proposed HDS hardware using a well-stirred reactor model. The model treats the recirculation 
zone of the vee-gutter flameholder as a well-stirred reactor. The govenung equations of the 
model, which incorporate chemical reaction kinetics arul the conservation equations, were 
solved using LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code.I 
When the analysis is conduded, the flame stability limits predicted by LSENS will be compared 
to predictions based on DeZubay's data arid to the experimental results from the HDS 
technology development test. 

2.0 HDS Technology Development Test 

AEDC personnel designed the HDS technology development test to investigate the 
flanunabiiity of hydrogen-air mixtures over the range of test conditions expected in ASTF where 
flammability is uncertain. The test fadlity, flow simulation requirements and test procedures 
are described below. 

2.1 Test Facility Desaiption 

AEDC is conducting the HDS secticm test in RlAl, a research test cell located in 
the Engine Test Facility (ETF). Figure 2 is a schematic of the HDS technology 
development test setup in RlAl. The primary airstream, supplied by atmospheric or 
high pressure air, flows downstream. The injector rings introduce gaseous hydrogen 
into the primary flow. The mixing duct provides a homogeneous, turbulent mixture of 
hydrogen and air at the flameholders. The secondary air supply produces a layer of air 
around the primary flow to prevent flashback upstream through the boundary layer. 

The two outer flameholders are centered over the primary and secondary flow and 
provide the proper boimdary conditions for the center flameholder. The center 
flameholder is immersed in tite primary flow and is the focus of study. Hydrogen-air 
pilot torches situated behind each flameholder provide ignition of the cold flow. 

Fig\rre 3 shows a detailed view of the test section. Upstream of the 
flameholders, the cross section of the primary flow is 4 in x 11 in. The spacing of the 
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Figure 2: HDS Technology Development Test Setup in R1A1 
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Figure 3: Detailed View of HDS Test Section 
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flameholders produces thirty percent geometric blockage at the flameholder minimum 
area. The geometry of the combusticm chamber provides a flow separation anchor and 
approximately cor\stant static pressure immediately dowitstream of the flameholders. 
Diagnostic systems include an infrared wavelengflt radiometer to determine flame light 
off and blowoff and a white light shadowgraph system for flame visualization. 

The flameholder geometry, shown in Figure 4, was desigrted as a typical vee- 
gutter flameholder in a gas turbine afterburner. 



Figure 4; Geometry of the Vee<5utter Flameholders 

The flameholder has a widfo of 15 in with a 30° included angle. The flameholder 
extends across the duct; however, the length of the flamdiolder exposed to tt\e primary 
flow is 4 in. The water cooling duct transfers heat from the flameholder surface. 

The exhaust system is designed to cool the hot combustion products, reduce foe 
volumetric flow rate, inert foe exhaust gas with air to insure foe hydrogen level is below 
flammability limits, and direct foe exhaust gas to the ETF plant exhaust system.1^1 

22 Flaw Simulation Requirements 

The HDS section test will investigate the effect of velocity, pressure, 
temperature and fuel-air equivalence ratio on foe flame stability performance of foe vee- 
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gutters. Table I outlmes the range of experimental flow conditiorts, defined upstream of 
the flameholders, for the section test. 

Table 1: Experimental Upstream Row Coruiitiims for HDS Section Test 

_ PARAMETER RANGE 

Approach Mach No. 0.18 - 0.48 
Static Pressure (psia) 13 - 8.0 

F/A Equivalence Ratio 0.2 —2.0 

Total Temperature (OR) 530, 930 

The mirumum velocity required to prevent flashback upstream and the 
maximum velocity to produce choked flow at flameholder minimum area bound the 
approach Mach number range. The minimum static pressure in RlAl is limited to 13 
psia, and the maximum static pressure is 8.0 psia because flammability uncertainty 
occurs at low pressures. The equivalence ratio limits are 0.2 and 2.0. The total 
temperature of the primary flow is 530 ®R or 930 °R.l^l 


23 Test Procedures 

The planned test procedures involve defiiting light off and blowoff limits of lean 
and rich hydrogen-air mixtures. Figure 5 is a diagram of the test procedure for lean 
light off and blowoff. 



TIME 


Figure 5: Diagram of Test Procedure for Lean Light Off and Blowoff I 


2-9 






Flow conditions for a given air mass flow rate. mgif.. are set and held constant. Once 
mair established, the hydrogen mass flow rate. mH2/ increases along a linear path. 
The point of light off (point 2 in Figure 5) is defined, and the flame is stabilized (point 3 
in Figure 5) at a higher flow rate for a short time. The point of lean blowoff is defined 
by decreasing mH2 on a linear path until flame extinction ocoirs. 

Figure 6 is a diagram of the test procedure for rich light off ai\d blowoff. 



TIME 

Figure 6: Diagram of Test Prcxedure for Rich light Off and Blowoff 


The hydrogen-air mixture is ignited following the same procedure outlined for lean light 
off. Once the flame is stabilized at the desired fuel-air equivalence ratio (point 3 in 
Figure 6), the point of blowoff (point 5 in Figure 6) is defined by decreasing the pressure, 
P, on a linear path until flame extinction cx:curs. 


3.0 Theoretical Considerations in Model Development 

Before detailing the mcxlel used to predict the performance of the proposed HDS 
hardware, t^»eoretical considerations in the mcxlel development are discussed. Theoretical 
considerations involve: (1) flame stabilization on a bluff bcxly and (2) the well-stirred reactor 
theory. 


3.1 Flame Stabilization nn a Bluff Body 

Inserting a bluff bcxly into a high-speed gas stream creates a region where the 
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flow velocity is less than the burning velocity of the mixture. Figure 7 is a schematic of a 
stabilized flame on a vee-gutter. 


Flame Initiated in 



|4 L 


Figure 7; Flame StabiUzatiiXk on a Vee-Gutter Flameholder 

The flameholder produces in its wake a region of recirculating flow. After ignition of 
the mixture, the redrculation zone acts as a pilot, igniting the fresh mixture in the shear 
layer and anchoring the flame to the flameholder. 

3.1.1 Definition of Flame Stability 

Expression of the flame stabilizing performance of a bluff body usually 
takes ^he form of: (1) the range of equivalence ratios over which combustion is 
achieved, and (2) the maximum flow velocity that a system can tolerate before 
flame extinction.!^ SiiKe the advent of the gas turbine engine, flame 
stabilization on bluff bodies has been studied extensively. 

DeZubay conducted flame stabilization experiments of hydrogen- 

air mixtures on small, ax ymmetrical disks for a range of pressures. DeZubay’s 
study resulted in a correlation between foe fuel-air equivalence ratio, foe 
velocity and pressure, defined at the nuiumum-area. Figure 8 shows DeZubay's 
flame stability parameter, 0, at blowoft as a fuiKtion of equivalence ratio. 

Figure 8 also includes a sixfo-order polynomial fit of DeZubay's flame stability 
parameter. 
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Figure 8: DeZubay's Rame Stability Parameter for H 2 with a Sixth-Order Poljmomial Fit 
For hydrogen-air mixtures, DeZubay developed the following correlation: 


P = U / (pO-61 dO-74) 


where U is the velocity at the disk (ft/sec), p is the static pressure at the disk 
(psia), and d is the diameter of tihe disk (in). For a given velocity, pressure and 
diameter, P is calculated. If P is less than th« blowoff value determined from the 
equivalence ratio, the flame is stable; if P is greater than or equal to the blowoff 
value, flame extinction occurs. Zukoski determined that the blowoff velocity 
for a twcHiimensional flameholder is one-half the blowoff velocity of an 
axisymmetrical flameholder with the same characteristic dimension. 

3.1.2 Theoretical Aspects of Flame Stabilization 

The theoretical study of the mechanism of flaune stabilization on a bluff 
body takes two basic approaches:!^ 
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(1) Following Longweil and associates, the near wake, or 
recirculation zone, of the flameholder is the focus of study. 
Artalysis of die flame stabilizing mechanism involves viewing the 
recirculation zone as a homogeneous chemical reactor. Flame 
blowoff occurs when the residence time of combustion gases in the 
recirculation zone is less than the time necessary to raise the 
temperature of the gases to the ignition temperature. 

(2) Following Zukoski and Marble, the free shear layer 
surrounding the recirculation zone is the focus of study. Fresh 
mixture from the freestieam is entrained in the shear layer and 
ignited by turbulent mixing with previously igruted hot combustion 
products from the redrculaticm zone. The burning mixture in the 
shear layer propagates downstream igniting the neighboring fresh 
mixture in the freestream. At the end of the wake, some of the 
burning mixture is entrained into the recirculation zone while the 
remainder continues downstream. The recirculating flow conveys 
the burning mixture upstreaun where it ignites the fresh mixture in 
the shear layer. This cyclic process is the mechanism of flame 
stabilization. Flame blowoff occurs when the resideiKe time of 
fresh mixture in die shear layer is less than the delay time for 
ignition of the fresh mixture by hot combustion products from the 
recirculation zone. 

On the basis of model simplicity, the author selected die mechanism 
proposed by Longweil to model flame stabilization on the center vee-gutter. 

Also on the basis of model simplicity and historical precedence,f^^' the 
author selected the well-stirred reactor theory to model the recirculation zcwie of 
the flameholder as a homogeneous chemical reactor. 

32 Weil-Stirred Reactor 

The well-stirred reactor, developed by Longweil and Weiss,^^! is a one¬ 
dimensional, mathematicai model for steady, homogeneous combustion. In the well- 
stirred reactor, the rate of mixing is assumed to be infinite; therefore, the combustion 
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gases in the reactor are homogeneous with unitorm temperature. Figure 9 shows a 
conceptual schematic of a well*stirred reactor. 


Premixed, Homogeneous 
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Figure 9: Conceptual Schematic of a Well-Stirred Reactor 

The premixed, homogeneous reactants enter the reactor volume, V, at a 
specified mass flow rate, m, and uniform inlet temperature, Tq. The concentration of 
eadt reactant species entering the reactor is 0 |*. The reactor contains combustion gases 
at a constant pressure, p, and uniform temperature, Tf. Heat, AQ, may be added to or 
extracted from the reactor volume, or the reactor may be adiabatic. For steady 
operation, the combustion products exit the reactor at a mass flow rate, m, and uniform 
temperature, Tf. The concentration of each product species leaving the reactor is 

4.0 Well-Stirred Reactor Model of Proposed HDS Hardware 

The recirculation zone of the center flameholder was modeled as a well-stirred reactor. 

The model assumes the following conditions: 

(1) primary flow from upstream to the flameholder minimum area is steady, 
compressible, invisdd and isentropic; 

(2) the flow experiences no body forces; 

(3) total pressure and total temperature remain unchanged between upstream and 
minimum area conditions; 

(4) the recirculation zone is an adiabatic, constant pressure reactor, 

(5) the hydrogen-air mixture is an ideal gas; 

(6) air is standard composition with argon and carbon dioxide inert, and 
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(7) the hydrogen injected into the wake by the pilot torches is r\eglected. 
Application of the conservation equatioi\s and the chetnical kinetics of hydrogen-air 
combustion results in a set of non-linear, algebraic equations that can be solved numerically 
The author selected a general chemical kinetics program entitled the Lewis General Chemical 
Kinetics and Sensitivity Analysis Code (LSENS) to solve the governing equations of the nnodel. 
LSENS was developed at the NASA Lewis Research Center.l^' LSENS requires the 
following set of input: 

(1) the chemical kinetics of hydrogen-air combustion with Ar and CO 2 mert; 

(2) the fuel-air equivalence ratio of the mixture; 

(3) the mitial static pressure and temperature of the mixture; 

(4) the reactor volume (i.e., the recirculation zone volume), and 

(5) the irutial and maximum mass flow rates (i.e., mass entrained into the wake) 
and a mass flow irKrement. 

The volume of the recirculation zone was determined by correlations for vee-gutter 
flameholders given by Baxter aiHl Lefebvrei^^) The mass entrainment rate into the 
recirculation zone was determined as a fraction of the mass flow rate through the flameholder 
projected area at upstream conditions.I^^l An entrainment rate of 5% was used for the initial 
mass flow rate, and a rate of 20% entraiiunent was used for the maximum mass flow rate. 
LSENS solves the well-stirred reactor equaticms for each mass flow iteration and outputs either 
residence time and outlet conditions (e.g., pressure, temperature, composition, etc) for a 
predicted stable flame or an 'apparent blowout' message for a no-flame condition. 

5.0 Project Continuation 

HDS technology development testing in RlAl is in progress at AEDC. The author 
began ruruting the LSENS software for the expected test conditions in RlAl on August 11,1993. 
The flame stability predictions by LSENS is being generated and analyzed. The flame stability 
limits predicted by LSENS will be compared to DeZubay's database and the HDS technology 
development database when testing is complete. Results of the comparison of the analytical 
predictions and the experimental data will be reported when the analysis is concluded. The 
project will continue through December 1993. 
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SIMULATION TECHNIQUES FOR PRESSURE AND TEMPERATURE DISTORTION 


Anthony L. DiPietro Jr. 

Pb D Graduate Student 
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Abstract 

This is a presentation of simulation techniques utilized to generate pressure and temperature distoctioa in 
ground gas turbine engine test facilities. It reveals what is currently available for generating pressure and 
temperature distortion, ramps, and transients. It also reveals some operational problems, advantages, and 
disadvantages associated with each simulator. It menuons that as a result irom 'ms vonc. aew concepts were 
conceived for producing pressure and temperature distottioo. These concepts were developed to possibly improve 
upoc what was already available in the Arnold Engineering Development Center gas turbine engine test cells and 
other ground gas nirbine engine test facilities. It also mentions that this work is the result of a continuing effort at 
Virginia Polytechnic Institute and State University to investigate axial flow fan and compressor response to 
unsteady pressure and temperature distortioa and combined pressure and temperature distortion. 
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SIMULATION TECHNIQUES FOR PRESSURE AND TEMPERATURE DISTORTION 

Anthony L. DiPietro Jr. 

Introduction 

Throughout an aircraft's flight envelope, the inlet and engine must remain aerodynamically 
compatible. Situations which lead to cominessor stall, engine surge, or other malfunctioning of the 
propulsion system must be avoided or at least reduced to a tolerably low frequency of occurrence. Such 
situations are produced by a departure of airflow, as delivered firom the inlet to the engine, from the ideal 
of a flow that is uniform in gaseous content, temperature, pressure, and uniformly axial in direction. No 
real flow ever achieves this ideal however, in the pursuit of high inlet pressure recovery for on design, 
flows are sufflciently close to it so as to pose no compatibility problem.s. 

.Aircratt engine manufacturers are concerned with engine operation under non inlet flow conditions. 
Specifically, the interest is in engine perfonnance with distorted inlet temperature and pressure profiles. 
Inlet distortion includes pressure distortion associated with aircraft pitch, roll, and yaw rates and boundary 
layer separation from the fuselage and inlet duct Temperature distortioa is produced by the ingestioa of 
exhaust gases from armament firings such as guns, rockets, and missiles. Also, gas turbine engine exhaust 
gas re ingestioo that occurs in helicopter and VTOL/STOL aircraft operation. Aircraft engine 
manutacturers as well as engine test facilities are makmg an etfoit to develop devices that can simulate 
these pressure and temperature distortions for direct connect tests on gas turbine engines. These devices 
produce actual flow conditioos at the engine inlet that may exist during an aircraft maneuver or armament 
firing. Data regarding engine perfonnance and stability under distorted inlet flow conditions, generated by 
these devices, can now be obtained. This information will be useful in making remedial measures to the 
propulsion system or to redesign the propulsion system to tolerate the distoruon. 
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Methcxtology 


The oietbod of generating a pressure distortion at the face of a gas turbine engine inkt is to 
produce a physical blocicage in a sector of the inlet Tbis physical biocicage locally reduces the axial 
velocity of the airstieain directly downstream from the blockage. This reduction in axial velocity translates 
mto a local reduction in total pressure behind the blockage. This physical blockage may be a screen, a 
mechanical vane, or some other physical obstruction. Another method of generating pressure distortion is 
to mount a series of aiijets in the inlet duct near the engine face. These jets are mounted such that the 
nozzle exits of all the jets ore in the same crossectional plane and point directly into the inlet flow. These 
mnets lischarae high pressure v;omnres.sed lir ;nto rhc incoming inlet 'low at high velocities. A 
momentum exchange occurs between the locally injected high velocity compressed air from the aiijet 
nozzles and the oncoming inlet flow to produce a local reducdoa in fluid axial velocity. This reduction in 
fluid axial velocity translates into a reduction in total pressure. The method of generating a temperanire 
distortion at the engine face i > to mount some type of heater or burner in the inlet dua . Tbis beater ' 
burner may be subdivided into sectors that operate independently of each other to locally increase v 
temperature in sectors of the inlet duct Another method of generating temperature distortion at an engine 
face is to actually mount and burn a missile or rocket motor in front of the inlet It is desirable that the 
pressure and temperature distortion devices have the ability to change patterns- while a test is tunning. 
Tbis ability allows transients to be simulated at the engine inlet that may exist during an aircraft 
maneuver or armament firing. 

Results 


The Naval Air Propulsion Test Center (NAPTQ utilized a distortion valve to produce steady 
state and transient pressure distortion. The distortion valve consisted of a square frame with a series of 
parallel flat vanes. Tbese vanes were connected to operate mechanically as a "Venetian blind” to block 
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quadrants of the inlet The valve was subdivided into four independently controlled square quadrants. The 
distortion valve produced good spatial pressure distortion however, it could not produce complex 
distortion patterns. Another disadvantage to the distortion valve is that the crossectioa of the inlet had to 
change from square to circular for direct connect tests. This transition would alter the pattern as it left the 
distortion valve. 

The National Aeronautics and Space Administration (NASA) Lewis Research Center utilized an 
aiijet distortion generator to produce transient pressure distortion, it consisted of 54 aiijets mounted in the 
inlet duct of an engine test cell. Six main control valves mounted outside the inlet duct controlled the 
discharge of compressed air to the jets. Each valve was installed at 60 degree intervals around the inlet 
duct and supplied air for nine jets. The generator operated by producing a momenmm e.tchange between 
the injected compressed air and the inlet duct airflow to obtain a reduction in total pressure. The device 
produced good controlled spatial pressure distortion and could produce o-ansients during a test However, 
it had a coniplex valve, air transport, air conditioning, and control system. During tests, the system 
experienced problems with the control valves. NASA also utilized a rotating screen assembly to generate 
steady state pressure distortion. It consisted of a mount that could support precut screens of varying 
solidity. This mount could rotate ± 30 degrees and was anchored in the inlet duct of an engine test cell. 
The device produced good spatial pressure distortion and was relatively simple and inexpensive. A 
disadvantage with the screen assembly was it could not produce transient while a test was rutming. It was 
limiting to the types of patterns that could be produced. It also bad a problem with the screens breaking 
apart during operating conditions with high mass flow. 

The Arnold Engineering Development Center (AEDQ also developed an airjet distordoo 
generator to generate pressure distortion for engine tests. It consisted of 56 au- injector ports that were all 
independently controlled and mounted in radial, aerodynamically. >ireamiined. supporung struts. This 
strut assembly was anchored in the inlet duct of an engine test cell with the airjet nozzle exits pointing 
into the oncoming inlet airflow. Utilizing a momentum exchange between the primary air in the inlet duct 
and secondary air injected by the nozzles, the generator produced controlled spatial pressure distortion 
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and pressure transients during a test. It bad (be abiUty to produce m(»e complex patterns than ibe NASA 
aiijet because each injector port was independently controlled. It bad a complex valve and control system 
and possessed some of tbe same disadvantages as tbe NASA airjet However, it is speculated that tbe 
AEOC airjet generates less freestream turbulence than tbe NASA airjet because tbe aiiyets are mounted in 
aerodynamicaUy streamlined struts. 

A pressure distortion device developed by General Dectric in Cincinnau Ohio was ibeir P3G 
wbicb is a planar pressure pulse discrete frequency generator. It consisted of a rotor driven by an electric 
motor between two stators. It was a relatively simple device and produce a good planar pressure pulse. 
Tbe frequency of tbe pulse was controlled by tbe rotor speed and the amplitude oi ihe pulse was controlled 
by tbe spacing between rotor and the stators. It was limiting because it could only nnxluce one narticular 
type of distortion pattern, a planar pressure pulse, and could not produce transients. 

The NAPTC had two techniques for producing temperature distortion. One was a bmkr plate 
rocket motor gas generator mounted direcdy in front of tbe inlet duct of an engine. It was designed to bum 
specially fabricated grains of solid rocket propellant. Tbe other simulator utilized a cylinder that burned 
soUd propellant inside two separate powder chambers outside the inlet duo. A circular manifold 
channeled tbe combustion gases around tbe inlet dua and injected it into tbe inlet ftow through the inlet 
dua wall. Both simulaaon devices produced correa gas chemical composidon and temperature ramps. 
Peak temperature ramps up to 10,000 F/s were adiieved. An advantage to tbe rocket motor that burned 
propellant outside the inlet duct and cbanneled tbe combustiOQ gases into tbe inlet was that it bad two 
powder chambers. These chambers could be ignited at dififerent times to simulate multiple rocket firings. 
Both simulation devices had the disadvantage of being unable to produce controlled spatial temperaoire 
distortion. Other disadvantages are that tbe powder chambers bad to be retilled alter each tiring and that 
conosive agents HCL and A1203 solid parucuiaies were produced and attached to engine internal 
component. 

Tbe NASA Lewis Research Center designed, fabricated, and operated two gaseous fuel hydrogen 
burners to simulate temperature distortion. One thermal distortion generator was designed for lame 
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engines sucb as tuibojets and turbofans. The other was a smaller tbennal distortion generator for use on 
turboshaft engines. The burner used for the turbojet and turbofan engines was a circular, four quadrant, 
segmented, hydrogen fueled burner. It utilized fuel injectors and V-gutter tlamebolders with swirl can 
ignitors. Each quadrant was independendy controlled to produce steady state or transient temperature 
distoruon. The turboshaft engine thermal distortion generator was circular and subdivided into eight 4S 
degree sectors. It had three swtri cup combustors in each sector for a total of 24 swirl cups. Both burners 
produced good spatial temperature distortioo and temperature ramps. Temperature ramps up to 17,750 R/s 
were achieved with the large distortioa generator. Both burners produced inconrea gas chemical 
composiuon due to the hydrogen and air combustion reaction. Hydrogen was used as the fuel because of 
us low tljunmability limit, fast ignition, high temperature rise, and steep temperature ramps. A problem 
with both designs was uniform ignition in each sector/quadranL Experimental tests on a turbofan engine 
indicated that the hottest most rapidly changing temperatures occurred at the center of the burner and 
produced core compressor stalls but not fan stalls. 

AEDC also developed two thennal disttmion simulation techniques. One simulator was a 
modification to their airjet distortion generate. Instead of injecting ambient temperature compressed air, 
the air was heated. The injection of high temperature compressed air produced the thermal distortion. 
This technique was easily accomplished with minor hardware modifications to the aiijet system. It 
produced good spatial temperature disttmioo however, pressure distortion accompanied the temperature 
distortion. It had the advantage of not having to stabilize a combustion reaction inside the inlet duct where 
the high velocity air was flowing. Some disadvantages were the aiijet produced low temperature ramps 
and required large volumes of hot air to achieve a desired temperature rise. The other temperature 
distortion sunulauon device at AEDC was a tumsule mount that supported four small solid rocket motors 
dll at 90 ilegrees to each other. This tumsule mount bad the capability of being posiuoned .it several 
distances from the inlet. This simulator produced good temperature ramps and correct gas chemical 
composition. It did have some of the same disadvantages as the rocket gas generators used at the NAPTC. 
Specifically, it produced corrosive agents HCL and A1203 solid particulates that attached to and damaged 
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engine internal components. 

An effort to investigate tan and compressor response to unsteady pressure and tenqierature 
distortion and combined pressure and temperature distortion was begun at Virginia Polytechnic Institute 
and State University in their Aerospace and Mechanical Engineering departments. The effon initiated 
with the design, fabrication, and operatioo of a gaseous fuel burner to simulate unsteady temperature 
distortions associated with gas turbine or rocket engine eiihaust gas mgestion into axial flow fans and 
compressors. The result was a cylindrical four quadrant propane burner that led to a Masters thesis. The 
device consisied of an 18 inch diameter cylinder, three feet in length, that was subdivided into four 
quadrants by use of steel plates. Occupymg each quadrant was a cenierbody that contained a fuel injector 
and ignition plug. The device wa.s specifically designed to produce temperature transients and un.siead\ 
temperature distortion with multiple quadrant combination patterns. However, it did have the ability' to 
produce steady state spatial temperature distortion. The device produced good spatial temperature 
distortions however the temperature ramps produced were low. It was speculated that if gaseous hydrogen 
was used as the fuel, steeper ramps would be generated. An effort is currently underway to develop a 
similar device that will generate an unsteady pressure distortion and the goal is to use the two devices in 
combination on a multi-stage transonic fan. This work wiU lead to a Doctoral dissertation at VPI. 

ConclusioBS 

Two techniques exist for producing a temperature distortion at the face of a gas turbine engine. 
One technique utilizes an actual rocket mottx’ that is mounted and burned at the engine inlet The other 
techmque involves the mounung of a burner or heater in the inlet duct that is capable of locally raising the 
lemperature ui sectors of the inlet Several devices exist for producmg a pressure distoraon at the lace of a 
gas turbine engine. One class of these devices operates by producing a physical obstruction in a sector of 
the inlet The velocity downstream from this obstruction is reduced and this produces a total pressure drop 
behind the obstruction. The other class of pressure distortion devices operates by utilizing an axial 
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momentum exchange between the pnmary miet duct auflow and secondary air that is injected directly in 
the counter direction of the primary airflow. Again, a veloaty reducuon in a sector of the inlet duct is 
obtained by a reduction in airflow axial momennim of the primary airstream in that particular sector. This 
translates into a total pressure drop in that particular seaor. 

.A large amount of published data, engine test data, demonstrates that pure temperature ramps at 
the engine face with rates on the order of 2000-10,000 R/s can induce engine surge or combustor blow 
out. Tests indicated that the m^inide of the temperature ramp had a mote pronounced effect than either 
emulated Ntendy state spatial temperature distortion or the absolute level of temperature obtained during 
(he transient. 

A plot of inlet total pressure amplitude verses pressure pulse duration data obtained by the 
NASA Lewis Research Center aiijet system indicated that the longer the pulse, the lower the total 
pressure. StaU occurred during long pulses only which indicates that stall is likely to occur over low 
frequency ranges. A possRile explanation is that there is a fmite time interval required for a rotor blade to 
stall in response to a total pressure distortion. During shorter durahoos. or higher frequencies, the 
distortion is too brief for the rotor blade to stalL Results indicate that the stall tolerance is a function of not 
only the instantaneous distortion level but also the rate of change of inlet pressure and dwelling time of 
the fan or compressor rotor blading in an engine inlet distortion. 

This work indicates that there are several devices available for producing pressure and 
temperature distortion. The review also revealed that each device has certain operational problems, 
advantages, and disadvantages associated with its design. The type of distortion research and experiments 
that will be conducted will dictate which simulation device is appropriate for that particular applic?*'on. 
This work is a conunuauon of the etfon being conducted at Virginia Polytechnic Insutute anu i:ate 
University to invesdgate fan and compressOT response to unsteady pressure and temperature distortion and 
combined pressure and temperature distortion. This work is also a contribution to the effort currently 
underway at the Arnold Engineering Development Center to establish the capability of pressure and 





teoiperature distortion tesung to the AEDC gas turbine engine test cells. Resulting from tins work came 
new concepts for generating pressure and temperature distortion that may be better tailored for the AEDC 
gas turbine engine test cells and VPI fan and compressor rig research. Some of these concepts may be an 
unprovement of what is already available at AEDC and other ground gas turbine engine test facilities. 
These concepts may be investigated further at AEDC and VPI. 
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Design ot an Aluminum Combusuon Simulator 

Robert]. Geierman 
Graduate Research Assistant 
Mechanical Engineenng 
University of Tennessee Space Insutute 

Abstract 

This research presents the design of a combusuon device which will simulate the aspects of aluminum 
panicle combustion inside a solid rocket motor. At the same time this combustor will allow optical data to 
be collected during the combustion process. 

The first step in this combustor design was to examine past expenments and methods which have been 
used to examine the combustion of aluminum particles. Next, the fuels and oxidizers, which will be used in 
the first planned combusuon expenments, were selected. The flame temperatures and uanspon properties 
of these fuels and oxidizers were calculated using a version of a NASA chemical equilibrium code. The 
third step in this combustor design was to define the physical dimensions and maximum operating times 
which would be allowed. These parameters were defined through simplified stress analysis and heat 
uansfer calculations, respectively. The final steps in this combustor design process were to define the 
procedures and instrumentation to be used during the operation of this combustor. 


4-2 






Design ot an Aluminum Combustion Simulator 


Robcn J. Geierman 


InffQduciiQn 

For many years, metals have been used as additives in solid rocket motor propellants. The reasons 
for the use of these additives include their ability to help eliminate combustion insDbilities and their ability 
to increase specific impulse. By far, the most popular metal additive is aluminum. The popularity of 
aluminum usage is due to its high heat of formation (of its oxides) and its relatively low cost. 

Although aluminum is widely used as a propellant additive, much of its combustion characteristics 
and opucal properties, especially during solid propellant burning, have not been quantified. By obtaining a 
better understanding of aluminum combustion, one could possibly increase the specific impulse of solid 
rocket motors without degrading combustion stability. The increase in specific impulse could be 
accomplished by either increasing the aluminum content to near its stoichiometric limit, while maintaining 
a high combustion efficiency, or by minimizing the two-phase flow losses. Similarly, an adequate 
understanding of the optical properties of aluminum is important to predict the radiation heat transfer to the 
propellant surface and the spacecraft itself, from the exhaust gases. In order to gain a better understanding 
of the combustion of aluminum and its properties, it becomes necessary to observe the process in the 
environment which is found inside a solid rocket motor. Since it is nearly impossible to observe the 
phenomena inside an actual solid rocket motor, a device is necessary which closely simulates these 
conditions while at the same time allows for observations of exhaust or combustion process to be made. 
This report contains an overview of the design process for a combustor which will simulate the combustion 
of aluminum under solid rocket motor conditions. For a more detailed report see (1 ]. 
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Methodology 

In order to design this combustor, it was necessary to first examine the past experiments which 
have investigated the aluminum combustion phenomenon. From an extensive examination of many reports, 
it was found that there is still a lack of information dealing with the aluminum :ombustion process. Never 
the less, the previous experiments were divided into three groups. These groups consisted of burning of 
propellant samples, particle/flame burners, and exotic methods. Burning propellant samples was 
accomplished by simply igniting propellant samples in some sort of combustion bomb and observing the 
process (mainly through the use of high speed cameras). Particle/flame burners merely consisted of 
injecting aluminum particles into a flame or torch. The drawback of these experiments was due to the 
combustion process in the flames occurred at atmospheric pressure rather than the high pressures found in 
typical solid rocket motors. The exotic methods consisted of the use of shock tubes, ignition of large 
particles with lasers, and electrical ignition of aluminum wires. These exotic methods usually burned 
aluminum in shapes and sizes which were not typically found in solid rocket motors. 

All of these experiments observed the combustion of aluminum in one way or another. However, 
these experiments often failed to bum the aluminum under conditions which are found in solid rocket 
motors. Additionally, these experiments rarely collected any optical data from the burning aluminum and 
aluminum oxide particles. For these reasons, it was decided to design a combustor which would allow for 
observations of aluminum combustion under a simulated solid rocket motor environment. This combustor 
should also be versatile enough to observe the effects of individual species on the aluminum combustion 
process. This combustor concept should be inexpensive to build and operate. For these reasons, it was 
decided to start with a combustor design which would be as simple as possible. The original combustor 
should start as simple as possible and gradually, as more knowledge is gained, work towards more complex 
experiments. Keeping this in mind, the original combustor was designed to use a coaxial injector with the 
aluminum particles injected with the gaseous fuels (see Figure 1). This original combustor concept was 
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later modified. The reasons for these modifications are discussed in the conclusions and recommendations 
section. 


Since the goal of the onginal combustor concept was to be as simple as possible, it was decided to 
use both gaseous fuels and oxidizers in the original design (and small mass fractions of aluminum panicles 
(0.1 to 5 percent)). Because hydrogen and water vapor are prevalent products of combustion in solid rocket 
motors, it was decided to use gaseous hydrogen as the fuel to be used in the first combustion experiments. 

It was also decided to use methane as the second fuel in this combustor. In this way the effects of carbon 
monoxide and carbon dioxide can be studied along with the effects of hydrogen and water vapor. 
Eventually, this combustor will be capable of simulating the majority of species produced in a solid rocket 
motor. 

Before the specific design characteristics of this combustor could be determined, it was necessary 
to perform a combustion analysis of the fuels to be used. This combustion analysis was performed by using 
the NASA equilibrium code (SP-273 [2]) to determine the chemical species which will be present in this 
combustor. It was found from these species that a significant mole fraction of aluminum hydroxide formed 
during the combustion of low mass frar-Jons of aluminum. As the aluminum mass fraction was increased 
to beyond 1 percent of the total fuel and oxidizers, the mole fraction of aluminum oxide produced was 
greater than the mole fraction of aluminum hydroxide. It was therefore recommended that the combustor 
be operated at aluminum mass fractions greater than 1 percent of the total flow. In addition to the 
chemical species present, this code also predicted the flame temperature and transport properties of the 
combusuon products. These properties were then used to predict the heat transfer rates and physical size of 
this combustor. 

After the combustion characteristics were calculated, it was possible to determine the physical 
dimensions of the combustor. Since this initial combustor design should be as inexpensive as possible, it 
was deemed necessary to make this combustor as small as possible. This not only makes the design less 
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costly (due to the use of less material), but it also makes for a safer design. A smaller volume at high 
pressure should cause less damage, than a large volume, in the event of failure. In any event, by starting 
with a small inexpensive combustor concept, it should be possible to scale to larger combustors if desired 
later m this research program. 

Since this initial combustor will be as small as possible, it was felt that a central fuel injector with 
a 0.25 in. outside diameter would be sufficient for the first two fuels to be used. By neglecting the two- 
phase flow effects of the aluminum particles (due to their low mass fraction of the total flow), it was 
possible to use the 0/F ratios (corresponding to the maximum flow rates) to size the other injector pons. 
Again, to minimize the cost of this combustor, it was decided to use standard size stainless steel tubing for 
the injectors. The lengths of these injector tubes were designed to create a constant velocity at any radial 
location. Also, the thickness of the injection tubes was found using a simplified pressure vessel equation 
and a factor of safety of 2. The dimensions and layout of the injector face are shown in Figures 2 and 3, 
respectively. 


Injection Tube 

Inside Diameter 

(in.) 

Outside Diameter (in.) 

Tube Length (in.) 

Fuel 

0.238 

0.25 

4.5 

Oxygen 

0.710 

0.75 

4.75 

Niuoeen 

0.930 

1.0 

4.75 


Figure 2: Injector Dimensions 

I 
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Igniter Wires 


Flow 

Recife uhtors 


Metal to Glass Seal 



Figure 3: Injector Face Layout 
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The next step in this combustor design was to determine the dimensions of the nozzle. In order to 
minimize the cost of this combustor, the original combustor concept was designed with a variable area 
throat (see Figure 4). The nozzle was to be cooled with water and the throat area was chosen to provide for 
choked flow at the maximum expected flow rates. Heal transfer calcubtions were performed to determine 
the necessary flow rates for the coolant water. The dimensions of the nozzle were determined with the 
same pressure vessel equation used for the injector calculations. However, as the temperature increases, the 
strength of the (standard size stainless steel tubing) nozzle decreases dramatically. Therefore, the required 
thickness of the nozzle walls was far greater than the injector wall thickness. The dimensions of the nozzle 
are listed in Figure 5. 

Unlike the nozzle and injector dimensions, the length of the quartz tube is determined by the flame 
dimensions. Since the objective of this combustor is to observe burning aluminum panicles, it is necessary 
that the length of the tube be sufficient for complete combustion of the particles to occur. Therefore, the 
maximum ignition and bum times were calculated for 25 jim oarticles. Using these times and the 
maximum fuel injection velocity, it was possible to determine the length required for the panicles to ignite 
and bum. An additional length of quaru tubing was required for the inviscid core which occurs at the 
injector face (see Figure 6). With an added margin of error, the total quartz lube was determined to be 10 
in. long. Since the inside diameter c nz tube must be the same as the inside diameter of the 

nitrogen injector lube (to eliminate recirculation zones near the quartz tube wall), an inside diameter of 0.93 
in. was selected for the quaru tube. However, due to the high heat transfer rates and the poor strength of 
quaru in tension, it was necessary to design a two wall construction for the quaru combustion chamber 
(see Figure 7). The reasons for this two wall construction are discussed in the conclusions and 
recommendauons section. 
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Pigiire 4: Nozzle Design Conrigiiration 













igiire S: Nozzle Dimensions 
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Figure 7: iiijeclur Munifulrt Design 







Conclusions and Recom mendations 

As stated earlier, the design of this combustor is intended to start as simple ei possible a: d 
gradually work towards more complex designs that would more accurately simubie the combustion of 
aluminum in a solid rocket motor environment. This design work was mainly based upon simplified 
pressure vessel equations, equilibrium chemistry calculations, and a simplified heat transfer analysis. (Each 
of these design enterion are discussed below.) 

Although the simplified pressure vessel equations are by no means a finite element analysis, they 
should be sufficiently accurate for the intended operations of this combustor. However, the dimensions of 
the central injector were arbitrarily chosen. In reality, the area and flow rate of each injection tube could 
drastically change the stability of the flame and the combustion products which are formed. In addition, 
this design assumed that the nitrogen and oxygen would not intermix. In reality, this is not the case. This 
is due to the fact that there will be some vonicity produced between the oxygen and nitrogen interface. 

Any substantial mixing between these two gases could cause the flame to expand to the quartz tube wail 
resulting in excessive heat transfer rates. For these reasons, it is imperative that the flow rates and flame 
stability be verified and all necessary redesign work be completed before the initial test operations begin. 
These redesigns may include changing the dimensions of the injectors, the flow re-circulator locations, and 
the flow rates of the gases. 

The use of equilibrium chemistry for this combustor should provide a fairly accurate overview of 
the species, transport properues, and temperatures which will be produced in this combustor. The accuracy 
of the NASA equilibrium chemistry code has been confirmed in numerous tests. However, since these 
chemistry calculations assume infinite reaction rates and reactions which continue until completion, the 
actual species produced in this combustor will probably be slightly different from those which were 
predicted. In reality, the rate of reaction will be mainly controlled by the turbulent mixing of the species at 
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ihe injector face. Also, the reaction of the aluminum panicles (with oxygen) may be difficult due to the 
(possible) lack of oxygen along the center of the combustor. Further complications which may affect the 
aluminum particle combustion include the formation of a protective oxide on the panicles and the collision 
of molten aluminum particles which may form large agglomerates (and therefore, cause substanually 
longer bum times). All of these factors must also be experimentally verified and any necessary design 
changes must be incorporated into the combustor before actual tests con begin. 

Another problem which may occur during the operation of this combustor is the erosion of the 
nozzle throat. Although the heat transfer calculations should insure that the nozzle will be maintained 
below Its melung temperature, it does not necessarily guarantee that the nozzle throat will not erode. The 
high fiow velociues of the reactive and corrosive combustion gases can still erode (and possibly react) with 
the stainless steel nozzle. Since the velocities should be their maximum at the nozzle throat, it is 
recommended that the nozzle design be incorporated with graphite inserts in the nozzle throat. Although 
the graphite insert will still probably erode, it should not threaten the integrity of the nozzle itself. This 
should, in the long run, decrease the cost and increase the safety of the nozzle design. Additionally, since 
the ceramic needle may cause the nozzle to plug with aluminum slag (resulting in catastrophic over 
pressures), it is recommended that the variable throat area design be discarded in favor of several different 
nozzles. Each of these nozzles would have graphite inserts with the throat area corresponding to choked 
flow for the flow rate of the particular experiment. Although this may be more expensive, it is necessary to 
insure the safe operation of this device. 

Due to the high heat transfer rates it was necessary to use a two wall construction for the quartz 
observation windows (Figure 4 and Figure 6). This was also necessary due to the poor strength of quartz in 
tension (especially at high temperatures). If a single glass tube were used for this combustor, it would 
require a tube which would be approximately 0.5 in. thick for a 1.0 in. inside diameter tube. This thick tube 
would act as a (powerful) lens which would make data acquisition difficult. Alternatively, the combustor 



could use a single chamber with a thick recessed quartz window. In order to cool ihis recessed window, it 
would be necessary to purge this window with a coolant gas. This may make it difficult to observe the 
combusuon process along the enure length of the combustor. Also, this altemauve design may result in the 
formation of recirculation zones of hot combusuon gases near the combustor walls (unlike the two wall 
design) resulting in cauisaophic heat transfer rates. Not only may this altemauve design be more 
dangerous, but it could also be more expensive to develop the purge system for the recessed window 
design. For these reasons, it is recommended to develop the two wall design configuration. 

Obviously, the combustion of aluminum particles in a flame may not accurately simulate the 
heterogeneous reacuons which occur within an actual solid rocket motor. However, this combustor concept 
is designed to be ca:?able of observing the effects of the individual species on the aluminum combustion 
process. This design should also be able to observe the effects of individual species on the opucal 
propemes of aluminum and its oxides. In order to accurately simulate the solid rocket motor species, it is 
necessary to operate this combustor in the fuel rich regime for the first two suggested fuels (hydrogen and 
methane). (The two main species produced in a typical rocket motor are free hydrogen and carbon 
monoxide.) Fuels which will be used in later experiments should produce hydrochloric acid and free 
nitrogen (two other species produced in solid rocket motors). It is recommended that each of these species 
and their effects on aluminum be studied individually (if at ail possible). Eventually, after numerous 
limited specie tests, this combustor will be able to observe aluminum combustion in an environment which 
contains all of the prevalent species found in a solid rocket motor. 

Once this combustor is built (after all of the necessary redesigns have been incorporated), it should 
provide a versatile research tool. Not only should this combustor provide insight into the aluminum 
/mbustion process, but it could also be a valuable tool for the prediction of coaxial jet mixing and 
combustion. This combustor may even be able to suggest new fuels and additives for solid rocket motors. 

In conclusion, this combustor (or its derivatives) should be a valuable research tool for years to come. 
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Condiuons m Solid Rocket MotorsA Thesis Presented for the Master of Science Degree, August, 1993. 

2. Gordon and McBride, "Computer Program for Calcubtion of Complex Equilibrium Compositions, 
Rocket Performance, Incident and Reflected Shocks and Chapman - Jouquet Detonations", NASA SP-273, 
February, 1973. 
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VALIDATION OF A RECENTLY DEVELOPED SOLID PARTICULATE DISPERSION 

MODEL 

William Kvasitak 

D(>partiii(‘iit uf Mechanical and Aeronautical Engineering 
Clarkson University. Potsdam. NY 13699 


Al>stracf 

X mIkI.ii ii>u III iraiispori iiuxld for prcdiffmg, solid particulate concentrations is presented. 

( Iiiilidi'iiri- ID lilt- model i.s supporled liy a coinpari.son of numerical computations from tlie L'l’MSD code 
il.-v,.lolled at .-M IX' with analytical ami experimental data. Kxact solutions of the etpiations of motion 
ire compareil with numerical results in order to validate Brownian motion ami gravitational sedimentation. 
.\ comparison with experimental <lala for particle deposition in a channel is u.sed to validate the turbulent 
dispersion model The instantaneous fluctuations are simulated as continuous (iaussian random vector fields. 
Ml** l-agraiigian e<(nafioiis of motion provirle particle velocities induced by Stokes drag. Brownian dilTiision, 
•iml uravii v An a.ssumptioii of one way fluid particle interaction and no particle - particle interaction places 
a maximum mi the particle concent rat-on of less than U.U2VJ by vohime. Three validation test cases are 
l•r•■s)'nl>■d for <<imparisoii wiili analytical and ••xperimeiital results. The e.xact solutions and data are in 
'■xecllent agreement wiili the nnmerical results. 
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Iiitroiluc-tioii 


rciilniiic.il prolili'iiis in(lnrt‘(l l>y polliilioii liave n^rpiiily heroine a priority to the Air Force, and 
inankiinl in u* ll•Tal Ddficiilties ransed liy and rain. snioR, oil spills, and forest fires, have caused the 
I MMiniiini iiial I’loti i tion Ai;enry (1'1‘A) to place strict regnlalioiis on the total amount of emissions allowed 
'\iilim (ontanuneiit areas. .M-’nt' has always operated within EPA (guidelines, but as EPA regulations 
l eiiinie more ainl more strict the need for atmospheric plume dispersion predictions become imperative to 
the \ED(' lest laciliiy permit process I'luler the 1992 (Jratlnate Student Research Program the LTM3D 
'■>lid particulate dispersion code was developeil to help fulfill this need. The motivation and methodology 
n-ed in ilie l.TM.ID model are detailed in Kvasnak (1]. This paper details the next step in the develojjment. 

I ■■ . the \alidation of the ro<le for simple test probleias. 

Hecanse o( I he importance of soli<l particle (|is|)ersion in numerous industrial and environmental procesties. 
lunch rt'search has been done in the area of aero.sol |>arlicle dynamics. (’on.se<iuently, a wealth of knowledge 
• xist- on the snl>|eci Euchs [2). lli<ly ami Ikock (3j. van de Hulst [4]. Twomey [.')]. and Cadle [ti] are just a few 
■ 4 I he classic text hooks on I he subject. Part icle dispersion in turbulent flows was simulated by Ahmadi [7], Li 
Old .Ahmadi [.s], and ()nms I'l al [9] The overall goal of t his project is to provide a large-scale simulation of 
itmosphenc dispersion of fine particles I he u.se of a supercomputer facilitates large-scale simulations which 
may provide the nee<led insight into the processes that dominate the motion and dispersion of fine particles 
III tli.> atmosphere The report of Kvasnak [1] detailed the development of an efficient algorithm for this 
pnrpo.se: however, the ( rank - Nicholson discretization has been replaced by a fully implicit backward-Euler 
scheme 1 his report details the validatioti and improvement of the LTM3D rode in order to provide a useful 
\alidated ami calibrated |>rodiict for AEDC 


Roviow of the LTM3D Motliurlology 

Kvasnak [1] .lescribes in detail the met liorlologx behind the LTM3D code. Here only a brief review is 
piesfiited III the sniiiilai ion a steady-state flow field is taken from an incompressible Navier-Slokes equation 
-olv.r I he bauraiiRian equation of motion for a small particle is then solved in time for a release of particles. 
I hi- concent rat loii at every location in the flow is taken as a mass percent of particles in a given region, also 
died a C. II. I he distribution of particles is calculated on a local and global level. This method provides a 
Miin’-amirai*’ ioi» cf tiu* fiispersioii of an 

Morhaiiirs of A^'ro^olft 


M,\ d. Iimii..n. an aei.,sol is a siispi-nsion of soIkI or liquid particles in a gas. Common aerosols found 
III n.iinre iiichide dust, smoke, log. haze, and smog .Aerosols typically range in size from U.OOl to 100 
microns m diameter fo place these dimensions in perspective note that the mean free path of air is about 
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II 1)7 Miicroii'i ;iiiil iliai I In' li of visiltU* is lielwppii 1) 4 to U.7 microns. Also, a diameter of 0.01 

mirimis ('(nrpsiiDiid-i in iln- Iraii.siiton limii Itciween molecule ami particle. Particles larger tliaii 

lull mn’i'iiii'; tvpn allv •!'< no) remain sn.'^piMided in air for a significant time iluration. Tims, the behavior of 
p.iil nles is signilicant ly ali'ecled by Iheir size. 


.Aofosol Partielo Motion 

file l,.ini anniaii e(|n.iiK)n ol mol ion for a Stokesian aerosol particle can he written as 

till'' dT/n/ , 
m—— = —-—(iP - ii' ) + im/. 

(It ( ,. 

ivlieii- III IS the ni.iss of the particle and ( is the ('iiimingham correction factor given by 

■2A 


( 1 ) 


|)|\ldnin f.'l ill by . yield.s 


l/l/) 


U here 


for relatively larne parii'des, (a: I ami 


/ 1 l'/\' 

(‘2) 

- + tl'' = 11 ■ -)- TlJ 

(3) 

defineil by 


inC^ il-p>'('. 

dT//f/ 18/1 

(4) 

jrr/'* 

in = ——// . 

(> 

(5) 

_ r/V 

“ 18p 

(fi) 


Siiijiilatioii of a Turbulent Flow Fiebl 


1 he I nrbnieiit velocii y field is constructed by the superposition of a fluctuating velocity onto a steady 
velofiiy field In K\a.snak [Ij a zero correlation (I'aussian random variable with a strength «' was generated. 
I he Instatit aiieoiis velocity at the positioti of the particle was given as 

1/ = r 4- I/' (7) 

where / is obtaiiieil from the Navier-Stokes solver Then, a tdaussian random variable with strength equal 
I't the rtiis velocity was used as a fliictttating cotnponent. 


1 his procedure essentially models turbulence as a white-uoi.se process. Kraichnan (10] suggested a sim¬ 
ple tiu'thod ol geiieratitig a (laussian ratidotn fiehl which resembles isotropic ttirbulence. Accordingly, the 
iii‘'tant.iiienns lbi< tiiatinn vi locitv is ni\ee as 


- y Y ( ^ -c* + ^2 "Tl^iilsinllt:,, ■ -f-Wnt") j . 

V., = i / 


( 8 ) 
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Ill lilt's t“i|ii;<i ion. iho iinaiititips if) anil ir< arp ilpfinecl a.« 

»i(^T, HC, (II) 

i73(t„) = . il„, (10) 

n ll li 

■ lll(^*n) ' ■ ll'j(A.‘||) = 1). (ll) 

111 .iisiiri' iii'oinpii'.>ssil)ility. 1 he (|iiantities y, aiul ^ are (lans-siaii distribuled random variables witli a 
--laiiilaril i|>'\iaiioii of 1. l.ai'li roiiipoiieni of h is a (laussiaii random varial)le willi a standard deviation of 
1 '2. 1 ho i|ii:iiiiii \ .\ is I ho nnmlier ol terms in tlie series, in this analysis. A' was arbitrarily chosen as 100. 

1 he iioniiiiiieiisional ((iiantilies are defined by 


r 

j* = — 


r =r — 
in 


n here /,,, and are the integral length aiul time scales of turbulence, it" is the friction velocity, and u' is 
ilie rms fliiciiiating llnid velocity which Kraiclman [10). assumed to be isotropic. For this pseiulo-tnrbulent 
\o|ociiy fiehl ilie energy spectrum E(k) is given as 


K{k) = 10 ^ 1 ', 




I he analysis ol l.i an<l Ahniadi [8] |>rovides an improvement based on the well-known fact of anisotropy in 
.1 inrbiilonl shear layer In order to estimate the anisotropy, they suggested the modincatioii 

= ii'r.i.U). (Ifi) 

"here r,(v) are oin|)irically generated shape fniictions chosen to match the data of Kreplin and Ecklemann 
[II). reproduced in Figure |. 


D row Ilian Motion 

Small particles suspended in a fluid undergo random translational motion due to molecular collisions, 
rids phenomenon is referred to as Brownian motion. The Brownian motion leads to the diffusion of particles 
Ml arror<Ianrt* with Firk's |;nv 

, (ir 

'■ = -0— 117, 


T) - 





Ill dll'; fi'iiiiiila. (■ IS i|i«' particle concent rat ton. 7 is tlie particle Hnx. and D is tlie diffusion coefficient. The 

■ lilliisiv It \ IS ni\en as 

II 

(18) 


kJ 

I ) = -:< . 


1 In' iliHiision e(|iiaiinn can l e olilained Ky ilirecliy snhstilutiiiR Pick's law into the equation for mass conser- 
\ nil >11 I lie Ui'i'wiiian forces aciinji, on the particle may l>e motleled as a white-noise process. The Brownian 
mol ion ol particles is moo I' d liy 

// I* 

(ly) 


(iu 

— + Jl, = »(/) 
(It 


wit h 


i = 


:{«•//(/ I 


( \in r 

Also, nil] is a (ianssian white-noise process with a spectral intensity 

■2kTJ 


( 20 ) 


( 21 ) 


I here k IS the Holtzmann constant aiul T is the temperature. From Eq. (21) it follows that the power 
pecirnm of the particle velocity is j^iven as 


S.ii(- ) =1 //(-•) I"’ -S.nli. 


Here the system function. H(jj). is defined hy 


//(w) = r 


1 


I hits 


■^i«ii('*'’) — 


iui “t* 7 

2/(-7'7 


Tlll(w' -f 7-’) 


I he particle anlocorrelalion fniiction is defined as 


/?<...(’■) = + t ). 


( 22 ) 


(23) 


(24) 


(25) 


and may he found liy takiny; tlie inverse Fourier Iransforiu of the power spectrum, i.e., 

/?..„(r)= —(26) 
m 

The difTnsiviiy is then niven hy 

rk-TC 

n= //..„( r)7r = (27) 

./o .{JT/Id 

N'Idilion.al del.ads of the Hrownian motion analysis and the other a.spects of particle motion and transport 
may hi' fonii'l in /Mmiadi [7] 
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III ilii' Inllciwiiit 'ifi'iioii ilirt'i' it'st rases are oiilliiieil. tarli rase was rlioseii to validate one of the 
iiicliMdiial III -1 li.iiiisiiis involved with at iiios|)heric dispersion of solid particles. 


Tost Case 1: Toriiiinal V'oloeity 

file ei|nni inn i .f mol inn f>f n small part irle in a Kravitaiional field is given by e<piation (3) and is reiterated 

III !.• Inr rniiveiili’liri' 

r—— + I/' = + Tij. (2«) 

iff 

If the M rinr vi’lnrity of the tlnid is identirally zero, (i.e.. stagnant), then Eq. (28) reduces to 

(29) 


(/«»■ _ 
~ r “ 


<J 


I i|naiinn I'.'tM i« a linear first-order ordinary rliffereiitial erpiation with constant coefficients and a constant 
Iniriim fnii' iion .\n i iinaiion of this ly|ie has a simple solution, which in this case is given by. 


It' = rg I erp t — 1 


(30) 


where (('■ is ifie particle velocity, r is the particle rela.xation time, g is the acceleration due to gravity, and 
I IS lime As the solution approaches steady state the velocity of the particle “relaxes" to a steady velocity 
eipial to r;/ 1 his is commonly called the terminal velocity. Figure 2 shows a comparison of the computed 

niimenral and analytical nondimensional velocity as a function of time. In this graph the velocity of the 
particle has been normalized by the terminal velocity. The curves show the excellent agreement between the 
iheoreiii-al and tin- computed numerical resnlts. Three curves are shown to illustrate the effect of size on the 
ii'rmmal \i lociiy, and the relaxation time of the particle. Here it is clearly seen that a 1/nii particle relaxes 
III lermin.d velocity in less than 2 x 10“'’ .seconds, at which time it reaches a maximum velocity of about 
I-’ ■ I** 'in/.s A lli/tm particle takes slightly longer, approximately 4 x 10“^ seconds, and has a maximum 
falling rate of (I I) Im/.s. A 100pm jiarticle will not reach terminal velocity until 0.07 seconds and will fall at 
.1 rale of 0 ipim/.s From these results, it is easy to see why a Ip/n or a lOptii particle will stay suspended in 
the atmosphere for long periods of time, whereas a 100pm particle will most likely precipitate out quickly 
under the force of gravity. 


Test Case 2; Drowiiiaii Diffusion in a Rectangular Box 


In this test case the Brownian dispersion model is validated. The discrete particles are replaced by a 
■ oniinnons fmn iion representing the dispersed matter. The Eulerian equation of conservation of species is 
as 

'It ~ \(^ CV ) 


(31) 
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Ii !<; ,1 ili.ii th- soliiiioM of Kci (;U) aiul the l..iij,raiigian e(iua!ious of motion should give P(|uivaleii( 


i. -.lilts f..i II liirn' iniml"'r of iiiiimtely small (lartirles. 

for a point soiiir.- in a rertannular hox. the lime history of the dispersion process is given l>y the 
I' >ns<r\ai ton of spet i.-s e<|nation vvitli delta f\mcti<i« forcing, terms Sperifically. 


— - /> ( -—^ - •< ■li'lV- - t„) 

,)l \;h - <)y- i>:- } 

I >1 all liomon' iieons homulary conditions this eipiation has an exact Fourier series solution given by 

-EZZjxr"""—'’"“"T’ 

n »»j V * 

/ ?/■ T* nr T‘ /»• ?r* \ 

( “77" +—p—+ “TT ) IITJ- . mst/ . 

. rp ^ ' ■' —)sin(—-^)si 


_ _^ turr 

sm( —) sin( —j—) sin( ) 


(3.3) 


\li.i.- /, /., and / are the lennihs of each side of I he fx.x in tin' X-. z-direction respectively, and J',.. 
are llie \ t- and Z-posilimis of the release, respeclnely 


In older to approximate this solniion as closely as pcrssible. the simulation is done using an extremely fine 
mesh 1 he mesh is fell annular and lie- spacinn between nodes is 10//in. It should lie noted that although 
I he mesh spacinn is small, the (tarticle diameter of 0 tfl/im is still 1.000 times smaller than the mesh spacing. 

I he clmid resnitinn from an instantaneous release of d.ttOO particles after -’lOO x 10“^ seconds is shown 
in I inure ;i alonn "Hh the conipniational mesh In this test case the x-. y-. and z-directions are chosen 
IIIhitranly I inures 1 tlironnh Ii detail the particle dispersion in time. Figure 4 shows the cloud dispersion 
III the x-direrii(iii, I'xcelh'iil agreement is shown between ilie analytical and computational results at small 
I lines- howexer, .is lime increases the iiiiiiieriral sclieiiie oxerpredicis tlie analytical .solution slight ly. Figure 5 
shows the ihspeise.l , loiid concentration in the y-direclion. .Again, the computed numerical solution slightly 
"\eresiim.iies ihe dispersion of the cloud. Figure li shows the trend of overpredictioii in the z-directiou of 
I he chuid dispersion. 1 he o\erpredictiolis are on the order ofin the worst case. 


I he coini'iiie.l immerical results have heeii sliowii to maicii well with the analytical results. If an infinite 
iitimher of massless particles were released, and if small disturbances could reach infinity in.stantly, then 
I he 1 wo results wolild iiiaicli exactly Frrors m this calculation arise from the method used t- '-stimate the 
■ oiiceiiirat ion, anil the inaccuracy of the continmim assumption implicit in the F.ulerian relerence frame. 
Iiicideiitallv. the o\er-esiiiiiation of the dispersion gives a more con.servafive estimation of poihitant dosage. 


Ti'st Chs»> 3; Den«>siti«)i» Rate in a Tni'hiiU>iit Chauiud 

In the prexioiis lest cases, exact soiiition.s to t he eipiations that describe the gravitational sedimentation, 
and Hrownian dispersion were given No roni|>arisons with experimental data have been presented. Tlie 

.') - ,s 








\ .'ll II |;i I ion I il I 111' III r 1)11 It'll I ilis|)i‘rsioii ran no) so|v«*il aiialy i irally. There is also a lark of ex|>erin»eulal tla(a 
'li'iailing siilnl parin li' ilispi-rsion in a Iiirlnilent How. However, (here exists some experimental data relating 
I III- ili iiosii 11)11 rail'' of |iariii |es in iiirhiileut cliamiel flow Tlierefore. a test case was chosen to valulate the 
I. IM:U> I oiii> wiili I'xpenmi nial flata ami test the lurhnlent ilispersion model used. 


riie (low lii'ld iis.'d in (he test rase is a three-dimensional representation of a two-<limensional rliannel. 
\ ‘JI ■ 2!' • li lei iannular mesh, in the sireamwise. vertical, and spanwise directions, respectively, was useil 
III this tesi rase I lie mean flow tield in a turbulent channel wa.s given by Ahmadi et al. [12] as 


11 

(I,. 


1 . 2 co.sh(^' I/) , A-' cosh A' 

‘ -t sinhA' 

I + iL tosl) <)•' 7 1 

‘ ^ F sinliA' 


(:H) 


where t; = I - 1 / li is the iiondimeiisional di.stance from the centerline. /» is the channel half-height. A', and 
A' .III' IIIiiidimensioiial parameters deflned as 


(IotiToiu/f" 

1 -II 71/fr'’'"3' 




,, 1 -II TlZ/r" 

' “ 0.71 


Here. II. Is the mean centerline ve|ocii\ ami is given by 


(36) 


r 

- 0,71//,•)■>!;« 


(37) 


where I IS the mean channel velocity ami 




(38) 


IS the He\Holds number, which in this case was set e(|ual to 17,700. The quantities p and p are the fluid 
d' lisiiN and dvnamic viscosity, respectively. Figure 7 shows the mean velocity profile for a 2r).4imii channel 
with a mean velocity of-'mi/.s At this Reynolds number the friction velocity is 0.3>»/s. The anisotropy 
I I the I III liiileiit iliici tint Ions of the channel flow that was used in this ti'st ca.se is shown in Figure 8. The 
normal coniponeni of the flint nations near the wall have a profound effect on the deposition rate of particles. 
It IS wi'll known that i' ha.s a quadratic variation at short distances from the wall, i.e.. 


t' = M' (.39) 

w here .1 IS a constant empirically chosen to be tj.()12 for this study. 

f igiire M shows t he nomlimensional ileposition velocity of particles as a function of nondimensional particle 
I' laxalion imw for a niuform concentration ol ( „ near the wall of the channel, the deposition velocitv is 
||<•fined as 

u,i = .//<"„ (40) 
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I IS tlir i.mukI.- tinx i" lli»' "all por null lime. The iioiuliiiiensioiial il^positiou velocity w given as 


II* = ii^/ ii~. (•! 1 ) 

1 li.- iK'iKlinn iisional parlid.- i>laxali<>ii liiiu' is defined lo lie 

^ (42) 

IXl/- 

lii I inure '.I, the present simulation results are compared with the accumulation of experimental data as 
produced III Pap.-uergos and Hedlev [13] Although there is snhstantial scatter in tiie data, the present com¬ 
pute.1 results agree wtihiti the scatter of the data. Fhe predictions are on the high edge of the experimental 
■ lata .'iisuritig a .niiservait\• estimate of deposition rate. 


C'oiicliisious 

1 h.' coiiipui. .1 results pr. ill. te.l hy the L I .\l3Dsolid paritciilale dispersion rode ha,s been compared to the 
i. siilts from three ive||-|,elia\ei| test ra.ses I he lermnial velocity. Brownian motion, and turhnlent dispersion 
models ha\e l>. - II li sted using two analytical and one experimental tests. The terminal velocity test case harl 
• xcell.'iil agre. iii.'iii li.'iweiMi the analytical atul crunputed numerical results. An error of less than 1*^' was 
shown h. iwi'en the two solutions The Brownian diffusion test case had slightly more error, approximately 
a'a'. lull it loo oHi-red sigiiilicaiil ronfiiieiice in the method. I'he deposition velocity calculation was shown 
to 111 in good agreement with the experimental <lata. although the data has .significant scatter. These results 
;in‘ very I'ncouiaging ami pim ide a higher degree of confidence in the code's predictive capability. 

Suggostoil Futuia' Work 

I he next step in the d.’velopment of the predictive capability should be a study of the spatial and 
temporal resolution that is i.’.ptired to resolve the st»cond-orrler rlifferences in the particle motion versus the 
passiv.' scalar lueiho.l I he inclusion of temperature and density variations shonhl be integrated into the 
mo.lel. .\uoilii r lonn al miprovenieni fi>r future capability is to integrate the LTM3D coile and the Lawrence 
l.iv.-rmore N.aiional Laboratories FF,M3A/13 code lo solve for ga.seons, lirpthl. aiul solid particulate dispersion 
■«miult aneoiisly Ihe combustion aiul chemical reaction of solid particles also may need to be accounted for. 
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A MODEL BASED PROGRAMMING ENVIRONMENT 
FOR PARALLEL IMAGE PROCESSING 


Michael S. Moore 
PhD. Candidate 

Department of Electrical Engineering 
Vanderbilt University 


Abstract 

This report presents the results of the 1993 AFOSR-RDL Graduate Student Summer 
Research Program. The eventual goal of this continued research project is to create a 
high performance, model-based parallel image processing system with a graphical model¬ 
ing environment. The main components needed to implement such a system are a modeling 
paradigm (language), a graphical model building environment, a model interpreter, a par¬ 
allel run-time environment, and an image processing application library. 

Parallel Processing Modeling Language (PPML) and the PPML interpreter provide a 
modeling environment that is adequate for modeling image processing algorithms. Khoros, a 
widely used image and signal processing system, has a graphical programming environment 
which can be used as a graphical model building environment. Khoros also provides a 
large, well developed library of image processing algorithms which can be re-used. In order 
to achieve the desired performance, we chose to use Multigraph as the parallel run-time 
environment. Multigraph is a parallel system integration tool developed at Vanderbilt 
University that allows the building of complex parallel applications from simpler processing 
blocks. By using the image processing libraries and gr'.phical programming environment 
of Khoros in conjunction with the PPML modeling paradigm and the Multigraph parallel 
execution environment, we have created a model-based parallel image processing system 
which features a graphical programming environment and rich image processing library due 
to Khoros and concurrent execution due to Multigraph. The system has been tested for a 
large variety of algorithms and data-flows, and is considered a success. 
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Introduction: 


At Arnold Engineering Development Center, there is a need for an image processing 
system with the high performance of parallelism as well as the flexibility and simplicity 
of visual programming. Usually, these requirements tend to produce contradicting results. 
W hile parallel implementations of algorithms exibit superior run-time performance, they 
are generally much more complex and require more programming time. It would be highly 
desirable to create an image processing programming environment which allows the user 
to create parallel algorithms while insulating him (or her) from the complexity involved in 
writing parallel programs. W'e have chosen to take a model-based approach to creating just 
such a programming environment. The various aspects of the algorithm are described by 
the user in the form of graphical pictures, or models. In the next stage, these models pass 
to an interpreter, which synthesizes the software necessary to implement the algorithm. 
The application is then sub-divided into parallel tasks, which are allocated to a hardware 
network and executed under the control of a parallel execution environment. The maiin 
components necessary for this system are (see figure 1) 

• graphical model building environment 

• .An image processing modeling language and model interpreter 

• A parallel execution environment 

• A library of image processing functions 

Khoros. a widely used image and signal processing system, has a graphical programming 
environment called Cantata in which the user creates image and signal processing appli¬ 
cations by drawing pictures that model the signal flow of the algorithm. This graphical 
programming environment is easy to use and has become very popular. For this reason, 
we chose to use Cantata as the graphical modeling builder for the new system. Also, we 
chose to take advantage of the large Khoros image processing library, which contains a 
rich assortment of commonly used algorithms. Re-using the Khoros programs and libraries 
has saved huge amounts of development time. Since Khoros does not contain facilities to 
efficiently parallelize algorithms, its run-time performance is not adequate. Thus, we chose 
Multigraph as the parallel execution environment. Multigraph, a parallel system integra¬ 
tion tool developed at Vanderbilt University, allows the building of complex algorithms from 
simpler processing blocks. It provides the underlying control and communication necessary 
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System Components 




Figure 1: A Model-Baised Parallel Image Processing System 

for parallel execution, and is capable of efficiently distributing processes across a non- 
homogeneous network containing a variety of computer architectures. It has been shown 
that with the use of Multigraph, the natural parallelism inherent to image processing can 
be exploited, increasing dramatically the run-time performance of many algorithms [4] [5] 
[8] [6] [7], Parallel Processing Modeling Language (PP.ML) and the Multigraph compatible 
PPML interpreter, which were developed at Vanderbilt University, were used as the image 
processing modeling environment for the system. 

To achieve a hybrid Khoros-Multigraph system, a two level interface was built. Level 1 is 
the model level. A model translator interfaces the Cantata visual programming environment 
to the PPML model interpreter. Level 2 is the program level. A program converter acts as 
an interface between the Khoros image processing programs and the Multigraph execution 
environment. 

The working system consists of the following components (refer to figure2) : 

1. The Cantata Visual Editor 
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2. The Cantata to PPML model translator 

3. The PPML model interpreter 

4. The Khoros program to Multigraph actor converter 

5. The Khoros image processing library 

6. The MGK parallel execution environment 


Khoros-Mgk System 


I Cantata Visual 
Editor 


I Khoros Imago 
Processirtg Libraries 


I Stand-Alone 
Khoros Programs 


Khoros to PPDL 
Translator 



Khoros Program to 
Kactor Converter 




Figure 2: Khoros-Multigraph Parallel Image Processing System 

By interfacing Khoros with Multigraph, we have created an image processing system that 
has the run-time performance of Multigraph as well as the rich image processing library 
and graphical programming environment of Khoros. 

In the following text, an explanation of the objectives of the summer research is given. 
Then follow descriptions of the image processing modeling paradigm and the model inter¬ 
preter. Next, the details of the Khoros-Multigraph interface are described, including the 
model level interface and the program level interface. Lastly, the results are discussed and 
recommendations for future research are made. 
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Objectives Of The Research Effort; 


The basic research objectives were: 

• Create a program level interface between Khoros and the Multigraph execution eni^- 
ronment. 

• Create a model level interface between the Cantata visual programming environment 
and the PPML model interpreter. 

• Integrate the system to create a model-based parallel image processing programming 
environment. 


The Program Level Khoros-Multigraph Interface: 

The reason for interfacing the Khoros image processing functions to the Multigraph 
kernel was to enable us to take advantage of the vast image processing libraries of Khoros 
without having to rewrite the programs. By re-using the Khoros programs, we have avoided 
wasting time and resources to “re-invent the wheel". Before describing the program con¬ 
verter, it is necessary to understand some of the details of both the khoros image processing 
environment and the Multigraph kernel. 

The Khoros Image Processing Environment 

The Khoros image processing environment includes a large array of commonly used 
image processing algorithms. Each of these algorithms is implemented by a stand-alone C 
program which can be executed from the command line or from within the Cantata visual 
programming environment. With either execution method, arguments are passed to the 
programs through a command line string. The Khoros programs read and write image files, 
the names of which appear after keywords in this command string. In Cantata, the image 
files provide the inter -program communication necessary to implement the data-flow. Each 
connection in the data-flow graph corresponds to an image file with a unique filename. 

The structure of the Khoros programs is very regular. Each program has a customized 
argument parser which constructs a program specific argument data structure from the 
command string. This structure contains algorithmic parameters, flags, and input/output 
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filenames. The program checks the validity of the arguments and ensures that the input 
files e.xist. The input images are read from the files and checked for data type compatibility 
with the algorithm. The program then makes a call to the appropriate image processing 
library function with the algorithmic parameters and flags. On return, it writes the the 
output image(s) to a file(s). Thus, each program is merely an argument parsing and error 
checking mechanism that serves a.s an interface to an image processing library function. 

The Multigraph Environment 

The computational unit in the Multigraph environment is called an actor. Each actor 
has an I/O interface, a script, and a context (see [1]). The I/O interface consists of two kinds 
of ports, input and output ports, through which the actor can receive or propagate data. 
The port communication is implemented by an underlying message passing mechanism of 
the kernel. The script is the code that will be executed when the actor is activated, and is 
usually on the subroutine level. The execution of an actor script may have either an IFANY 
or an IFALL triggerring mode. An actor with IFALL triggering mode will be triggered only 
when data is present at all of its input ports, while an actor with IFANY triggering mode 
will be activated if data is present at any of one its input ports. The actor’s context is a 
static memory area available to the script which is protected by the kernel. It can be used 
for passing parameters to the script or for storing state information. 

The Kactor Program Interface 

My approach was to create a protocaU by which any khoros stand-alone program could 
be converted to a Multigraph actor. Pains were taken to make the converter work without 
any required editing of the khoros program sources, so as to minimize the depth of under¬ 
standing and effort required to install and maintain the system. The program conversion 
utility is called “kactor”, for khoros-actor. Conversion is performed by recompiling the 
un-edited khoros source programs with the macros in “kactor.h” included and then Unking 
to “Ubkactor.a” (in addition to the normal khoros support Ubraries). To include the macros 
without editing the source files, I copied the files “vinclude.h” and “xvinclude.h” to the 
kactor directory and added the foUowing Unes to them; 

#ifder MGK.VERSIOH 
tinclude "kactor.h" 
tendil 

Since either vinclude.h or xvinclude.h is included by every khoros program, this technique 
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does the trick. The kactor macros perform the following: 

1. '‘main" is replaced by a standard IFAL*. actor script which first calls initialization 
fuction called ’’kactor-init’’. kactor-init parses the strings contained in the context 
structure into the argument vector used by the khoros program and creates tables used 
by mgk-getimage and mgk-send Jmage to map filenames to I/O port identifiers. The 
standard actor context structure includes three strings: (1) a command string (2) an 
input file keyword list (3) an output file keyword list. The original program is changed 
into a subroutine, which is called by the actor script. The names of the script and the 
subroutine are created via macros that append “.script** and “subroutine” onto the 
base of the filename. For example: vconvolve.c contains the script “vconvolvescript” 
and the subroutine “vconvolvesubroutine". 

2. The file I/O is changed to port I/O by replacing each call to “readimage” with a call 
to “mgk.getJmage”, and each call to “writeimage” with a call to “mgksendJmage”. 
The mgk.get Jmage and mgksendJmage routines receive/propagate images from/to 
the correct port by referring to the forementioned filename to I/O port conversion 
tables. 

The kactor script behaves exactly as the stand-alone program did. The argument parsing, 
input checking, and error handling facilities are left intact. The main differences are (1) 
instead of reading and writing files, the kactors receive and propagate image data using the 
Multigraph message passing mechanism. (2) The program is executed as an actor under 
the control of the Multigraph kernel instead of as a stand-alone program under the Unix 
shell. 

The following khoros programs were successfully tested as kactors. 

• vconvert (converts image data types) 

• vconvolve (convolves an image with an arbitrary kernel) 

• venhance (performs pixel re-mapping based on local statitistics) 

• vnormal (normalizes an image) 

• vsub (subtracts two images) 

• vhstr (performs a histogram stretch) 
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• putimage (displays an image with update) 

• vgauss (creates a gaussian image) 

• vgcirc (creates the image of a filled circle) 

• vrtoc (converts reaJ data to complex) 

• vbandcomb (combines bands from m images into one multiband image) 

• vtranspos (transposes an image) 

• vpad (pads a data sequence with a constant) 

• vbandsptl (extracts a specified data band from a multiband image) 

• vfft (performs a 2-D fft on an image) 

• vbrf (applies a band reject filter) 

• vraoredge (performs a morphological edge extraction) 

• vabs (computes the absolute value of an image) 

The kactor conversion technique is applicable to any khoros program, so the new system 
can take advantage of the entire image processing library of Khoros. 


The PPML Modeling Language: 

PPML stands for Parallel Processing Modeling Language. It was developed at Vander¬ 
bilt University by Dr. Gabor Karsai [2]. PPML is a declarative language that includes the 
facilities to model three aspects of a processing network: 

• The processing (data flow) aspect 

• The hardware aspect 

• The allocation of sub-tasks to the hardware. 

The PPML models are made up definitions of the following types; 

• primitive 
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• compound 

• node 

• network 

• configuration 

Primitive and compound definitions form the processing aspect model. Node and network 
declarations comprise the hardware aspect model. The configuration definition models the 
system as a whole, including the sub-task to processor allocation. 

.\ primitive is the elemental processing block, and is described by a PPML “defprimitive” 
declaration. This declaration contains a primitive name, a triggering mode, an I/O interface 
description, a context description, and a script. Each of the primitive definitions represents 
a separate kind of process, or actor, that can be used in the signal flow. A compound is 
a group of primitive instances which are interconnected *. Compounds are declared by a 
“defcompound*’ block in PPML. A defcompound block contains a name, an I/O interface 
description, a list of signals (data nodes), and a structure definition. In the structure 
definition the internal connectivity of the compound is described, and the context of each 
primitive instance is set. The node and network models describe the hardware. A defnode 
declaration describes a particular type of hardware element (or node), such as a Sun spare 
station or a TIC40 processor. The defnet declaration describes a hardware network which 
is made up of instances of the defined node types. The configuration model is represented 
textually in a defconfig declaration, which includes an instance of a compound (processing), 
an instance of a network (processors), and an allocation (a mapping of the processing blocks 
to the network nodes). A more complete description of PPML will be given in a user manual 
in the future. 

.4n example PPML script is shown below. This script was created using the Khoros- 
PPML translator that is described in the next section. It represents a signal flow that was 
created in Cantata. Note that the first part, up until the defnode declaration, represents 
the signal flow. The signal flow is made up of 10 actors of 7 varieties, and 7 signals which 
connect them together. The defnode and defnet blocks define a hardware network made up 
of two "’spare’’ nodes, called "’prod” and ’’proc2’’. The defconfig declaration is a mapping 

'Since PPML is a hierarchical language, the structure can also include instances of other compounds. 
The top-level compound is distinguishable by the fact that it will have no inputs or outputs in its I/O 
interface description. 
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of the primitive instances to the processing nodes in the network. Note that all of the actors 
in the graph are allocated to prod. This is the default of the translator. At this time, the 

allocation of parallel tasks to hardware must be done manually. 

; **♦**♦*♦***•♦•***•♦•♦•••***•♦**•**••*•*•**•*•••*♦•♦*•♦*••♦*•**•••♦•♦ 

(defpriaitive »convert_typ«l ifall ((woidp inputl) -> (voidp oatputi)) 

((string comnand) (string input.keys) (string output.keys)) 

(script *• extern void vconvert.scriptO; vconvert.scriptfkcoBBand); ")) 

(defpriaitive putinage.typel ilall ((voidp inputl) -> ) 

((string command) (string input.keys) (string output.keys)) 

(script "extern void putimage.scriptO; putimage.script(kconmaad); ’*)) 

(defprinitive vconvolve.typel ilall ((voidp inputl)(voidp input2) -> (voidp outputl)) 
((string command) (string input.keys) (string output.keys)) 

(script " extern void vconvolve.scriptO; vconvolve.scriptfkcommand); ")) 

(delprimitive vsub.typel ifall ((voidp inputl)(voidp input2) -> (voidp outputl)) 
((string command) (string input.keys) (string output.keys)) 

(script " extern void vsub.scriptO; vsub.scriptfkcommaad); ")) 

(defprimitive vnormal.typel ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.keys) (string output.keys)) 

(script " extern void vnormal.script(); vnormal.script(kcommand); ")) 

(defprimitive readimage.typel ifall ( -> (voidp outputl)) 

((string command) (string input.keys) (string output.keys)) 

(script " extern void readimage.scriptO; readimage.script(kcommand); ")) 

(defprimitive vriteimage.typel ifall ((voidp inputl) -> ) 

((string command) (string input.keys) (string output.keys)) 

(script " extern void vriteimage.scriptO; vriteimage.script(kcommand); ")) 

(defcompound ultra (->) () () 

(signals (voidp dnl) (voidp dn2) (voidp dn3) (voidp dn4) 

(voidp dnS) (voidp dn6) (voidp dn7)) 

(struct 
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(vconv«rt_typ*l_nu«l (vconvert_typ«l ( dnl -> <in2) 

('*vconT«rt -i ./ultrasound.xv -o ./AAA -t \"lloat\" -b 0 ** "-i” ’*-o"))) 
(putiaags.typsl.nual (putinage.typsl ( dnl -> ) 

("putinage -i ./ultraaound.xv -updats 2 " "-i" ”"))) 
(vconvol»s_typsl_nual (vconvolva_typ«l { dn2 dn3 -> dn4) 

("vconvolve -il ./AAA -12 ./adgakam -o ./AAB " "-11 -12" "-o"))) 
(vsub_typal_nuBl (vsub_typsl ( dn2 dn4 -> dnS) 

("vsub -11 ./AAA -12 ./AAB -o ./AAC " "-11 -12" "-o"))) 
(vnonBal.typal.nuBl (vnonal.typal ( dnS -> dn6) 

("vnonial -1 ./AAC -o ./AAD -n 2SS " "-1" "-o"))) 

(vconvert_typ«l_nun2 (vconvsrt.typal ( dn6 -> dn7) 

("vconvart -1 ./AAD -o ./output.xv -t \"byta\" -b 0 " "-1" "-o"))) 
(putlmags_typ«l_nuBi2 (putlmaga.typsl ( dn7 -> ) 

("putImage -1 ./output.xv -update 2 " "-1.’))) 

(readlmaga_typel_numl (readlmage_typel ( -> dnl) 

("readlmage -o ./ultraaound.xv.-o"))) 

(readlmage_typel_nun2 (readlmage.typel ( -> dn3) 

("readlmage -o ./edgekem" "" "-o"))) 

(vrltelmage_typel_numl (vrlteimage.typel ( dn7 -> ) 

("writeimage -t ./output.xv" "-i.))) 

)) 

(defnode spare () 

(speed 1.0 0.001) 

(sizes (char 1 1) (uchar 1 1) 

(short 2 2) (ushort 22) 

(Int 4 4) (ulnt 4 4) 

(float 4 4) (double 8 8) 

(string 4 4) (voldp 44)) 

(queueing Icfs) 

(cc "gcc -g -c -I./ *a") 

(Id "bulld-khoros 'a serv") 

(sh "(./serv 8000 t); sleep 1") 

(Ident (values "alex" 0))) 

(dsfnet proc () 

(Internals ethers) 

(Init ) 
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(parts 

(prod (spare ())) 
(proc2 (spare ())))) 


(defconfig test () 

(proeassing (ultra (ultra ()))) 
(proeessors (spa (proc ()))) 
(alloeatloa 

(ultra readiaaga_typ«l_aual => 
(ultra raadiBiaga_t 7 psl_aum 2 => 
(ultra vritsiBaga.typal.numl => 
(ultra veonvert_typ«l_numl => 
(ultra Vconvert_typel_num2 => 
(ultra vconvolva_typal_nuiBl => 
(ultra vsub.typel.nual => 

(ultra vnonBal_typel_numl => 
(ultra putimags.typel.nuBl => 
(ultra putinage_typal_num2 => 


spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 
spa prod) 


)) 


; *•••♦♦♦•**•**♦••*•♦•**♦♦*♦♦•♦*♦**♦♦*♦♦•*•*••********•** 


The Model Level Khoros-Multigraph Interface: 

The model level Khoros-Multigraph iterface consists of a model translator that creates 
a PPML declaration from a Cantata data-flow model. The Cantata data-flow model is 
obtained by capturing the standard output of Cantata as the graph is executed, which 
consists of the exact commands that the user would type if he (or she) were executing the 
data-flow step by step from the command line. This information is adequate to reconstruct 
the signal flow of the system, since it contains the program names, the algorithmic parame¬ 
ters and flags, and the filenames (the filenames specify the interconnectivity of the graph). 
The translator is a stand-alone C program which reads in the logged standard output from 
Cantata and writes the corresponding PPML declaration. An example of the operation of 
the translator will be given in the next section. 
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Results: 


The Khoros-Multigraph image processing system was tested for a number of different 
data-flows. Each of the khoros functions mentioned above were tested in at least one graph. 
An example is given below that demonstrates the operation of the system. Note first the 
Cantata data-flow shown in figureS. This graph, which is an example released with the 
standard distribution of Khoros^, perfoms an edge detection on an artificially generated an 
image and displays the result. As the graph in figureS was executed, the standard output 
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Figure 3: Example Cantata Data-Flow 

of Cantata was captured and piped to a log file, which was then translated into a PPML 
declaration. Both the log file and the PPML models are shown below. 


^This workspace can be found in KHOROSJIOME/repos/workspace/vipLwksp/vipl.Z 
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««*««•**«««««*«**«*««««•««**««««•«*««•*«««*««««*•«*«•««•*«*«««•«**«« 


; This is ths standard out Irom Cantata 

; m******************************************************************* 

uncomprsss -c . ./khoros/repast^rkspacs/vipl.Bksp/Tipl.Z 

vgauss -o ./vgausAAAaOOlss -c 150 -r 128 -x 64 -y 64 -?x 30 -vy 30 -ra 0 -aap 1 -n 0 
vgcirc -o ./vgcirCAAaOOlss -d 50 -x 64 -y 64 -t "byte" -b 0 -f 255 -r 128 -c 150 
vrtoc -il ./vgausAAAaOOlss -o ./vrtocBAAaOOlss -t 1 

vconvert -i ./vrtocBAAaOOlss -o ./vconvDAAaOOlss -n 255 -t "byte" -b 0 
vbandcomb -il ./vgcirCAAaOOlss -i2 ./vconvDAAaOOlss -o ./vbandLAAaOOlss -b 0 -c 0 
vtranspos -i ./vbandLAAaOOlss -o ./vtranNAAaOOlss 

vpad -i ./vtranMAAaOOlss -o ./vpadXIAAaOOlss -r 256 -c 256 -d 0 -e 0 -j 0 -k 0 
vbandsptl -i ./vpadXIAAaOOlss -o ./vbandOAAaOOlss -b 0 
vhstr -i ./vbandOAAaOOlss -o ./vhstrEAAaOOlss -p 0 

vconvert -i ./vhstrEAAaOOlss -o ./vconvGAAaOOlss -n 255 -t "float" -b 0 

vfft -il ./vconvGAAaOOlss -ol ./vfftFFAAaOOlss -d 0 

vbrf -i ./vfftFFAAaOOlss -o ./vbrfNHAAaOOlss -n 1 -c 0.5 -h 10 

vfft -il ./vbrfMHAAaOOlss -ol ./vfftllAAaOOlss -d 1 

vconvert -i ./vfftllAAaOOlss -o ./vconvKAAaOOlss -n 255 -t "byte" -b 0 
vBoredge -i ./vconvKAAaOOlss -k ./cell20.Bor -o ./vaoreJAAaOOlss 
putiaage -i ./vaoreJAAaOOlss -x 700 -y 0 -update 2 

•^^*^i,*m***m****************m*^************************** 

; This is the translated PPNL declaration 

(defpriaitive vgauss.typel ifall ( -> (voidp outputl)} 

((string coaaand) (string input.keys) (string output.keys)) 

(script "vgauss.script")) 

(defpriaitive vgcirc.typel ifeill ( -> (voidp outputl)) 

((string coaaand) (string input.keys) (string output.keys)) 

(script "vgcirc.script")) 

(defpriaitive vrtoc.typel ifall ((voidp inputl) -> (voidp outputl)) 

((string coaaand) (string input.keys) (string output.keys)) 

(script "vrtoc.script")) 
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(dsfpriaitiva Tconv«rt_typ«l ilall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.koys) (string output.keys)) 

(script "vconvert_script")) 

(dsfprimitivs vbandcomb.typsl ifall ((voidp inputl)(voidp input2) -> (voidp outputl)) 
((string command) (string input.kays) (string output.keys)) 

(script "vbandcomb.script")) 

(delprimitive vtranspos.typsl ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.keys) (string output.ksys)) 

(script "vtranspos.script")) 

(defprimitive vpad.typsl ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.keys) (string output.keys)) 

(script "vpad.script")) 

(dafprimitivs vbandsptl.typal ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.ksys) (string output.kays)) 

(script "vbandsptl.script")) 

(dafprimitiva vhstr.typal ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input.kays) (string output.kays)) 

(script "vhstr_script")) 

(dafprimitiva vfft_typal ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input_kays) (string output.kays)) 

(script "vfft.script")) 

(dafprimitiva vbrf.typal ifall ((voidp inputl) -> (voidp outputl)) 

((string command) (string input_kays) (string output.kays)) 

(script "vbrf.script")) 

(dafprimitiva vmoredga.typal ifall ((voidp inputl)(voidp input2) -> (voidp outputl)) 
((string command) (string input.keys) (string output.kays)) 

(script "vmoredge.script")) 

(defprimitive putimage.typal ifall ((voidp inputl) -> ) 

((string command) (string input.keys) (string output.kays)) 
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(script ''putiaags.script")) 


(defpriaitive readiaage.typsl Kail ( -> (voidp outputl)) 

((string couiand) (string input.ksys) (string output.ksys)) 

(script "rsadiBags_script“)) 

(dsfpriaitiva vritaimags.typsl ifall ((voidp inpntl) -> ) 

((string coBusand) (string input.ksys) (string output.ksys)) 

(script "sritsiBags.script")) 

(dsfcompound tost (->) () () 

(signals (voidp dnl) (voidp dn2) (voidp dn3) (voidp dn4) 

(voidp dnS) (voidp dn6) (voidp dn7) (voidp dnS) 

(voidp dn9) (voidp dnlO) (voidp dnll) (voidp dnl2) 

(voidp dnl3) (voidp dnl4) (voidp dnlS) (voidp dnl6)) 

(struct 

(vgauss.typsl.nuBl (vgauss.typsl ( -> dnl) 

("vgauss -o ./vgausA&AaOOlss -c 150 -r 128 -r 64 -y 64 -vx 30 

-vy 30 -ra 0 -aap 1 -n 0 .-o*'))) 

(vgcirc.typsl.nual (vgcirc.typsl ( -> dn2) 

("vgcirc -0 ./vgcirCAAaOOlss -d 50 -x 64 -y 64 -t byts -b 0 

-f 256 *r 128 -c 150 . "-o"))) 

(vrtoc.typsl.nuBl (vrtoc.typsl ( dnl -> dn3) 

("vrtoc -il ./vgausAAAaOOlss -o ./vrtocBAAaOOlss -t 1 ” "-il" "-o"))) 
(vconvsrt.typsl.nuBl (vconvsrt.typsl ( dn3 -> dn4) 

("vconvort -i ./vrtocBAAaOOlss -o ./vconvDAAaOOlss -n 255 

-t byts -b 0 " "-i" "-o"))) 

(vbandcoBb.typsl.nuBl (vbandcomb.typsl ( dn2 dn4 -> dnS) 

("vbandcoBb -il ./vgcirCAAaOOlss -i2 ./vconvDAAaOOlss -o 

./vbandLAAaOOlss -■ 0 -c 0 " "-il -i2" "-o"))) 
(vtranspos.typsl.nual (vtranspos.typsl ( dn5 -> dn6) 

("vtranspos -i ./vbudLAAaOOlss -o ./vtranNAAaOOlss " "-i" "-o"))) 
(vpad.typsl.nual (vpad.typsl ( dn6 -> dn7) 

("vpad -i ./vtranNAAaOOlss -o ./vpadXIAAaOOlss -r 256 -c 256 

-d 0 -a 0 -j 0 -k 0 " "-i" "-o"))) 
(vbandsptl.typsl.numl (vbandsptl.typsl ( dn7 -> dn8) 

("vbandsptl -i ./vpadXIAAaOOlss -o ./vbandOAAaOOlss -b 0 “ "-i" "-o"))) 
(vhstr.typsl.nuBl (vhstr.typsl ( dn8 -> dn9) 
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("vhstr -i ./vbandOAAaOOlas -o ./vhstrEAAaOOlss -p 0 " "-i" "-o"))) 
(vconv«rt_typel_nuB2 (vconv«rt_typ«l ( dn9 -> dnlO) 

("vconvart -i ./vhstrEAAaOOlss ~o ./vco&vGAAaOOlss -n 25S 

-t float -b 0 •' "-i" ••-o"))) 

(vfft_typel_nual (vlft_typal ( dnlO -> dnll) 

("vtft -il ./vconvGAAaOOlss -ol ./vfftFFAAaOOlss -d 0 " "-il" "-ol"))) 
(vbrf_typel_nual (vbrf.typel ( dnll -> dnl2) 

("vbrf -i ./vfftFFAAaOOlss -o ./vbrfNHAAaOOlss -n 1 -c 0.6 

-9 10 •• "-i" "-0"))) 

(vfft_typal_nu«2 (vfft_typel ( dnl2 -> dnl3) 

("vfft -il ./vbrfMHAAaOOlss -ol ./vfftllAAaOOlss -d 1 " "-il" "-ol"))) 
(vconv6rt_typal_nun3 (vconvart_typel ( dnl3 -> dnl4) 

("vconvart -i ./vfftllAAaOOlss -o ./vconvKAAaOOlss -n 266 -t byte 

-b 0 " "-i" "-0"))) 

(vBoredga_typal_nual (vmoredge.typel ( dnl4 dnl6 -> dnl6) 

("VBoredge -i ./vconvKAAaOOlss -k ./call20.Bor -o ./vaoreJAAaOOlss " 

"-i -k" "-0"))) 

(patiaage.typel.ntial (putiaage.typel ( dnld -> ) 

("putiaage -i ./vaoreJAAaOOlss -x 700 -y 0 -update 2 ” "-i.))) 

(readiaage.typel.nual (readiaa.ge_typel ( -> dnl6) 

("readiaage -o ./cell20.aor.-o"))) 

(writeiaage.typel.noal (vriteiaage.typel ( dnl6 -> ) 

("ariteiaage -i ./vaoreJAAaOOlss" "-i" ""))) 

)) 

(defconfig test.config () 

(processing (test (test ()))) 

(processors (spa (proc ()))) 

(allocation 

(test vgauss_typal_nuBl => spa prod) 

(test vgcirc.typel.nual => spa prod) 

(test vrtoc.typel.nual =:> spa prod) 

(test vconvert.typel.nuBl => spa prod) 

(test vbandcomb_typel_nual => spa prod) 

(test vtranspos_typel_nual => spa prod) 

(test vpad_typel_nual => spa prod) 

(test vbandsptl.typel.nual => spa prod) 

(test vhstr_typel_nual => spa prod) 
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(t«st vcon»*rt_typ«l_num2 => spa prod) 

(tsst vfft_typ«l_iiuail => spa prod) 

(tsst vbri_typ«l_auml => spa prod) 

(test vfft_typ«l_nu«2 => spa prod) 

(test vconvert_typel_nun3 => spa prod) 

(test v«oredge_typel_nuiBl => spa proci) 

(test putimage.typel.nunl => spa prod) 

(test readiaage.typel.nuBl => spa prod) 

(test sriteiaage.typel.nual => spa prod) 

)) 

;••«•«*«««*««*«*••«*«•**«*«««««««•«««*«••«««*«***««**««* 

The PPML models shown above were loaded into the model interpreter and the graph 
was executed by the Khoros-Multigraph run-time system. Note that each of the Khoros 
programs used in this Cantata workspace had previously been converted to Kactors and 
linked with the .Muitigraph kernel. The resulting output from this execution was identical 
to the output of the Cantata execution. The exe'-ution times were also very similar. On a 
Sun Sparc 10 running Solaris 5.2, the Cantata version took 31 seconds, and the Multigraph 
version took 29 seconds. This speedup is minimal, as expected. Only one processor was 
used in each case. 

By using Multigraph, we have gained the ability to easily parallelize the data-flow. By 
altering the allocation section of the configuration model, any actor of the graph can be 
executed on any processor in the network. Thus, it would be possible to distribute the 
example graph over a network of 10 processors without changing the result. It is important 
to note that this type of parallelization will not decrease the execution time of a graph 
unless the graph has a naturally concurrent structure. A serial data-flow wiU remain serial 
even if the sub-tasks are executed on different nodes. 

In order to guarantee an execution time speedup, a different approach is necessary. The 
image can be split into ‘“p” pieces, each of which can be processed concurrently by a separate 
processor. The processed pieces can then be merged back together to form the output. This 
technique is called “split-and-merge” parallelization, and has been shown to be effective for 
reducing the execution times of many of the commonly used image processing algorithms 
[4] [5] [6] [6]. The split-and-raerge paradigm can easily be added to the present system, 
which will greatly enhance its performance. 
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The new system will work as follows: The user will first build and test a data-flow in 
Cantata, as before. As an added step to the modeling process, he (or she) will designate 
split-and-merge processing blocks. The translator will then create a data parallel PPML 
data-flow model, leaving the number of parallel channels as an unspecified variable. During 
interpretation, the model interpreter will determine the number of parallel channels based 
on relavent factors (see [6]), and will automatically allocate the parallel channels to the 
available hardware. The execution time (latency) of the resulting spbt-and-merge data¬ 
flow graph will be much less than the execution time of the corresponding Cantata data-flow. 
Thus, with the addition of the split-and-merge processing paradigm, it is expected that 
the Khoros-Multigraph image processing system will greatly out-perform Khoros, but will 
be every much as easy to use 
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RECOMMENDATION’S FOR FUTURE RESEARCH: 

Although the Khoros-Multigraph system was successfully implemented, the real goals 
of its implementation have not yet been met. The execution time under the Multigraph 
network was marginally better while processing on one node. However, the goal is a high- 
performance system. The desired performance can be attained with the use of the split- 
and-merge processing paradigm, as has been shown by previous results [4] [5] [6] [6). The 
split-and-merge processing technique was not include merely due to the limited time span 
of the RDL summer research program. In the near future, facilities should be added to the 
system to include the this feature. 
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Adaptive Engine 
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Abstract 


The concept of system identification was discussed, as well as a specific 
identification program called SCIDNT (Systems Control Identification). An 
engine simulation program was also reviewed, specifically ATEST (Advanced 

Turbine Engine Simulation Technique). An adaptive simulation routine, 
whose origins were from a mainframe operating environment, was modified 

for use on a DOS operating system. The adaptive engine simulation routine 
was then evaluated using turbine engine test data from an engine test cell. 

It was concluded firom the evaluation that rn engine model trimmed using a 
system identification code results in an overall improvement in the prediction 
of engine measurements. 


Introduction 


Over the past several years, engineers at Arnold Engineering 
Development Center (AEDC) have been developing mathematical models to 
simulate various aircraft systems, including turbine engines. A primary use 
of these models is to examine engine response in operating regimes that are 
difficult or costly to examine in the engine test cell. Engineers are constantly 
establishing and evaluating new codes, as well as "tuning" existing codes to 
accurately reflect the specific component and overall engine performance. 

System identification is a technique that aids in the evaluation and 
"tuning" of mathematical models. System identification computer programs 
use a specified model and a given set of data, and then attempt to optimize 
the model to accurately reflect the data. The result is a set of parameters or 
"scale factors" which adjusts the model to reflect the data, or rejects the 
model as unfit to represent the physical system. 

The use of a system identification technique in identifying parameters 
within a turbine engine simulation model is explored within this report. 
Based upon actual engine data, the identification routine combined with the 
engine model is referred to as an adaptive engine simulation, due to the 
ability of the model to adjust to the current state of the engine. 

The following report begins with a background of the system 
identification code used, as well as a brief discussion of the mathematical 
model employed. Next, the objectives are defined and the approach taken to 
achieve the objectives stated. Finally, the results are presented, followed by 
recommendations for future work. 
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Background 

The system identification routine that is used is known as SCIDNT 
(Systems Control Identification), originally developed by Systems Control 
Technology Inc., under contract from AEDC. SCIDNT identifies, or 
estimates, model parameters to fit the model to experimental data based 
upon criteria specified by the user. 

The basics of system identification are examined in detail by Norton 
[1] and Eykhoff [2]. The system being modeled is an aircraft turbine engine. 
This type of system, when modeled, is referred to as a nonlinear, multiple- 
input-multiple-output system. Two common approaches to identif 3 dng the 
parameters of such a complex model are the equation error method and the 
output error method [3]. 

Both of the above mentioned methods estimate unknown parameters 
through the minimization of a performance index [3J. The equation error 
method, also called the least squares method, is effective at dealing with a 
d)rnaniic process. However, this estimation method does not deal well with 
disturbances in the measurements (i.e. measurement noise). Conversely, the 
output error method deals well with measurement noise, but does not deal 
well with process dynamics or "process noise". 

The SCIDNT identification routine utilizes an estimation method 
which combines the advantages of both methods. This method is a variation 
of the Gauss-Newton (GN) method known as the Levenberg-Marquardt (LM) 
algorithm. As with the equation error and output error methods, the LM 
algorithm within SCIDNT seeks to minimize a performance index. As 
previously mentioned, the LM algorithm is a variation of a GN method; the 
variation being that the algorithm is a combination of a gradient search and 
a GN method to search for a minimum. 
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The performance index that SCIDNT uses is a cost function associated 
with the error corresponding between the actual versus the simulated 
measurements. SCIDNT seeks to minimize the cost function through the 
modification of user specified model parameters. Model parameters 
corresponding to a minimized cost function are referred to as estimated or 
"optimized" parameters which produce a trimmed model. 

An important feature of SCIDNT that must be addressed is that 
weighting factors are associated with the measurements. Each measurement 
has a corresponding weighting factor that must be specified by the user. A 
modification of a particular weighting factor increases or decreases the 
impact of the corresponding measurement error on the overall cost function. 
A drop in the weighting factor corresponds to an increase in the impact; a 
nominal weighting factor setting is 1.0. Representative weighting factors, as 
well as the parameters and measurements, are shown in Tables A1 and A2. 

Setting all weighting factors to 1.0 does not imply that all 
measurements have the same importance in the determination of the overall 
cost function. With the current program design and nominal weighting 
factors, the thrust measurements have a dominant role in the determination 
of the cost function. This current design is simply a starting point and future 
modifications of the various weighting factors is anticipated. 

The turbine engine model that was used is known as ATEST 
(Advanced Turbine Engine Simulation Technique). This simulation code is a 
1-D, component-matching routine, capable of modeling arbitrary engine 
configurations. For the current research, a dual-spool, mixed-flow turbofan 
engine model is employed, illustrated in Figure 1. The engine model can 
simulate both steady-state and transient operation (0 - 20 Hz), but not 


7-5 






dynamic operation. A more detailed description of ATEST is given by 
Chappell [4], 
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Figure 1. Schematic of a duai-spool, mixed-flow turbofan engine. 


As a stand-alone program, ATEST will calculate a design point using 
user supplied input data. This design point establishes the scale factors for 
component performance maps and mixing plane areas. With a design point 
calculated, the model is in effect trimmed to a family of engines. Due to 
variations in engine characteristics (e.g. engine manufacturing, tolerances 
variations, and engine deterioration), the model output can not be expected 
to reflect exactly the performance for any particular engine. 
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Objectives 

The first objective was to convert the ATEST and SCIDNT FORTRAN 
programs from a mainframe operating environment to a DOS operating 
environment. The second objective was to modify the adaptive simulation 
code to establish a means of evaluating the code. The last objective was to 
perform a limited evaluation of the adaptive simulation code. 


Approach 

The ATEST model and the identification program SCIDNT were 
converted from a mainframe to a DOS environment. For discussion 
purposes, the combination of SCIDNT and ATEST, shown if Figure 2, will be 
referred to as the adaptive simulation process. Newly integrated code is 
shown in the figure with a single line border. 
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Previous versions of the ATEST program had been converted to a DOS 
environment, thus this conversion was used as a learning process for the 
more formidable task of conversion of the SCIDNT program. Due to the 
standardization of the FORTRAN programming language, the bulk of the 
effort in the conversion process was spent on tailoring the I/O (i.e. 
input/output) functions to the conventions of the particular compiler, in this 
case, the Lahey FORTRAN compiler. One major modification came in the 
method of transfer of the test data to the SCIDNT program. In the 
mainframe environment, the test data was transferred to the program in 
binary form. It was decided that in the DOS environment that the test data 
would be transferred in ASCII form, which allows easier data manipulation 
for the program user. 

Once the adaptive simulation code was operational, a modification was 
made to evaluate the performance of the code. The modification came in the 
form of an additional execution of the ATEST engine model with the 
estimated parameters from SCIDNT and the input test data. The resulting 
output from the additional execution are calculated measurements which can 
be compared to the actual test measurements. 

Finally, a limited evaluation was performed on the adaptive 
simulation program. This evaluation compared the actual measurements to 
the simulated measurements from the trimmed, engine simulation model. 
These comparisons are presented and discussed in the resvdts and 
conclusions sections. 
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Results 

A limited evaluation on the adaptive simulation code was performed. 
Data from a particular engine was used and comparisons were made 
between the trimmed model output and the actual measurements. Data &om 
the dual-spool, mixed-flow, turbofan engine are shown in Figures A1 to All. 

The evaluation began with an initial execution of ATEST, in stand¬ 
alone form, to calculate design point operation. This initial execution of 
ATEST creates a baseline from which the identification code can be evalu¬ 
ated and represents an engine model trimmed to a family of engines, but not 
to the particular engine tested. 

In Figures A12 & A13, a comparison of actual measurements versus 
the simulated measurements is shown. One can see that the measured and 
simulated net thrusts compare well, while the engine airflow measurements 
differ sigxuficantly. Assuming valid data, the reason for this difference is 
that the engine model is not trimmed to the particular engine firom which the 
data originated. 

The next step in the evaluation process is to examine the change in 
simulated measurements when SCIDNT determines the model parameters. 
This initial identification is the first step in creating an engine simulation 
trimmed to simulate a particular engine with its individual characteristics. 
Figures A14 & A15 show the comparison of the actual measurements with 
those simulated using the adaptive engine simulation model. For this 
comparison, the weighting factors are at there nominal settings. SCIDNT 
output from this execution of the adaptive engine simulation ; shown in 
Table AL. 

As shown in Figure A14, the simulated net thrust compares favorably 
with the actual measurements. This is due to the dominance of the thrust 


7-9 






measurements on the determination of the cost function using nominal 
weighting factors. The simulated engine airflow resulted in an improvement 
over the baseline, but an appreciable difference still exists. 

The next step in the evaluation was to examine the impact of 
modifying the weighting factors. Due to the limited evaluation, modification 
of only the airflow weighting factor could be evaluated. The modification in 
the weighting factor increased the importance of the airflow in the 
calculation of the overall cost function. In this execution, the airflow 
measurement is now on the same order of magnitude as the thrust 
measurements. Table A2 shows the SCIDNT output with the modified 
airflow weighting factor. 

The results of the modified airflow weighting factor show a significant 
improvement in the simulated engine airflow. This improvement is due to 
the fact that the increased importance of the airflow measurement "drove” 
the simulated values closer to the actual, as seen in Figure A17. However, 
agreement between simulated and measured net thrust was worse (Figure 
A16). The disagreement can be attributed to, in part, the decrease in 
importance of the thrust measurements in the determination of the overall 
cost function. 

Conclusions and Recommendations 

The ATEST design point calculation improved the engine model 
simulation. The use of the SCIDNT identification code resulted in better 
agreement between the model and the actual engine tested. All the 
weighting factors must be modified to trim the engine model. Modification of 
the airflow weighting factor resulted in closer agreement with the airflow 
measurement. However, it is also apparent from the decreased agreement in 
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net thrust, that all the weighting factors must be examined and modified to 
produce the desired simulation accuracy. 

Further examination including an in depth analysis of measurements 
and associated measurement uncertainties of the adaptive engine simulation 
is necessary to realize the full potential of this code. Initial evaluation of the 
system identification technique has proved it to be a powerful tool in tuning 
engine simulation mathematical models to more accurately predict as-tested 
engine performance. 
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FINAL V,\LL tS 


FTJNCTION \ ,U.L't 1829604 4 



RMS 

FUNCTION 


MEASUREMENT 

ERROR 

BUDGET WEIGHT 


FG 

140.454 

986372. 

1 00000 

FN 

172014 

0.147943E-07 

1.00000 

FRAM 

142.555 

0101609E-07 

1 00000 

pn 

1 37830 

949860 

1 00000 

p-< 

1 11314 

61 9539 

1.00000 

P3 

130021 

845278 

1 00000 

T13 

309582 

479206 

1.00000 

T25 

25 4611 

32413.5 

1.00000 

T3 

41 5665 

86388.7 

1.00000 

W2 

628874 

1977.41 

1 00000 


PARAMETER 

Value 

1 WAC2SCA 

9.6757847E-01 

2 EFF2SCA 

1.0217005E+00 

3 PR2SCA 

1.0135444E-^00 

4 XNC2SCA 

8 7436050E-01 

5 WAC25SCA 

1 0432723E+00 

6 EFF25SCA 

95254391E-0I 

7 PR25SCA 

9 97I3647E-01 

8 .\NC25SCA 

8.4555548t-01 

15 AML\16 

5.0647099E+01 

16 AMIX56 

7.1390942E+02 


Table Al. SCIDNT output with nominal weighting factors. (See reference 4 
for nomenclature) 


FINAL VALUES 

FUNCTION VALUE= 4175071.6 



RMS 

FUNCTION 


MEASLTIEMENT 

ERROR 

BUDGET WEIGHT 


FG 

207.230 

0.214721E+07 

1.00000 

FN 

245.647 

0.301713EH)7 

1.00000 

FRAM 

51.0057 

130079. 

1.00000 

P13 

1.24736 

77.7953 

1.00000 

P23 

1 00294 

50.2948 

1.00000 

P3 

23.6521 

27971.1 

1.00000 

T13 

45.4212 

103154. 

1.00000 

T25 

41.0924 

84429.4 

1.00000 

T3 

76.1612 

290027. 

1 00000 

W2 

2.25832 

0.235001E+07 

O.IOOOOOE-01 


PARAMETER 

VALUE 

1 WAC2SCA 

9.1950560E-01 

2EFF2SCA 

1.0175272E+00 

3 PR2SCA 

1.0l34316E<-00 

4XNC2SCA 

9 8229617E-01 

5 WAC23SCA 

1.0438070E+00 

6 EFF25SCA 

9.2382407E-0I 

7 PR25SCA 

1.1610364E*00 

8 XNC25SCA 

6.1839682E-01 

15 AMIX16 

2.0432941E+01 

16 AMIX56 

30939557E+02 


Table A2. SCIDNT output with modified airflow weighting factor.(See refer¬ 
ence 4 for nomenclature) 
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Figures Al, A2, & A3. Engine test data. 
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Figures AlO & All. Engine test data. 
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Figures A12 & A13. Comparison of actual measurements versus ATEST 

design point calculation. 
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Measured -Adaptive Simulation 


Figures A14 & A15. Comparison of actual measurements versus SCIDNT 

estimated parameters with nominal weighting factors. 
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Figures A16 & A17. Comparison of actual measurements versus SCIDNT 

estimated parameters with a modified airflow 
weighting factor (W = .01) 
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DESIGN AND TEST OF A DIGITAL CROSSBAR SWITCH 
FOR TMS320C40 DIGITAL SIGNAL PROCESSORS 


Bryant E. Sorensen 
Graduate Reseaxch Assistant 
Depeurtment of Electrical and Computer Engineering 
University of Tennessee Space Institute 

ABSTRACT 

A digital crossbar/crosspoint (Xbar) switch has been designed for the commu¬ 
nication ports icommports) of the Texas Instruments TMS320C40 Digital Signad 
Processor (DSP). The Xbar switch allows any one of eight commports attached to 
the switch to communicate with any of the other seven commports. with a majdmum 
of four connections in the switch at a time. The crosspoint switch allows any one 
port to broadcast data to as many of the other seven ports eis desired. This Xbau: 
design weis reviewed, laid out on a printed circuit board, built, and the crossbeu’ part 
of the switch tested during the 1993 AFOSR summer research tour. This report 
discusses the design of the switch and the software routines used to test and control 
the switch. These software routines form the basis for a complete Xbar switch 
control software paclcage. The feasibility of using this switch in a TMS320C40 
parallel processing system is demonstrated by these tests. 




DESIGN AND TEST OF A DIGITAL CROSSBAR SWITCH 
FOR TMS320C40 DIGITAL SIGNAL PROCESSORS 

Bryant E. Sorensen 

I. INTRODUCTION 

The Texas Instruments TMS320C40 Digital Signad Processor (DSP) has six 
high-speed communication ports (commports) on the chip. This presents a great 
advantage in parallel processing systems, ais each processor can transmit and receive 
data from up to six other processors quickly. However, in larger systems that handle 
data acquisition and analysis, it is advantageous to be able to communicate with 
more than six processors for data comparison (such as correlations) or processing 
speed-up. A crossbar switch (Xbar) has been designed which allows any one of 
eight DSP's connected to it to communicate with any one other DSP through a 
commport on each DSP. with a maximum of four connections in the switch at a 
time. .A. crosspoint switch has also been integrated in the system which provides the 
ability for one DSP to broadcast data from its commport to as meiny of the other 
seven ports as desired. This report reviews the design and details the testing of the 
crossbar part of the switch. 

1. Discussion of Problem 

Both the crossbar and crosspoint parts of the switch must be able to connect 
any two of the eight data busses in the system. Each connection must be completed 
in minimum time so the system does not become bottlenecked. 

The crossbar switch must handle the handshaking to chajige the data bus 
direction. Each commport has a bus request line which it asserts when it wants 
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ownership of the bus in order to transmit data. Each port has a bus grant signal 
which it asserts when it has finished transmitting data and is ready to receive data. 
These signals must be handled so the switch can change the bus direction to be 
ready when the commports change their direction cind avoid electrical damage to 
the DSP's. 

The crosspoint switch does not allow direction changes, but instead must 
logically combine the multiple data handshaking signals (the 'data ready’ lines) from 
the receiving ports to create one signal at the broadcasting port. This combination 
must only use those signals from the ports that are receiving data and ignore einy 
ports which are not involved in the broadcast. 

2. Methodology 

The design was divided into five main parts: 

(1) The data bus connection for both the crossbar and crosspoint switch; 

(2) A bus direction handshaking arbiter for the crossbar switch; 

(3) .4 data bus handshaking arbiter for the crosspoint switch; 

(4) .Address decoding for the different parts eind system interconnects; 

(o) A control unit for the entire device with softwaire. 

This report examines the design and metisured performance of the Xbax switch, 
specifically the data bus connections, the bus direction heindshaking arbiter, and 
the control unit with address decoding. Time did not allow for testing of the data 
bus handsheiking arbiter and the crosspoint part of the switch. The results show the 
feasibility of using a crossbar/crosspoint switch in a C40 parallel processing system. 
Data transfer speeds through the switch were also measured. 
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II. THE TI TMS320C40 DIGITAL SIGNAL PROCESSOR’S 
COMMUNICATION PORTS (COMMPORTS) 

1. Communication Port Signals and Their Operation 

The Texas Instrument TMS320C40 (C40) is a 32-bit floating point digital signal 
processor (DSP). It htis the unique feature of six communication ports (commports) 
on the chip. Each commport consists of twelve lines - eight data lines, a data 
strobe (CSTRB), a data acknowledge {CRDY), a bu' request (CREQ), and a bus 
acknowledge (CACK). Each of these lines is bidirectional. The data, CSTRB, and 
CACK lines are driven when the commport is in output mode, with CRDY and 
CREQ being driven in the input mode. The CSTRB cind CRDY lines provide 
common handshaking for the data, with CSTRB being driven low by the output 
port when valid data is present, and CRDY being cisserted low by the input port 
when the data has been received. The CREQ and CACK lines are used to arbitrate 
data bus ownership between two connected ports. The input port will assert the 
CREQ line low when it has a non-empty output buffer and its output has not been 
halted by the DSP’s CPU. The output port will then drive CACK low when it hats 
finished transferring a 32-bit word (four consecutive byte transmissions) or if its 
output buffer is empty. The input port then detects the CACK and removes the 
CREQ signal. This drives the CACK line high at the output port and signals the 
other lines to switch directions. The input port then becomes the output port and 
vice versa. (For further detail on commport operation the reader is directed to [5], 
chapter S). 

On power-up. half of the C40’s commports are automatically configured as 
output ports (ports 0. 1. and 2) and half as input ports (ports 3, 4, and 5) ([5] pp. 
S-14 and S-15). No two input or two output ports should be directly connected as 
this can damage the 040. 
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III. COMMPORT LINE CONNECTIONS USING THE 
LSI LOGIC CROSSBAR SWITCH 

1. Hardware Features of the Crossbar Switch Part 

The actual connections in the crossbeir/crosspoint switch are implemented using 
the LSI Logic L64270 64x64 Crossbar Switch (Xbar) chip. This chip allows ziny one 
of 64 outputs to be connected to any one of 64 inputs, to be set to a constant value 
(high or low), or to be tri-stated. Each output has a 64-to-l mux and registers 
which select the input line to which it ’<vill connect. All outputs can be disabled at 
once also. The chip supports a bidirectional mode, in which each input is hardwired 
to a correspond! output. The Xbar allows for groups of ports to be treated as a 
single entity. These groups must be a power of two in size, ranging from width 1 to 
width 64. This design uses the Xbar in its bidirectional mode. 

To avoid confusion, each input and output on the Xbar will be referred to as 
a 'line' while the combination of a pair in the bidirectional mode will be called a 
■port'. A grouping of these ports will be called a ‘bus'. 

The Xbar chip has 19 control lines. The lines CL0-CI.5 are used in initiedization 
to indicate the bus width to be used. After initialization, they are used to select 
the input lines in the 64-to-l mux. CL6 is asserted in initialization to indicate 
the bidirectional mode, and afterwards is kept low unless the output is to be set 
at a constant or disabled. CI.T is not used in initialization and indicates pipeline 
mode otherwise. It is tied to ground since pipelined mode will not be used in 
this design. The lines REGADR.0-REGADR.5 select the output bus that is being 
configured. REG ADR.6 is asserted to indicate initialization of the Xbar. The 
collections CL0-CL6 and REGADR.0-REGADR.6 will be referred to cis Cl and 
REGADR respectively. BDCYCl is used in the bidirectional mode to indicate 
the configuration cycles (explained below). WE is the control line data strobe. 
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BNKLDI is a latch signal used to load the configuration data from a holding memory 
to the individual output control registers. 0~E is the global output enable signed 
I PI p. 1ST). 

2. Configuration of the Xbar chip 

To initialize the Xbar chip, REGADR.6 is asserted and the bus width to be 
used is specified on CI.0-CI.5 ([2] p. 1S9). Cl.6 is eisserted to indicate bidirectional 
mode is desired. REGADR.0-REGADR.5 and Cl.7 are don't cares during the 
initialization. WE is then pulled low for at least 10 nsec. (For more complete 
timing of the initialization and configuration timing, see [2] p. 192). 

Connections in the bidirectional mode take two cycles to complete. During the 
first cycle. BDCVCT is held high while specifying the output port or bus address on 
REG.A.DR and the input port on Cl. WE is pulled low. then returned high. This 
first cycle disables the output line of the input port so as to not damage the output 
line driver. In the next cycle. BDCYCl is taken low. REGADR and Cl remain 
the same, and WE is pulsed low again. This cycle makes the desired connection 
between the input line of the input port and the output line of the output port. If 
desired. BNKLDI can be held low during both cycles, then pulled high any time 
after completion to latch the configuration. This allows for multiple connections 
specified at the Xbar to be latched in at the same time: for instance, if all four 
connections were to be made in the Xbar before latching them in to start data 
transfers. ([2] pp. 1SS-1S9). 

The .Xbar chips were placed in parallel with respect to the control signals to 
speed up configuratioi; c f the switch. However, since the CRDY and CREQ lines 
run opposite to the other ten lines from the commport. the Xbar chip used to make 
their connections switched the Cl and REGADR from the other Xbar chips used 
to connect the Data. C.-ICA'. and CSTRB lines. 
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IV. BUS ARBITRATION HANDLING OF CREQ AND CACK 

The C40 commports use the CREQ and CACK lines to arbitrate ownership of 
each data bus connection. As explained before, each port starts eis either an output 
port with bus ownership or zis an input port without bus ownership at power-up. 
However, it may be desired to connect two ports which conflict in bus ownership. 
This can create electrical problems ([5] pp. 8-14,8-15) as well as data loss. The 
CREQ/CACK arbitration unit provides both electrical isolation to prevent harm 
to the C40’s as well «is a go-between to assure correct bus direction and ownership. 

The arbiter seeks to assure that the CREQ will be responded to quickly enough 
and that turnaround occurs quickly enough so that the commports, the Xbiirs, and 
the other chips are not damaged in the process. It also works to make changes of 
direction as quickly as possible to achieve maximum data throughput. This is all 
handled in hardware so the control unit does not need to be involved, which would 
slow down the process. The arbiter also handles reconfiguration requests from the 
control unit by assuring that data has stopped flowing before the OK is given to 
reconfigure the system. This should prevent data from being lost. 

1. Bus Transfer Sequence and Timing 

The bus transfer sequence is as follows. The input port of a connected pair of 
commports asserts CREQ when it has data to send. When the output port has 
completed its data transfer or is not sending data, it asserts CACK. The input 
port detects CACK and removes CREQ. This drives CACK back high, at which 
time the commport signals chzuige direction. 

The timing for the tremsfer of data bus ownership (called “token transfer”) is 
given in [5], pp. 14-33 to 14-36. Notice, however, that these are relative timings 
to the internal synchronizers of the C40’s. As far as this author knows, no solid 
absolute timing is known for the bus token transfers. The timings given do indicate 
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some idea of what two commports would do if directly connected and this design 
tries to mimic these as nr irly as possible. The times given will have to be stretched 
in many cases to allow the control unit to respond, however. 

2. State Machine Design for Bus Direction Switching 

A state machine for CREQ and CACK handling was written in Altera Hard¬ 
ware Description Language (AHDL) (see [4]). The corresponding state machine 
diagraun is included in Figure 1. The OUT signal to the state machine is created by 
decoding the address lines of the output port being connected. This address comes 
from the three upper bits of Cl at the Xbar chips. Similarly, the IN signal comes 
from the three upper bits of REGADR. The ‘O’ and ‘I’ inputs are data lines which 
give the output and input ports’ status, i.e., if the desired output port is already 
in the output state, O is asserted and likewise for ‘I’. These status bits are only 
used on the initial configuration sequence since the arbitration state machines keep 
track of each port’s state afterwards. CONFR denotes ‘configuration request’ and 
is 2 isserted to request a reconnection during operation (not at initial configuration 
of the switch). BC is the ‘broadcast’ signal, which is asserted when CREQ zind 
CACK eire to be ignored ^lfter the connection is made (crosspoint switch). LTCH 
is a latching signal which is the same as BNKLDI at the Xbars. It is handshaking 
from the control unit used to show that a connection heis been completed. 

CREQ and CACK are the commport signals into the arbitration unit, while 
CRX and CAX are the request and acknowledge signals fed to the Xbeu’s, re¬ 
spectively. The CREQ, CACK, CRX, and CAX signals all go through tri-state 
buffers from this chip so that they may be turned on and off. They are all pulled up 
through resistors to provide a default high state when disabled. These four outputs 
and their enabling signals are eisserted as shown in Table 1. By separating these 
signals and tri-stating them when necessary, the state machine can respond more 
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quickly to CREQ 2 uid CACK and provide electrical isolation. 



FIGURE 1 - STATE MACHINE DIAGRAM FOR CREQ/CACK ARBITRATION 
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Table 1 - Output Table for CREQfCACK State Machine 


CREQ 


CACK 


STATE CREQ ENA CACK ENA 


sO 

si 

s2 

s3 

s4 

s5 

s6 

s7 

s8 

s9 

slO 

sll 

sl2 

sl3 


1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 


0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

0 

1 

1 


0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 


CRX CAX 
CRX ENA CAX ENA 
10 10 
10 10 
10 11 
10 11 
10 0 1 

10 0 1 

10 10 
1110 
0 110 
1110 
10 10 
10 10 
10 10 
10 10 


OK is a handshaJcing signed to the control unit indicating that the commports 
eire ready to be connected as feu: eis the CREQ and CACK signeds are concerned. 
The OK signal is only asserted in State 1 and State 6, when the state meichine is 
looking for the LTCH signal indicating that the connection is made by the control 
unit. The control unit will be able to read these OK lines to determine when the 
data has stopped flowing and when the ports are in the desired states. 

Each port heis its own state machine to handle its CREQ and CACK signzds. 
All eight state machines are then connected with proper address decoding going to 
the OUT and IN state machine lines. 

4. Software Control of the Arbitration Unit 

A connection on this chip is configured at the exact same time as the Xbars. 
This is made possible by the overlap in addresses at the decoder (see Sec. V). 
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The exact seime addresses are used and the same WE, BDCYCL, and BNKLDI 


signals are used. The WE a-t the arbitration chip enables the input and output 
port addresses at the D flip-flops so that the chip may have time to make necessary 
changes before allowing the connection to be completed. The control unit may make 
multiple connections without receiving the necessary OK signals to complete them, 
but none of these connection will be latched into the Xbar and this arbitration chip 
by BNKLDI going high until all the OK signals the controller is looking for have 
been received. The control unit will assert a ‘read OK’ (RDOK) signal, and then 
will have access to the OK’s from each port on the upper byte of its data bus. The 
BRDCST signal may also be cisserted if broadcast mode is desired. 

When a reconfiguration is desired, the CONFR line is asserted with the ad¬ 
dresses of the two ports involved being latched in as well. The broadcast status 
is also latched in at the BC line. The control unit then asserts RDOK. When the 
ports have stopped transferring data and are going the correct direction (i.e., output 
or input state), the two ports specified will have their OK signals asserted. The 
control unit then reads ail the OK lines, checks for the two from the ports to be 
reconfigured, then makes the changes at the Xbeirs when both OK’s are received. 
This assures that no data will be lost and that improper changes are not made. 

V. ADDRESS DECODING FOR CROSSBAR SWITCH 

Address lines from the control unit are used to address the Xbax chips, the 
CRDY arbiter, the CREQfCACK handler, the CREQ buffer, and provide other 
control signals to the CRDY arbiter, CREQfCACK handler, and the Xbeu: chips. 
The design uses address lines instead of data lines to provide these control signals 
because the address lines are 2 Llways in a determined state and not tri-stated sis the 
data lines are between data transfers. 

Nine address lines (A0-A8) from the control unit are used both to create control 
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signzds and for decoding to access the diiferent parts of the switch. AO is used to 
create the BDCYCL signal to the Xbcirs and A1 becomes BNKLDI. A2 is used for 
CONFR and A3 for BRDCST for the CREQjCACK Eirbitration unit. A4, A5, and 
A6 are used for ADRO, ADRl, and ADR2, respectively, by the CRDY arbiter (this 
report does not cover the CRDY arbiter). A7 and A8 axe aO and al. respectively, to 
the address decoder explziined below. The address decoder zdso uses A2 (CONFR). 
The decoder uses the lOE signal as the data strobe, and the W line as R/W (read 
asserted high, write asserted low) signal. 

The address decoder needs to provide 6 signals. They are RDCRQ (read the 
CREQ lines), WER (write to CRDY arbiter for broadcast setup), WEQ (write 
to CREQICACK handler), RDOK (read the OK signals from the CREQ/CACK 
handler), OEX (enable/disable the CRDY' Xbax chip), and WEX (write to the 
Xbar chips). The decoding for these signals and the memory map of the Xbar 
switch are shown in the figures below. 


PIN 

2 

CONFR 

INPUT 

PIN 

5 

AO 

INPUT 

PIN 

7 

RW 

INPUT 

PIN 

9 

STRB 

INPUT 

PIN 

11 

Al 

INPUT 

PIN 

15 

SET 

CMBFBK 

PIN 

16 

RST 

CMBFBK 

PIN 

17 

WEX 

COMBINATORIAL 

PIN 

18 

OEX 

CMBFBK 

PIN 

19 

ROOK 

COMBINATORIAL 

PIN 

20 

WEQ 

COMBINATORIAL 

PIN 

21 

WER 

COMBINATORIAL 

PIN 

22 

RDCRQ 

COMBINATORIAL 

; logic equations 


EQUATIONS 



/WEX 

- /RW 

* /STRB 

• /Al • /CONFR 

/WER 

- /RW 

• /STRB 

* /AO • Al 

/WEQ 

- /RW 

• /STRB 

* AO • /Al 


+ RW * 

/STRB • 

/AO • /Al 

ROOK 

- RW » 

/STRB • 

/AO • Al 

SET - 

■ /RW » 

/STRB • 

AO * Al 

RST - RW • 

/STRB ♦ 

AO ’Al 

OEX - /SET 

• RST + 

/SET » OEX 

/RDCRQ - RW 

• /STRB 

• AO * /Al 


; configure request 
; low address line 
; Read/Write line 
; Data strobe 
; high address line 
; set for OEX 
; reset for OEX 
; write to Xbars 

; enable/disable Xbar for /CRDY,/CREQ 
; read the OK lines at /CREQ,/CACK chip 
; write to /CREQ,/CACK chip 
; write to /CRDY chip 
; read the /CREQ lines 


FIGURE 2 - .4DDRESS DECODING EQUATIONS FOR XBAR SWITCH 
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FIGURE 3 - MEMORY MAP OF XBAR SWITCH IN CONTROL UNIT 
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VI. software routines for control 

The crossbar switch was Izdd out on a VME-size board, ordered, and put 
together this summer. The design of the Xbar switch was edso reviewed, although 
the main pzu't of the design took place before this AFOSR research tour. The 
board was then tested with the five test software routines given in the following 
pages. These algorithms were tested using a Spectrum Signzd Processing PC-host 
C40 board with a DSP'Link interface (for detziiled information on DSP^Link, the 
reader is referred to [3]) and the accompanying debugger. DSP'Link is a hardware 
interface protocol which can expand memory on the Spectrum board and/or provide 
communication between other DSP elements. DSP'Link resides in the C40’s global 
memory starting at address 9000 8000 Hex, as shown in Figure 4 of the preceding 
section, and provides 16 data lines. 

These routines show the basic commands necessairy to configure the XbM switch 
as well as check the data throughput and integrity through the Xbar. These routines 
are merely rough sketches of the necessary parts of control code for the Xbar and 
need to be integrated into a more complete software package to be useful in a parallel 
processing systems. Each routine hats comments about its use and the part of the 
Xbair it tests. 


The meaning of each configuration word is given in the template below. 


OSQ/CSK cMp 


12 

11 

10 




IS 


02 

01 

00 




08 

Xbars 

RG6 

RGS 

ROM 

RGS 

RGZ 

RGl 

RGO 


CL6 

CLS 

CLH 

CLS 

CL2 

CLI 

CLO 


CR09 chip 









U7 

U6 

US 


US 

U2 

U1 

UO 

Raads 

CK7 

OKS 

OKS 

OKH 

OKS 

0K2 

OKI 

OKS 

CR7 

CM 

CHS 

cm 

oa 

CR2 

CRt 

CRO 


RGk ■ REGAORx 

IS ■ BTAT Onput Port 

OS ■ OSTAT <Outptn Port Si««u>> 

Uk • Data Unas - Indicatsa Port n will rocaiua broadcast 
OKs • OK Unas to Raad tnm QSBrCiS! dUs 
Otx • QHCQx for Raads from OKBSufrar Chip 

FIGURE 4 - CONFIGURATION WORDS TEMPLATE 
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?R0G1.ASM 

ta« connaccion at two co-ioipocts tSrouqn c.-.a XEar 
. data 

XBARI word 9000800ah . addrass za tn 2 .cia.:.:o xdar cnips 

:jkTCH word 9Q00803:h ; aodrass za latcn ra coar.ecrrras 

CRDYEN word 90008190h ; addrasa for enaoixa^ /C?rf ;t=ar 

RDOK .word 90008100h ; addraas to raaa cde OK '.zr.as 

OYCIXl .word 90008081fi ; lat oycia of conaecfion wtita, non-broadcaac 

3ASEO word OOOffSOOh ; baaa addraaa far acarr of coda 

3LOBINT .word OlOfSOlOh ; giobai mamory control word 

OONNEC7 .word QOOOeiZlh . donnactxon to ba aiaca 

doen^ortO cti .word Q0100040h; Port 0 control 

coftn oortO rap .word 0OXQ0O4Ih; Input 

conn oortO otp .word 00100042h: output 

conn^__poct3 ctl .word QOIOOOdOh; Port 3 control 

coan^ort3 lap .word OOlOOOOXh.* xaput 

conn_pott3_atp .word 00X00072h; output 


. Z9XZ 


LDP 

8BASECD,0P 

LOA 

8BASECD,SP 

LDHI 

OOXOh,ARO 

LDI 

SGLOaiNT,RO 

STI 

RO,'ARC 

create 

the OK maak 

•JDI 

3CONNECT,R3 

IDl 

-4, R4 

LDt 

-12, R5 

Lwl 

0X00h,ARX 

LSI 

070h.B6 

Lo: 

070CQh,a7 

AND 

R3,R6,R0 

AND 

R3,a7,RX 


R4, RO 

ISH 

R5,RX 


RO,ARX,RO 

L3H 

RX,ARX,RX 

OR 

R0,R1,R6 

load control ragxatera for 

lOA 

aXBARI,AR0 

IDA 

9ROOK,AR2 

IDA 

8CRDYEN,AR3 

IDA 

9UiTCH,AR4 

IDA 

9CfCl,El,ARS 

' initialize the Xbar chipa 

STI 

RO, ‘ARS 

ID I 

9E8Eh,Ra 

STI 

RO, *AR0 

make one~one connectiona 

371 

R3, 'ARS—(X) 

STI 

R3, 'ARSw-dRlI 

:escok: 


AND 

R6,*AR2,RX 

SUB! 

R6,RX 

3NZ 

teatok 

IDI 

•AR4,R4 

IDA 

acomm_pottO_ctl,ARO 

IDA 

3comm_port0_ocp,ARX 

:-DA 

3 c otiin_po r 13 _c 1 1, AR2 

DDA 

3com_poct3 inp,AR3 

'JDHI 

5555h,RX 

ZR 

;5a5h,RX 

AND 

ODFFFh,3T 

IDI 

0100h,RO 

ID I 

XE00h,R2 

D’ST: TST3 'ARO, RO 

3NZ 

CUT 

5t; 

Rla ’ARl 

:N: and 

•AR2,R2aR3 

3Z 

:n 


•AR3,R3 

5'JBI 

Rla R3,R4 

3Z 

CUT 


.' aat atacK and data poxatera 

.’ Xoad addraaa for global m.nory xntarfaca rag 
; Xoad global .namocy xntarfaca control word 
; aat tha global mm xntarf rag 

; aat up connact addraaa 
; ahxft countar, gat 2nd nxbbla 
; anift countar, gat 4th nxsbla 
; numbar to ba ahxfted 
; maaK for 2nd nxbbla 

; maatc for 4th nxbbla 

: gat 2nd nxbbla 
; gat 4tn nibble 
; anift down to rignt 4 bita 
; ahxft down rxgnt X2 bxta 
; ahxft X laft by RO bxta (mx 7) 

; ahxft 1 left by R1 bita l.mx 7) 

; OR reaulta and atora xn .R6 
Xbar awitch 

; Xbar initialization addraaa 
; addraaa to read OK'a 
; addraaa to enable/diaabla /CRCY xbar 
latch aignal for connactxona 
; Xbar addraaa for Xat cycle of configure write 

; dunay write to enable /CRSY Xbar 
: control word indicating f-wide porta at Xbar 
; initxallza the Xbara 

; nuilce the firat connection cycle 
; aecond connection cycle keeping BNKLDI low 

; read OK Ixnea (AR2->OK addraaa,RE^OK aaaX) 

; check if OK'a are aaaerted 
; keep checking if not both aaaerted 
: if both OK'a aaaerted, latch xn connection 
; by changing 3NKLOI to nig.n I dummy read! 

; load comm port 0 cont reg addr 
; output FIFO addr 

: load conra port 3 control reg. addreaa 
; input FIFO addreaa 

; Output word aetup 


CPU global interrupt diaaoled 
Maak for bit 8 of CPCR, port 0 
Maak for i.nput level of port 3 
Check if O/P FIFO ia full 

If yea, check again 

write value to O/P, port X 

c-neck input level of port 3 

keep checking if input FIFO ia empty 

Xoad value from I/P, port 3 

Check that the output ana input valued are 

if they ate, keep outputting 

otherwiae atop 
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PR0G2.ASH 

Test o£ coflwport daca sp^^<i by polling tnathod 


COW PORTO CTI. .word 00100040n 
^CMM”'?ORTO"daT .word 00100042h 


CC>W^?CRT3”0At 


.aotd Q0100071h 
.word 12S00000 
COMM_?ORTO_CTL 
3C0MM PORTO^CTL.ARl 


port 0 control 
pore 0 oucpuc riro 
port 3 input rCtO 
nuxBcor or woras to sane 


; bytes eacn) 



IDA 

3COMM PORTO 3AT.AR2 


IDA 

3C0MM PORT3_DAT.AR3 


ISA 

aLOOP COUNT.AR4 


Ajro3 

OEFh,'ABl 

.R9 set up port 0 rcc oucpuc 


st: 

R9.-ARl 



tOHI 

55AAh.R0 

; word to be output 


OR 

55AAh.R0 



STI 

RO.•AR2 

; output 3 words to fill port 0'a output FIFO 


STI 

RO.•AR2 



st: 

RO. •AR2 



STI 

RO.‘ARa 



STI 

RO.•AR2 



st: 

RO.•AR2 



STI 

RO. •AR2 



STI 

RO.•AR2 


4AIT1 

ISH 

-4,*ARI,R1 ; vait for FtrO to o« 1/2 full 

LOOP 

3C 

WAITI 



STI 

RO, •AR2 

; Store 4 more words 


STI 

RO. 'AR2 



SIX 

RO. -ARZ 



st: 

RO. 'ARZ 



IDl 

•AR3.Ri 

read the transferred words frer*. input rirO 


’.n: 

V^R3.R1 

of port 3 


030 

AR4 loop 

; decrement the count 


IDI 

•AR3,r; 



101 

•AR3.R1 



ISH 

-a,-ARl, 

R1 ; check output FIFO level agai.h 

iTCP 

3R 

STOP 

: enaless loop when done (breakpt) 

?ROG3 

.ASiM 



Check 

State 

mac.hines 

ability to change cotwwsort oonnectith iirecticn 

AECONN 

. word 

00002 

060h ; switch direction of ports i.-.put->oucput. 


SAME AS PROGl UP TO HERE 
10r 3R£CCNN,R7 
lOA 3coflW_pocc0^inp,ARl 
lOA 3coinm^ort3_ocp/AR3 
IDI ;,R3 


.*naJce change in direction by outputting to port 3 


load connecti^on change^ OK .r.asR is same 
changing over from output to input 
changing input to output 
token Showing transfer successful 


OUT3: TS7B •.^R2,R0 


SNZ C 
3TI f 
3TI f 
3TI } 
testok: M 
SUB I i 
asz 1 

IDI 

:n0: and 

3Z 
101 
SUBI 
3NZ 1 
STOP: SR 


OUT 3 
R3.*AR3 
R7, ♦ARS—(IJ 
R7,*AR5f+ 
and <AR6.R6.RI 
R6, R1 
testok 
•AR4,R4 
I VARQ.RZ^Rl 

:no 

• ARl, R1 
1,R1 
ERROR 
JR STOP 


ERROR: BR 


Check if o/p buffer full 
if so, check again 

write out to poet 3, initiate cnange 

tnake the first connection cycle 

second connection cycle ■<ceeping 3NKL0I low 

read OK lines <ARS->CK address, R6wOK masK) 

check if OK's ace asserted 

keep checking if not botn asserted 

latch in connection 

check input level of port 3 

keep checking if input fzro is empty 

load value from I/?, port 0 

check for correct value 

go to error if not • 1 

otherwise stop 

reaches here on an error 


PROG4.ASM 

Cherlc for Seconfigur* abiiicy oi state aac.'tines 
CCNFR .woed 9000608Sh reconfigure reoues* 

RECONN word 00002031!! ; reconnection 

coii«n_oortl_otl .word QOlOOOSOh ; commporc I control 

coB»B_poctl_otp .word 00100052h ; portl output FIFO 

create toe OK mask for the 2nd connection 
iOI 3R£CONN,a7 ; set up connect address 

(same procedure as in ?ROGl) 

OR R0.r 3.R6 ; OR results and store in 

Make connection, check the data flow as in PRCGl *•* 


101 3COKFR.AR7 

lOK 3coasa_port0_inp, ARl 

lOA 3ooni«_portl_cci( AR2 

iOA 3coiiiB_pottl_otp, AR3 


address to request reconfigure 
changing over from oucpuc to i.-.put 
goto port 1 
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I —Iw th« caeontlgur* nquaac 
SIX K7,*AJR7 
tMCOk: AMD *kRS,R6,IU 
SUBI M.Rl 
am tMtok 
STX 

arx R7,»juu~(i) 

LSI •AR4,M 
LDI •AXC.RS 


*** ClMck iMif colUMctloo aa in PROGl ••• 


■aka UM taquaac — pace I wluia ba iaie'od 
chock for cha OK to coaa back 
aaa if it ■atcbaa cha aaak 

■aka Cha eonaaction 

caad Co latch in connacclon 
caad (chock OK' a) to 90c latch o«f 


PROGS.ASM 

; Chock acaco aachlna'a ability to cum porea around at initial coim.«i»» 
••• SAMX AS PROQl EXCKPT FOR CONMKCTIOM KORO SHOW! BElSw ^ connaction 

.dAC4 

COIOIKCT .woed OOOOSOCOh ; NOTK THAI BOTH PORTS HAVK TO RX TtJRMED AROOMD 


VII. RESULTS 

The board was first tested the ensure that the parallelization of the crossbar 
chips would work. This proved to be the czise. The three chips can all be pro¬ 
grammed for a certain connection in two or three cycles (two if the connection 
flows through instead of being latched in adter it is written). These cycles must 
be accompanied by the handshaking necessary to ensure the state machines are 
responding, speciflcally by checking the OK lines from the CREQICACK chip. 
With this test, it was tdso verified that the data, CACK, and CSTRB lines 
can flow opposite the CREQ and CRDY lines, thus proving the concept of the 
parallelization of the LSI Logic Crossbar chips, at least in the bidirectionad mode. 

Next, the data flow rate was checked using am efficient polling method with 
the commports. The polling only took place every four words aind is compairable 
to DMA for speed. In multiple tests, the crossbar achieved data throughput on the 
range of 7.2 MBytes/sec to 7.5 Mbytes/sec. This compaires to the approximately 
12 MBytes/sec that similar tests provide on straight commport-to-commport con¬ 
nections with the approximately 1 meter of cable used. The dropoff in speed comes 
from the propagation delay of the signals through the Xbar chips, which cam reach 
25 nsec each way. These data rates are encouraging since they do not even cut in 
half the originail speed of the link in this test. 


8-18 




The state machines for CREQ/CACK heindling were then further tested, 
specifically to check for direction changing between two connected commports. 
Again, the results were favorable. Although the TI data book shows relative timing 
for these bus ownership transfers, it was discovered that much more leeway exists 
for the CREQ and CACK signals that one might expect from studying the data 
book. This is very encouraging as it indicates the ability of slower PLD’s to keep 
up with the C40. The timing for the CRX and CAA'signals which come from the 
CREQtCACK chip and go to the Xbars can be improved, and hopefully with the 
next design revision this will be corrected. 

The next test demonstrated favorably that the Xbar can reconfigure during 
(operation (on-the-fly). Two ports were connected, then one was requested to 
reconnect to another port after turning around. The data was once again able 
to go across the Xbar without corruption after both the original and the second 
connections were made. 

Finally, the state machine demonstrated that it can turn around ports at 
initialization. It was requested to change both an input port to an output, tmd 
an output port to an input. The connection wzis made successfully and data wzis 
again sent uncorrupteci through the crossbar. 

VIII. CONCLUSIONS AND FUTURE WORK 

The crossbar switch has been proven in a limited manner to perform satis¬ 
factorily in allowing arbitrary connections between any two commports attached 
to the switch. Data transfers occur uncorrupted at a maximum tested speed of 
7.5 MBytes/sec for polled commport trcinsfers, eis compared to 12 Mbytes/sec for 
direct conimport connections under the same algorithm. The switch should become 
a useful tool in parallel processing systems using the TI TMS320C40, allowing 
topologies to be changed in the system during processing. 
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Time did not allow for testing of the crosspoint switch part of the design, but 
this is a much simpler part of the switch; this part will be tested soon in the future. 

Much remains to be done in developing the Xbar for further use. The software 
algorithms need to be integrated into a program which can run independently so the 
controlling C40 can automatically monitor the switch and maice the changes. The 
software also needs to be optimized for speed to make the switch even more useful. 
A host PC int«-rface is being developed now so that the switch can be controlled 
in a user-friendly environment. The switch can also be expanded to 16, 32, or 64 
ports. Work has already been done in the design of these larger switches, but the 
concepts proved in this smaller S-port version carry over quite nicely to the larger 
switches. The actual hardware implementation, debugging and checking of Izirger 
switches remains for the future. 
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INCOMPRESSIBLE NAVIER^TOKES FLOW ANALYSIS 
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Michael A. Weaver, MSAE 
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Abstract 

An incompressible Navier-Stokes flow solutirm was obtained for the Arnold Engineering 
Development Center HEAT-Hl Test Unit nozzle coding passage. This flow solution is one component of 
a conjugate heat transfer algorithm develt^ied by ttte author to predict heat load failure in the HEAT-HI 
Test Unit. The HEAT-H1 Test Unit is a high enthalpy, arc-heated, free jet test facility. A backside water- 
cooled nozzle accelerates hi^ temperature air to supersonic parallel flow at the exit. The research flow 
solver INS3D-UP (horn NASA Ames Research Center) was used to solve the steady, incompressible 
Navier-Stokes equations for water flow on the backside of the nozzle. INS3D-UP uses an i^rwinded 
discretization scheme with the pseudo-compressibility solution method. An axisymmetric solution was 
obtained by solving the full three-dimensional system over three azimuthally rotated planes, while cnftxcing 
symmetry boundary cmditions on the two outer planes. The results indicate the presence of recirculating 
flow and a stagnation point in a critical region of the cooled nozzle wall Figures show velocity vectors and 
the pressure distribution for this critical region. The flow solution does not consider two-phase, subooded 
boiling effects, which are expected to be present 
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INCOMPRESSIBLE NAVIER-STOKES FLOW ANALYSIS 
FOR THE ARNOLD ENGINEERING DEVELOPMENT CENTER 
HEAT-Hl TEST UNIT NOZZLE COOLING PASSAGE 

Michael A. Weaver 

IntroductiQQ 

The Arnold Engineering Development Center (AEDQ HEAT-Hl Test Unit is an arc-heated, free 
jet test facility, providing extremely high enthalpy air flow. Flows with enthalpies rangiiig from 2,000 to 
8,500 Btu/lb„ and pressures ranging from 20 to 1 IS atm are routinely {Moduced [1]. Mach numbers 
ranging from 1.8 to 3.5 are achieved with interchangeable nozzles. These nozzles are actively cooled with 
backside coflowing water undergoing subcooled nucleate boiling. 

Proposed future applications for the HEAT-Hl require operating conditions beyond the current 
survival envelope. These applications need pressures of up to 200 atm. while nozzle wall heat load failures 
have been observed at pressures in the range 120 to 130 atm. The ability to identify and predict the heat 
load failure modes of the nozzle wall is crucial to extending operation of the HEAT-Hl beyond its current 
capabilities. 

Heat load in the HEAT-Hl nozzle wall can be categorized as a coupled or conjugate heat transfer 
phenometwa The conjugate heat transfer domain comixises coflowing high temperature air and subcooled 
water coolant on opposite sides of a copper-zirconium, converging/divetging nozzle. Heat transfer ditough 
the nozzle wall is characterized by solid body conduction with convection boundary condititms along the air 
side arxl water side of the nozzle wall. The air side heat transfer is characterized by forced convection with 
a turbulent boundary layer. The water side heat transfer is characterized by forced convection, subcooled, 
nucleate boiling. Convective heat transfer coefficients on each side of the nozzle wall are functions of die 
wall temperature and the respective flow properties, thus coupling the three regions of the domaia 

The AEDC operating contraaor for the HEAT-Hl (Calspan Corporation/AEDC Operatitms) has 
conducted a heat load failure analysis by assuming one-dimensional conjugate heat transfer through the 
nozzle wall [2, 3]. This approach gives good results for the thin-walled region of the nozzle, but failure 
analysis for the complete nozzle must coisider multidimensional conjugate heat transfer. With sponscMship 
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from the Air Force Office of Scientific Research (AFOSR), the author pursued such a multidimensianal 
approach while paiticipating in the 1992 AFOSR Summer Research Program [4]. This multidimensianal 
heat load failure analysis subsequently received additional funding through the 1993 AFOSR Summer 
Research Extension Program, and the 1993 AFOSR Summer Re^iearch Program. 

The multidimensicnal approach requires numerical models (based on eitter analytic medrods or 
engineering mathematical correlations) for the unsteady solid body heat transfer in the copper-tdictxiium 
nozzle, the flow properties and heat transfer coefficients on the air side of die nozzle, and die flow 
properties and heat transfer coefficients on the water side of the nozzle. During the 1992 AFOSR Summer 
Research Program, the author assembled the numerical models, developed a flilly-ooupled soludon 
algorithm, and demonstrated the feasibility of the solution scheme. The 1993 AFOSR Summer Research 
Extension Program provided nine months of funding for refining the numerical models and scdution 
algorithm. 

The 1993 AFOSR Suimner Research Program provides an additional three months of funding for 
improvement of the flow properties model on the wa^r side of the nozzle. This effort seeks to obtain a 
complete incompressible, Navier-Stokes flow solution for the water side of the HEAT-Hl nozzle. 


Discussion of 




In the HEAT-Hl Test Unit, tangentially injected air flows between an anode and a cadxxfe for 
heating by an electric arc. The tangentially injected air creates swiriing flow v^ch serves to keep the 
electric arc from making contact with the wall in the regicxi between the anode and cathode. lliB high 
temperature, subsonic air passes from the cathode, through a conical contraction, to die nozzle inlet, then 
expands through the nozzle to parallel, supersonic flow at the nozzle exit The HEAT-Hl arc-heater 
configuradon is shown schemadcaily in Fig. 1. 

This and the previous studies focus on the HEAT-Hl Mach 1.8 nozzle. This axisymmetiic noole 
is 2.3 inches long, with a throat diameter of 0.9 inches. The material composidon is cof^r-zirconium. 
The nozzle and surrounding cooling jacket are approximated in Fig. 2. Water flows through the coding 
passage between the nozzle and cooling jacket providing backside cooling for the nozzle wall. Away fitxn 
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the end regions, heat transfer through the nozzle wall is approximately one-dimensianaL Near the flanged 
ends of the nozzle, the solid body heat transfer becomes axisymmetric in nature. 


ANOM (-*-) 


ILKTWAK 


UTHODE (-) 



TANGENTUU. INJCaiON Of AIR 


SfCnOi A4 


Rg. 1. Schematic repreaentationoflheHEAT*Hl are-heated wind tunaeL 



coouNa jAcnrr 


flOZZLE 
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noramifonn, dissociated, equilibrium air, with decreasing stagnation temperature ai^iroaching die wall 
Transport properties for the hi^ temperature, high pressure air are strong functions of temperature and 
weak huictions of pressure. 

Water flow through the cooling passage is assumed steady, incompressible, viscous, and turbulent 
The water enters and exits the active cooling region in the subcooled state. Subcooled nucleate boiling is 
assumed to occur at the nozzle wall, while the outer coding jacket wall is assumed adiabatic. Transport 
properties for the subcooied water are functions of temperature and pressure. Due to high flow oonvectian 
rates, the momentum properties of the flow (density, velocity, and pressure) and the flow temperature 
distribution are decoupled. 

Method of Analysis 

In the original approach [4], the flow properties of water through the armular cooling passage were 
determined by application of the steady-state conservation laws to a one-dimensional control volume with 
area change. An elemental control volume was assumed with the known quantities, 

Aj, = inlet area. 

- wetted surface area, 

= equivalent hydraulic diameter. 
m = mass flow rate. 


= inlet static pressure, 
p. = absolute viscosity, 

p = amstant mass density. 

Then from conservation of mass in the coitrol volume, the inlet and exit bulk velocities were 


,, m . m 

^OMlP 


The expression for conservation of momentum in the control volume was given by 
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Here the turbulent wall shear stress x was detennined with Prandtl's universal law of friction for smooth 
inpes [S], and the known values D^, V^, V^, ami p. This left exit pressure as the only uidmown m the 
momentum equation. 

The improved approach uses the incompressible fonn of the full three-dimensional Navier-Stokes 
equations to obtain turbulent, steady-state velocity and pressure distiibutiais of water through the coding 
passage. The governing equations are formulated in the research flow solver INS3D-UP from NASA 
Ames Research Center [6]. This flow solver uses the pseudo-compressibility method, first introduced by 
Chorin [7], and successfully applied by Kwak et al. [8], Rogers and Kwak [9], and others. 

Briefly summarizing this approach [10], the continuity equaticxi for the incompressiUe Navier- 
Stokes equations reduces to 

VV=0. 

The steady-state momentum conservation equation, without body forces, beccunes 


v(v®v)=lv.(-p/+f), 


where, V 
P 
P 
1 

X 


= Cartesian velocity vector (u,v,w) = (V,,V' 2 ,V,), 
= constam mass density, 

= static pressure, 

= unit tensor. 


= shear stress tensor = x = p 


——- 


dx. dx 


jj 


p = absolute viscosity. 

The pseudo-compressibility method then replaces the continuity equation with the time-dependent 

^|*pv.v=o. 

The parameter ^ represents a pseudo-speed of sound for the transformed system, with pseudo-pressure 
waves propagating at finite speed. Thereby, the elliptic governing equations have taken on a hyperbolic 
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character. The pseudo-speed of sound ^ must be adjusted for optimum oonveigoice to the steady-state 
solution. The apparent transient behavior (tf the transfbnned system has no physical meaning until the 
divergenoe-hee condition, and thus conveigence, is reached. As implemented in INS3I>-UP, an t^rwind 
differencing scheme based on flua-diffeienoe splitting is used to compute the convective tenns, and second- 
order central differencing is used for the viscous terms [6, 9]. Also, the governing equations are 
nondimensionalized, and solved in a generalized curvilinear cooniinate system. 

An axisymmetric flow solution is obtained from INS3D-UP by solving the full three-dimensional 
problem on three azimuthally rotated planes. Axial symmetry boundary conditions are enftMced on the 
lateral walls of the resulting wedge-shaped configuration, thus giving an axisymmetric sohrtian on the 
center plane. For the HEAT-Hl problem, each of the two outer {danes are rotated by 1** from die center 
plane. 

The axisymmetric computational domain for the HEAT-Hl Mach 1.8 nozzle cooling passage is 



Hg. 3. Computational domain for the HEAT-Hl Mach 1.8 nozzle coding passage. 






shown in Fig. 3. No-slip velocity condititHis are enforced at all solid boundaiies. Flow enters the inflow 
txMindary with a ^)ecified total pressure p„, and exits the outflow boundary with a q)ecified static pressure 
p^. In the course of obtaining a flow solution, the pressure differential (p, - p^ must be adjusted to 
return the blown system mass flow rate m. Since experimental values for p, and m are available, the 
pressure differential should be adjusted exclusively by varying p^. 

Experimentally measured values for p^ and m are shown in Table 1 for four heat load failure test 
cases culled fn»n the HEAT-Hl Test Unit run history with the Mach 1.8 nozzle [2]. All four test cases 
have the identical value of 1,000 psia for the water inlet total pressure. Cases I and n have the identical 
value of 11.54 Ib^s for water mass flow rate. The water mass flow rates for Cases III and IV are nearly 
equivalent In the present arudysis, results will be shown only for Case I inlet total pressure arul mass flow 
rate conditions (also equivalent to Case II). 


Table 1. Heat Load Structural Failure Test Cases fortl^ HEAT-Hl Mach 1.8 Nozzle 


Test Case Number 

I 

n 

m 

IV 

Water inlet total pressure (psia) 

l.(X)0 

1,000 

1,000 

1,000 

Water mass flow rate Gb,ys) 

11.54 

11.54 

7.09 

7.063 

Water inlet temperature (R) 

557 

552 

520 

520 

Water temperature rise (R) 

25 

27 

34 

36 

Air total pressure (atm) 

126.5 

137 

104.4 

94.3 

Air total temperature (R) 

9.000 

9,432 

8,280 

9,180 

Air total enthalpy (Btu/lb^ 

3.480 

3,750 

3,062 

3.644 

Survival/Failure 

Survived 

Failed 

Survived 

Failed 


The computational grid for the HEAT-Hl cooling passage contains 243 nodes tangential to flie 
mean-flow direction, 61 iKxies normal to the mean-flow direction, and two azimuthal planes rotated 1° fiom 
the central axisynunetric plane. In all, the computational grid contains 44,469 ixxle points. Hyperbolic 
tangent node clustering is used in the normal mean-flow direction to capture boundary layer effects. In the 

tangential mean-flow direction, node points are tightly clustered in regions of anticipated large flow 
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gradients. For illustration. Hg. 4 shows a portion of the computational grid surrounding the first 
(upstream) Qow turn. Tight clustering can be observed near the solid walls, and in the vicinity of the flow 
constrictioa For clarity, node points m the normal mean-flow diiectian have been thinned by two-thiids. 



Fig. 4. Detail of the computatiotuil grid for the HEAT-H1 Mach 1.8 nozzle cooling passage. 

Heat load deformation and failure in the HEAT-H 1 nozzles has been observed in the region 
between the nozzle inlet and the nozzle throat. Nozzle wall failure has not been observed downstream of 
the nozzle throat. For this reason, the grid detail shown in Fig. 4, is also the area of primary interest for 
obtaining a flow solution. A recirculatitMi region, or a nozzle-side (hot side) stagnation point in this area of 
the flow field might be indicative of observed heat load failures. 
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Bssults 


INS3D-UP was used to obtain a converged flow solution for the HEAT-Hl Mach 1.8 nozzle 
cooling passage with Case I inlet conditions. For Case I. the inlet total pressure p, was 1.000 psia, and the 
mass flow rate m was 11.54 Ib^s. An outflow static pressure of 832.5 psia was required to maintain 
the value of m to within 0.2%. The solution used a constant mass density p of 0.035727 lb,gAn^, and a 
constant absolute viscosity p of 3.4422 x 10 lb,^ s. A value of 10’ for the parameter p (here, the 
square of the nondimensianal pseudo-speed of sound) gave a stable, well-behaved convergence histCHy for 
the solution. Convergence required 2,8(X) time steps, with a time step size of 10 (nondimensianal units). 
Unsuccessful anempts at using larger time step sizes produced quick solution divergence. 

Results from the flow solution indicate that large recirculation regions exist in the vicinity of the 
first flow turn. Velocity vectors scaled by magnitude in Fig. 5 show aorfib^ratinn as the flow ap pmarh^t 



Fig. 5. Velocity vectors in the HEAT-Hl Mach 1.8 nozzle cot^g passage (Case I)- 
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the constriction at the beginning of the turn. As the Qow negotiates the hairpin turn, two separation regions 
develop. The first separation region hugs the inside comer of the turn, and contributes to the elongated 
recirculation pattern carried partially downstream. The second separation region occurs on the outside 
edge of the turn, and leads to the large circular recirculation pattern. 

Of particular interest, a stagnation point lies on the nozzle-side (hot side) of the cooling passage 
wall at the edge of the second separation region. A close-up view of this s^nd region is shown in Hg. 6 
with constant length vectors. The stagnation point appears as a confluence of velocity vectors on the lower 
wall. A complementary stagnation point exists on the upper, adiabatic wall. 



Fig. 6. Recirculation patterns in the HEAT-Hl Mach 1.8 nozzle cooling passage (Case I)- 
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Pressure contours for the flow field correspond closely to the velocity solution. In Hg. 7. die 
stagnation region on the nozzle-side wail appears as a high pressure area of approximately 1,000 psia. 
Pressure values then drop to approximately 77S psia at the throat, as the flow accelerates downstream fiom 
the flow turn. The high pressure gradients across the constriction correspond to slope discontinuities in the 
computational grid. Also, the slope discontinuity at the sharp comer of the flow turn is most likely 
responsible for a small region of nonphysically negative static pressure (not apparent in the figure) near the 
comer. Any future work on flow analysis for the HEAT-Hl nozzle cooling passage must include 
improvement to the computational grid. 


Increment = 2S psia 



Fig. 7. Static pressure distribution in the HEAT-HI Mach 1.8 nozzle cooling passage (Case I). 
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CoaclusioD 


The existence of a rpciiculatian region and a stagnadtxi point on the hot wall reduces the efficiency 
of heat transfer fiiMn the wall. In the HEAT-Hl nozzle coding passage, heat is swqrt away fiom the 
nozzle wall by continuously flowing cool water. If the water flow is allowed to recirculate in any single 
region of the wall, the local water temperature will rise, thus removing less heat from that wall legioa 
Wall failure might occur when more heat is added to the air-side of the nozzle wall, fltan can be removed 
from the coding passage side of the nozzle wall. The low or vanishing water velocities in the vicinity of 
the stagnation point will similatly reduce the amount of heat removed from that region of the wall. For 
these reasons, the cooUng passage geometry should be redesigned to eliminate recirculation regions and 
stagnation points from the hot wall. 

The HEAT-Hl nozzle cooling passage flow soluticm is one component of a conjugate hea transfer 
algorithm developed by the author to predict hea load failure in the HEAT-Hl Test Unit Bulk flow 
velocities from the flow solution will be used with the energy conservation equation to obtain a bulk 
temperature distribution in the cooUng passage. The energy solution does not need to be obtained 
simultaneously with the flow solution, because the inoHnpressibility and high convection rates of the flow 
combine to decouple the continuity equation and momentum equations from tlx energy equation. Once 
determined, the bulk temperatures, aong with bulk velocities and wall static pressures from the flow 
solution, will be used with Shah's correlation for subcooled nucleate boiling heat transfer [11, 12]. The 
subcooled nucleate boiling correlation will predia heat transfer coefficients aong the hot wall of the 
HEAT-Hl nozzle cooling passage, for use in finding the axisymmetiic conduction solution for the nozzle 
wall. 

It is important to note that the flow solution obtained with INS3D-UP does txn consider two-phase 
flow effects. Two-phase flow arises from the presence of subcooled nucleate bdling on the hot wall of the 
nozzle cooling passage. Water vapor bubbles wodd act to decrease the effective flow area, and flius mi^ 
lead to higher flow velocities than would be expected in single-phase flow. Shah's correlation attempts to 
partially account for this effect 
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Abstract 

A combination experimental/numerical investigation of flow separation for a 
static airfoil has been conducted using Higher-Order Spectral Analysis (HOS) 
and the Volterra functional series approach. Input and output time series 

were recorded downstream of the airfoil at certain points in the flow to 
investigate and quantize the natural energy transfer mechanisms present. 
Time scries were recorded for two angles of attack, and hot-film anemometers 
were used to record input and output flow velocities downstream from the 
airfoil at various fractions of the chord length downstream. Linear and 

quadratic transfer functions were calculated for these trajectories. These 

functions were then examined for evidence of natural energy transfer 
mechanisms in the flow. 

Following this analysis of the unforced flow, a study was conducted using 
pulsed air as a mechanism for control of flow separation. Linear and quadratic 
transfer functions were calculated for these test cases, and insights were 

obtained on the mechanism by which the angle of incipient flow separation is 
increased. 
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Introduction 


A large body of previous work has been done on controlling growth of a two- 
dimensional mixing layer. In particular, the work of Oster and Wygnanski 
[OSTE] examined the use of a small flap to periodically force the flow and alter 
the growth of the flow instability downstream of the plate. Ho and Huerre [HO] 
discussed the use of acoustic as well as mechanical forcing for control of the 
mixing layer. More relevant to the problem of flow separation downstream of 
a static airfoil, Lovato [LOVA] investigated active control of the shear layer on 

j an airfoil using external acoustic forcing and internal-tangential pulsed air 

control. Their results showed that the airfoil shear layer responds to the 
active forcing akin to a free shear layer. 

j 

i 

I An extensive amount of research has been performed concerning the analysis 

of fluid flows using Higher-Order Spectral analysis (HOS). Of particular 
interest to the research discussed below is the work of Miksad, Jones et. al 
(MIKS). They studied the transition of a laminar two-dimensional wake to 
establish the role of amplitude and phase modulations in the spectral- 
broadening and energy-redistribution process. Later work by Ritz, Powers 
and Miksad et. al. (RITZ) concerned the spectral dynamics of transitioning flow 
in the wake of a flat plate. For this work, quadratic transfer functions from a 
Volterra functional approximation were calculated to examine the nonlinear 
spectral dynamics occurring between two sensor points. Finally the work by 
Hajj, Miksad and Powers [HAJJ] involving quantification of the possible 
nonlinear and parametric resonance mechanisms associated with the 
transition to turbulence in plane mixing layers was most influential in the 
following research. In this paper, HOS analysis techniques were used to 
investigate both nonlinear and parametric transfers of energy to lower 
frequencies. 

The work performed by the principal investigator this summer of 1993 builds 
on the above work by examining the flowfleld downstream from a static airfoil 
fixed at an angle of attack. Many comparisons have been formed with the flow 
downstream of an airfoil and both a mixing and shear layer [LOVA]. Hot film 
anemometry was used to record two time series simultaneously downstream 
from the airfoil at various relative positions to the airfoil. Quadratic transfer 
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functions were then calculated for the time series data using the approach 
developed by [RITZ], and these were analyzed to confirm fundamental flow 

frequencies and basic energy transfer mechanisms in the flow. In addition, a 

brief study was done on the energy transfer mechanisms involved in the 

dynamically pitching airfoil. A similar analysis was conducted as for the static 
case, with the hot film anemometers held at a fixed position relative to the 
pitching airfoil. This work will lead to future study of this important problem. 

Experimental Set-up 

The hot film anemometry experiments were conducted in the Frank J. Seiler 
Laboratory's open return, low speed wind tunnel at the U.S. Air Force 
Academy. The .91m x .91 m test section is designed for use with flow 
visualization and flow sensor measurements, and has free stream turbulence 
intensity levels below 0.5% A NACA-0015 airfoil with a 15.2 cm chord and a 
61.0 cm span was used for all tests. Chord Reynolds numbers were at 24,000 for 
the hot-film data. The two time series were recorded using dual 60 mil hot- 
film probes in conjunction with a TSI Model 1054A anemometer. The overheat 
ratio was set at 30%. All data acquisition utilized a MassComp MC-5500 micro¬ 
computer. The flow was sampled at a frequency of 80 Hz for use in the HOS and 
transfer function analysis. The signals were all filtered with a low-pass filter 
at 300 Hz. 

Flow modification was accomplished individually by internal tangential- 
pulsed air blowing. Pulsed air was directed through a downstream facing slot 
located on the upper airfoil surface. The slot was 1.6 mm wide and spanned the 
middle 15.0 cm of the airfoil, or 25% of the entire span, located at the leading 
edge of the airfoil. The system employed two solenoid valves capable of 
pulsing between zero and 45 Hz. The valves were controlled through a timing 
circuit using a square wave generator. Compressed air was supplied at 30 psig 
through tygon tubing into each end of the slot, resulting in a peak injection 
velocity of 8.5 m/sec. The blowing momentum coefficient, C^, was 0.07 for a 

chord Reynolds number of 24,000. 
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Transfer Function Calculations 


The method for calculating nonlinear transfer functions to second order is 
summarized in [RITZ] and [HAJJ]. The general idea behind such work is that a 
nonlinear response of a dynamical system can be condensed, up to second 
order, into additive responses consisting of a linear frequency response. 
L(fm). plus a quadratic frequency response. Q, as given by 



Basically, it involves solving two moment equations in the frequency domain, 
given by [R^TZ] as 

E(Y(f™)X-(fj] = UfJE(X(fJX-(f„)] Z Ql3,E(X(fJX(f3)X-(fJ] 

and 

E(Y(UX-(fJX‘(fj] = MfJE[X(fJX-(fJX'(fO] + Z Q!3jE(X(fJX(fj)X-(yX'(fJ] 

where 

fm = fi±fj = fk±fl 

where ^m is the output frequency of the model, E is the expectation operator 
and X(fj) and Y(fj) are the Fourier transforms of the input and the output 
respectively. Importantly, these two equations allow one to express both 
transfer functions in terms of the various HOS moments that are given by the 
upstream and downstream sensors of the probe (HAJJ). For example, the LHS of 
the first equation is the cross-spectrum, the RHS of the first equation 
associated with the linear transfer function is the power spectrum, and with 
the quadratic transfer function, the auto-bispectrum. For the second equation, 
the LHS is the cross-bispectrum, the term associated with the linear transfer 
function is the auto-bispectrum, and the term associated with the quadratic 
transfer function is the fourth-order spectral moment. Both these equations 
generate a square matrix with enough equations to solve for all the unknowns 
in both the linear and quadratic transfer functions. 

This summer, a large program was written that would assemble these two 
equations into a matrix and solve for both the linear and quadratic transfer 
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functions. The name of this program is SKETCHY. SKETCHY initially reads in 
the input and output time series recorded from the hot-film anemometry 
probe. Following this, the time series is broken up into a number of blocks, 
and the various averaged spectral moments are calculated, such as the auto- 
and cross-spectrum, the auto- and cross bispectrum and the fourth-order 
spectral moment. These moments are then used to construct a large two- 
dimensional square matrix associated with the unknowns, along with a column 
vector for the various terms in the transfer functions. This is solved using a 
lower/upper triangular solver. The unknowns are then reassembled into 
appropriate matrices representing the linear and the quadratic transfer 
functions. The magnitudes of these complex arrays are taken and the results 
are written to output files, which are displayed using a plotting application 
developed by W. Davis using the graphics language IDL (Interactive Display 
Language) by Research Systems Incorporated. To further understanding of 
the quadratic transfer function results, a non-dimensional quadratic transfer 
function is also calculated, defined as (HAJJj 



This nondimensionalized quadratic transfer function shows the magnitude of 
the energy transferred due to a given interaction, combining affinity for a 
given three-wave interaction to exchange energy multiplied by the actual 
energy present to be transferred. 

The program was verified for correctness of results by experimenting with 
sample cases where the linear and quadratic transfer functions were known. 

In particular, two input-output combinations, one to test for a pure linear 
response, and the other to test for a pure quadratic response were analyzed, 
and both cases yielded the correct answer. However, preliminary examination 
of the matrix used to solve for the linear and quadratic transfer functions 
indicate that the matrix is poorly conditioned. One of the highest priorities 
upon continuance of this work will be to use a sparse matrix solver designed 
for poorly conditioned matrices, as well as implement the entire code in double 
precision. Because of the size of the matrices manipulated, this will involve 
use of a supercomputer. 
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Results and Discussion 

Quadratic transfer functions were calculated using 16384 points collected by 
the Masscomp computer, subdivided into 64 point blocks for the averaging 
techniques. The sampling rate was chosen to be 80Hz, well above the 
necessary Nyquist frequency needed to examine the higher frequency 
content. The lowest frequency available for examination was 1.25Hz. This was 
deemed adequate for a preliminary study of the phenomena; however, the 
final study to be conducted this winter will use more data and a 128 point block 
size. 

For the static case, hot film pair measurements were made in a variety of 
positions downstream of the airfoil. Two angles of attack were chosen; 15° and 
20° . Six cases were analyzed for the 15° angle of attack and seven cases for 
the 20° angle of attack downstream of the airfoil. 

Figure 1 is a representative plot of the nondimensionalized quadratic transfer 
function at 15° angle of attack for the unforced case. The transfer function 
shows that the flow is self-organizing, with energy from all frequencies in 
the input being coupled to, or transferring energy to 5 Hz, the lowest 
subharmonic present in the flow, 20 Hz, the main harmonic, and 30 Hz, the 
first superharmonic. The primary frequency triads present in this plot are 
(15Hz, 5Hz, 20Hz), (lOHz, -5Hz, 5Hz), and (20Hz, lOHz, 30Hz). 

Figure 2 is the plot of the nondimensionalized quadratic transfer function at 
20° angle of attack. The flow is also self-organizing, transferring energy 
primarily to approximately 4 Hz, the lowest subharmonic, and 30 Hz, the first 
superharmonic in the flow, with large energy coupling mechanisms at (lOHz,- 
6Hz,4Hz) and (20Hz,10H2,30Hz). 

Because of the large amount of coupling observed between for the frequency 
values of 5Hz, lOHz and 20Hz in the natural flow case, indicating strong flow 
organization because of these frequencies, a pulsed air control at these 
frequencies was used to control flow separation at 20° angle of attack. For 5Hz 
pulsed air frequency, the flow quickly becomes organized and periodic, with 
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strong frequency content at lOHz, 15Hz, and 2SHz. Figure 3 shows the 
nondimensionalized quadratic transfer function for a typical trial. Strong 
coupling exists between SHz in the input and lOHz in the output, from 
frequency triad (SHz, SHz. lOHz). In addition, the pulsed air serves as a strong 
organizer of the flow, as all energy from all frequencies is directed into the 
lOHz. ISHz and 2S Hz frequencies in the output. 

Conclusions 

Preliminary work concerning flow control for separation from a static airfoil 
mounted at given angles of attack was completed this summer by the P.I. The 
Volterra functional series/ quadratic transfer function calculation method 
was implemented and used to obtain results for flow control. Initial results are 
promising; fundamental flow frequencies are successfully identified, and 
energy transfer mechanisms in the flow are highlighted. Future work will 
involve examining the dynamic pitching airfoil for energy transfer 
mechanisms so that a control scheme can be proposed for delaying the onset of 
dynamic stall, and enhancing airfoil performance in the post-stall regime. 
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Figure 1. Nondimensionalized quadratic transfer 
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Abstract 

The research performed pertains to the development and testing of synergistic active and 
passive control designs for the piezoelectric laminate beam experiment at FJSRL. Under the 
.\FOSR 1992 Summer Research Program, a multi-input, multi-output transfer matrix model for 
a slewing beam system with piezoelectric actuators and sensors, was experimentally verified. The 
piezoelectric beam system at FJSRL was used for this work. 

The transfer functions previously established for piezoelectric laminate beam experiment were 
used to develop and experimentally validate a simultaneously optimal active and passive damping 
design for the experimental system. However only the optimal active damping was verified here. 
A preliminary damping design has already been developed and tested at FJSRL. In this study. 
IMSC, (Independent Modal Space Control) was used for an optimal active control. An optimal 
active and passive design synergism has previously been developed. The work reported here 
supports ongoing the progress of contract #RDL 93-156, toward verification of both the active 
and passive control measures together. 
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ACTIVE CONTROL DESIGN FOR THE FJSRL SLEWING 
PIEZOELECTRIC LAMINATE BEAM 


Travis DuBois 


1 Introduction 

It is well known that the achievable performance of active control designs can be limited by a number of 
practical factors, including spillover effects [1] which could lead to instability. Indeed, in initial experiments 
with the JPL precision truss, Fanson, et.al. [2] found that simple controllers successfully attenuated lower 
frequency modes but inadvertently destabilized higher modes. The augmentation of active control with pas¬ 
sive damping provides a practical solution [3] to such problems. Viscoelastic passive damping treatments can 
be designed to perform well in a prescribed frequency range and are particularly well suited for the control 
i.f higlier modes of structural vibration. .Moreover, active control of higher modes is seldom required for 
purposes of, for example, shape control. Rather, enhanced damping of the higher modes is the desired prop¬ 
erty. .Active control measures, are normally viewed as being well suited only for lower frequency structural 
modes, due to actuator bandwidth and computational speed limitations. Considered together, these facts 
suggest that a synergistic approach employing active control for lower frequency modes and passive control 
for higher frequency modes deserves serious consideration. Passive damping can enhance the performance 
of active controllers while easing the burden of active control and providing enhanced system stability [3]. 
Desirable properties of passive control approaches include simplicity, reliability, low cost and light weight. 
The goal of this research was to develop and experimentally validate a simultaneously optimal active and 
passive damping design for the experimental system. The control designs for the beam experiment will be 
based on the models developed and verified during the summer of 1992 at FJSRL. The initial approach to 
optimal active/passive control synergism will employ the Independent Modal Space Control (IMSC) [4], [5], 
[6] method for active control. IMSC was chosen due to its simplicity and the intuitive feel that it lends to 
the active control design problem. 

2 Piezoelectric Laminate Electromechanical Relationships 

A typical section of the beam/piezoelectric layer laminate has the top piezoelectric layer, the base structure 
and the lower piezoelectric layer with the subscripts a, 6, and c respectively. 

In the present development, the a layer will serve as actuator and the c layer as sensor. A voltage 
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.ipplied across the actuating layer a induces a longitudinal stress <t; given by: 


aj(x.i) = E,~^Va(x,t) ( 1 ) 

where dai is the electric charge constant of the film (m/v), and Ea is the Young’s modulus of the film (N/m'). 
This stress in turn generates a bending moment Mja about the composite system’s neutral axis, given by 


A//.(x 




ajwydy = -Ead^ibita + tfc)Vo(i-0 = CaVa{x,t) 


where Ca = ^Ead 3 iu(<o + <6) is a geometric constant, the form of which is determined by the geometry of 
the composite beam-actuator-sensor system. The constant changes if, for example, both the a and c layers 
are used as actuators with opposite polarity, then Ca = £a<^ 3 i*^(<a + 4)- It is also worth noting that the 
expression for Ca is significantly simplified by the fact that the laminated beam's neutral axis coincides with 
Its center due to the symmetry of the added piezoelectric elements. 

.Next an expression is obtained for the voltage across the sensing layer. Note that the strain is related 
to the curvature of the-beam in the following way: 


<c(x,0 = ~ + <0^ 


Due to the piezoelectric effect this strain gives rise to a charge distribution per unit area, q(x.t), which is 


given by: 




where ksi is the piezoelectric electromagnetic coupling constant and gsi is the piezoelectric stress constant. 
The total charge developed on the sensing layer is obtained by integrating q{x,t) over the entire surface area 
of the piezo element, 

Q(t)= u'9(x,t)di=-u/ —-(-to - - - (5) 

A, \2 / (731 dx 

The piezoelectric material is similar to an electric capacitor and the voltage across the two layers is given by 
I lie following formula: 

i;«) = ^ , = c, ^ (61 

Cu'fij —ii) ox 

where C, is a constant incorporating various structural and piezoelectric constants in (5) and (6), C is the 
capacitance per unit area, and w(x 2 - ii) is the piezoelectric element’s surface area. In the next section 
a dynamic modelling method for piezoelectric laminate beams is presented using the relationships given in 
this section. 
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3 Dynamic Modeling 


The system under consideration is a generic model representing slewing flexible beams with piezoelectric 
actuators and sensors. The beam is mounted on rotary base, actuated by a torque motor, and rotating in 
the horizontal plane. The piezoelectric patches are attached to both sides of the beam. In this development, 
one side acts as an actuator and the other acts as a sensor. Thus the system inputs consist of voltage Va 
applied to the actuating layer and the torque T applied by the torque motor. Outputs include the voltage 
\'f sensed by the sensing layer, motor hub angle 0 and beam tip position Ytip. 

L sing the usual assumptions for technical beam theory, the differential equation of motion for the com- 
|iosite beam-film system can be expressed as a Bernoulli-Euler beam equation with an additional term due 
to the actuating layer [8]. The equation is given as 

" liere £ and I represent respectively the effective Young’s modulus and area moment of inertia for the 
composite system. .Notice that if the piezoelectric laminates do not cover the entire beam surface, then 
hotli El and pA are functions of x. In the laminated portions of the beam, the product El is given by 
El - Etila + Ekh + EcE, and the product of the effective mass density p with the composite system’s cross 
sectional area A is pA = paAa + PkAi, -i- PcAe- Since the piezoelectric layers are often thin by comparison 
to the base structure, the present analysis is simplified by assuming that El and pA are uniform over the 
length of the beam. The boundary conditions are: 


y(0,t) = 0 

(8) 

5-y(0,t) d^y{0,t) 

ai“8r 

(9) 

8z> + '• dt^dx - “ 

(10) 

ai> 

(11) 


To allow conversion to transfer functions, the differential equation (7) and boundary conditions (8)-(ll) are 
Laplace transformed. Letting = -pAs^/El, (7) becomes 

= ( 12 ) 


with boundary conditions: 


Y(0,s) = 0 

EIY"(0, s) - hs^Y'iO, s) + T(s) = 0 
ElY"(I,s) + Its^Y'(l,s) = 0 


12-5 


EIY"'(l,s)-M,s^Y{l.s) = 0 (16) 

where the primes indicate spatial derivatives. Together, equations (12)-(16) constitute a linear ordinary 
differential equation in x with mixed boundary conditions — two at x = 0 and two at i = /. The differential 
equation (12) is a Laplace transformed version of a classical Bernoulli-Euler beam equation with the second 
spatial derivative of 14(x,s) as a forcing input. Notice that when the spatial distribution of the actuating 
layer is uniform with respect to x, then the voltage term has a non-zero contribution only at the boundaries 
of the actuating layers. The actuating voltage Va is constant in the interval xi < x < xo, but undergoes a 
step change at each of the boundaries of this interval. Thus, the second spatial derivative of the actuating 
voltage function yields, 

v;"(x, s) = [6'{x - X,) - 6'(x - x.)]V,{s) (17) 

where 6'( ), the spatial derivative of the dirac delta, represents the unit dipole function. It can be seen that, 
in effect, the actuating layer applies moments to the beam at its endpoints xi and X 2 . For a system with N 
piezoelectric actuator/sensor pairs a state-space type representation can be formed: 


■y'(x,s) ■ 


y"(x.s) 


y'"(x,s) 


. y""(x.s) . 



0 10 0 
0 0 10 
0 0 0 1 
0 0 0 


K(x,s) 

y'(x,s) 

y''(x,s) 

V""(X,S) 


0 c<,ra(s) 


^y(x-x,)(-l)* + ' 


Defining 


Z(x,s)ify(x,s) y'(x,s) Y"(x,s) y"'(x.5))^ 


U{x,s) ^ - x,){-iy^^ 

•=i 

equation (IS) can be written in compact form with appropriate definitions of .4 and B as follows; 

Z'{x,s) = AZ(x,s)-i- BUix.s) (19) 

Then using the fact [9] that, 

r 6'’(x)^(x)(fx = (-l)>”(0) 

J — OO 

equation (19) can be expressed for ,V = 1, (i.e. only one piezoelectric actuator/sensor pair) as; 

Z(l, s) = e^'ZiO, s) + (20) 

El 

Now combining (20) with the boundary conditions (13)-( 16) yields 8 linear algebraic equations in 8 unknowns, 
which can be solved simultaneously for the elements of Z{0,s) and Z(l,s). Once the end conditions Z(0,s) 
and Z(t,s) are known a general expression for Z{x,s) can be written as follows; 

{ e^*'Z(0, s), 0 < X < xi 

e^'Z(0, s) -h xi < x < xj (21) 

e^'Z(0,s)-t->l[e^<'-"-'>-e^('-*»^]5^i^, X2 < x < I 
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The matrix exponential function corresponding to (18) can be expressed as 




r /"(*) nx) 

/"'(*) /"(r) 

I3*fix) rix) 
.0*nx) ^/'(X) /3^f{x) 


/(^) 1 
fix) 

/"(X) 

r(^) J 


(22) 


where, f{x) = ^[sinh(/?x) -sin(^x)]. 

The solution expressed above by (21) has been derived for a single piezoelectric element, however, because 
of the system's linearity, superposition allows the same analysis to be applied to systems with any number 
of piezoelectric elements. In the analysis which follows, transcendental expressions are developed for the 
transfer functions of interest. To reduce the complexity of the expressions, the tip mass Aft and inertia It 
.-ire set to zero for the remainder of the paper. 


4 Transfer Functions 


The system under consideration lias two discrete inputs, 14 and T, and three outputs, Ti,p, 9 and V’,. Thus, 
the transfer matrix can be expressed: 


■ 0(s) ■ 


G»,v.(s) 

<j».t(s) 

V.{s) 

= 

<JV..v,(s) 

Gv,.t(s) 

. l'(tr(s) . 



Gr_r(s) 



■ K(s) ■ 


. . 


(23) 


Tlie subscripts p and q of each element Gp,^{s) = Np_^{s)/D{s) identify the associated output and input 
variable respectively. Each of the six transfer functions for this system have the common denominator D(s) 
and are given in [11]. 


4.1 Rational Laplace Domain Transfer Functions 


The transfer functions presented above are exact relationships (for the model used) in terms of transcendental 
functions of 3. The desired form of the transfer functions (23) is in terms of ratios of polynomials in the 
Laplace operator s. For this purpose, the numerator and denominator transcendental functions can be 
rationalized as shown by Schmitz [10] using the Maclaurin series expansion, 

n3) = ^^-jP-{l + k,P^ + hl3^ + ...) (24) 


or in equivalent product form 


/(/?) = 


p! dpp f_l 



(25) 


The derivative term in the leading coefficient is obtained by differentiating the function f(l3) until a 

non-zero constant is obtained. The p is the index of that first non-zero derivative and N corresponds to the 
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number of modes to be retained in the rationalized representation. This procedure leads to the following 
rational expressions needed for this work from (23) in terms of s. 


(d - x.f - (I - X,)’) x,r)) n 

^ 1 = 1 


({s/yry + 1 ) 

((5/n.)- + i) 


(26) 


where the total rotary inertia Ij- is defined /t = fk + 


(27) 


5 IMSC 


Cliven a linear transfer function representation of a flexible structure system in the form: 


Yis) = G(s)F(s) 


(28) 


where: 

F(s) is a vector of control input forces. 
C{s) is a matrix of transfer functions. 
T(s) is a vector of outputs. 

The system output can be expressed as: 


>'(«) = '7.(5)+ ^C7(s) = 

1 = 1 


^ nB.s + C,)U(s) \ 
{s + c,r- + dj ) 


+ 6U{s) 


(29) 


from Appendix A (58). where 17 ,( 5 ) is the Laplace domain form of the ith modal coordinate. Equation (29) 
implies that in the time domain the ith modal differential equation of a system with proportional damping 
and multiple inputs can be written as 


'7'i(0 + 2C,Pir7,(t) + P?i7i(0 = + C,u(t) + fi[t) 


(30) 


where is the ith modal damping ratio. /,(t) the ith modal independent control force, u(t) the n x 1 distur¬ 
bance input vector, and 5, = [5,iS,2 ... Bi„Y the i'^ modal participation vector corresponding to the inputs 
derivative and C, being likewise for the inputs. Since each element /, of the modal control vector is designed 
to depend only upon the corresponding i'* modal displacement and velocity, n representing the number 
of actively controlled modes, then the closed loop system remains decoupled and the control is referred to 
as Independent Modal Space Control (IMSC). The IMSC method has several desirable characteristics [4], 
[5], [6] including simplicity in design and implementation, and relatively good robustness properties. It also 
has some undesirable features, namely, that the favorable robustness and closed-loop modal independence 
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properties cited above are lost when the self-adjointness requirement is not satisfied, and that the number 
of actuators required is equal to the number of controlled modes. 

If the modal control law is defined as 

/.(<) = (31) 

I lie )th modal equation of the closed loop system becomes 

^.(O + (2C.Pi + hi)tii(t) + (p? + gi)r]i(t) = Biu(t) + C.u(/) (32) 

riiis equation is much easier to solve if we assume very light damping and hence remove the term. 

Suppose that the inputs consist of n unit impulse disturbances, (32) can be solved using Laplace transfor¬ 
mation [12] as: 

T]i{t) = sin bit (33) 

bi 


c, = ^B._ 


n,Pi + IT 


= y pr(i - ) + 9i- - -f (36) 

For independent modal space control, the performance index can be expressed as a linear combination of 
modal performance indices as 

n 

(37) 

i=l 

where each modal index J, can be minimized independently. In the work presented here, the modal perfor¬ 
mance indices are defined as 




The first two terms in the integral in (38) represent the total mechanical energy corresponding to the »th 
mode of the closed loop system defined in (32), and the third term in the integral is the control effort w’ith 
a penalty factor a. Substituting (31) and (33) into (38) and performing the integration yields 


T o 


[ 2 q ^ 2a2(p ?-p jj)J 4C,pi-F 2A< 

The modal index J, can be minimized with respect to gi and />, by setting the corresponding first partial 
derivatives equal to zero. This procedure renders the following optimal gains: 


gr = 0 

= -2CiP.-F27<2p.?-PQ2 
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Examination of (31), (32), (40) ,and (41) indicates that in the optimal design presented here the control 
effort is used only for increasing the damping of the closed loop system, and no effort is expended to change 
the natural frequencies of the system. The total performance index for the optimally controlled structure, 
can be obtained by substituting (39), (40) and (41) into (37), i.e.. 


“ .=1 ^ 


This partially optimized index will be further minimized using a designed-in viscoelastic passive damping 


treatment. 


5.1 Implementation of IMSC 

For implementation we want to take the developed equations from above and represent them in a State-Space 
formulation. From (29) it is possible to convert the transfer function to a partial fraction expansion form; 

where n is the number of modes and m is the number of inputs and once again F, is a vector of input control 


forces. From here we can deduce that; 




In this particular experiment there are tw'o inputs, so to simplify this development, let j = 1,2 and 

_ {B,is -F Cii)Fi 4- (Bj^s 4- C,z)F 3 , 

(s-Hc,)2 + d? 

using Laplace Transform rules, this can be converted directly to the second order, time domain modally 
decoupled form; 

rji + 2c,rji + (c^ -(- = BnFi -F 5,2^2 + C,i A + C.eA (46) 

where c, comes from (34) and contains the fth modal damping ratio, B and C come from the pfe, and F and 
F being the input force and its time derivative, respectively. From (46), the independent modal forces are: 

/. = Bn A + fl.zA + C.,F, -(- C. 2 F 2 (47) 


.Solving (46) for ;j, and integrating twice: 


roo yOO rOO yOO r<x> 

7 i = -2ci / T],cH + Bn / Fidi + B^ / Fndt - (c‘+ d;) / / r;, 

Jo Jo Jo Jo Jo 

n OO yOO yOO 

Fidtdt + Ci2 / / F^dtdt 

Jo Jo 


dtdt + 


12 - 10 







12 - 11 






where x is an estimate of x, the state variables /I, B, and C are the plant’s -4, B,C matrices, K is a matrix 
of Kalman gains, and u and y are the plants input and output. We already know that the output and input 
vectors are 

y — Cx-\-Du u =—Cx (51) 

wliere G’ contains the optimal gains from (40) and (41) and D is the plant D matrix. Therefore (50) becomes 

i = {A-K(C + DG)-BC)i + Ky (52) 

substituting in the y from (51) and simplifying 

i = (A - KC- BG)i + KCx (53) 

Now representing the observer in State-Space form 


Aob — ~ BpiG)] Boh = 


So thus far we have the plant model, which has an output of modal displacements and an observer whicli 
lias an output of modal displacements and modal velocities, and j),. From here we have to get the control 
force For implementation reasons it will be much easier to assume small damping and drop the B,jFj terms 
III (47) Therefore 

/. =c.,fi+C.2F2 (54) 


Now we have to keep in mind that /, is the modal input to mode i and not the input to any one particular 
actuator. The individual actuator inputs are the Fj's. To make it clearer what has to be done, (54) rewritten 


for two modes is 


Cn C,2 Fi 

Gj) C22 F2 


The matrix containing the C,, terms will be denoted as C. By multiplying (55) by C 


Gw C\2 
G 21 C 22 


W'e also have to remember the feed-forward term from (43), 6. 


6 Experimental Results 

The experiment consists of a 2 m long uniform aluminium beam of rectangular cross section (76.2 mm x 
G.35 mm), hanging vertically from a very low friction binge. The hinge is in a knife edge arrangement 
and is situated in such a way that the root of the flexible beam is coincident with the hinge axis. This 
also allows PZT actuator/sensor pairs to be located with one edge coincident with the hinge axis. The 
hinge arrangement has effectively zero inertia, however, a hub inertia has been added in order to increase 




pole-zero separation. Because the angular displacements of the hinged joint are very small, it was possible 
to get an adequate measurement of hub angle by measuring translational displacement of a point on the 
arm extending from the hinge. A Kaman KD-2300-10CU non-contact displacement measuring system was 
used for this measurement. This sensor has a rated 3dB frequency response range of static to 50 kHz., a 
measurement range of 25.4 mm, and 0 254 mm midrange resolution. The piezoelectric ceramic elements are 
\erniiron PZT5A patches mounted in pairs on opposite faces of the beam. A 63.5 mm long by 38.1 mm 
wide piezoelectric material patch is attached to each side of the root of the beam so X\ = 0, 12 = 63.5 mm, 
and u’ = 38.1 mm. A second pair is located roughly midway down the beam at ii = 1.048 m, 12 = 1.111 
m. The two pairs of patches were tied together to effectively make one actuator at the root and one in the 
middle of the beam. 

.V signal generator was used as input to the patches for excitation. When the beam’s motion reached 
what appeared to be steady-state the generator was turned off and the controller, on. The controller was 
implemented on a Sun SP.4RCstation equiped with a Solflower input/output device, where the Kaman sensor 
-ignal was the input and the two control voltages for the patches was the output. A C-language program 
which read the discretized .4. B, C, D matrices from a file was used and operated at 0.6 kHz. The two signals 
from the Solflower were each sent to an amplifier and then to their respective piezoelectic patch-pairs. A 
Tektronix model 2630 Fourier Analyzer was used for data acquisition. The parameters for the experimental 
system are tabulated in Table 1. Table 2 contains the parameters of the piezoelectric material. 


Hub inertia, /h 


Beam length, 1 

2.0 m 

Volumetric mass density, p 


Cross sectional area, A 

g|:kl:'kgll.>ilWai 

Young's .Modulus, E 


Area moment of inertia, / 

IQQEBCUli 


Table 1: Experimental Beam Parameters 


Charge Constant, da 

BQQjjySBI 

Voltage Constant, 931 

-11.4 X 10“^ Vm/N 

Coupling Coefficient, kn 

-0.340 

Capacitance, C 


Surface area, w x (ij — ii) 


Thickness ta 

IIBlSQQNESHi 

Young's .Modulus. Ea 





Table 2: PZT5A Parameters 

On this beam there are two actuators so two modes can be simultaneously controlled using IMSC. Due 
to a lack of voltage amplification the controller was implemented for one mode at a time. The amplifiers 
used had a mziximum power output of ±35 volts. It isn’t until about ±25 volts that the these particular 
patches seem to respond with sufficient control authority. As it turned out for a particular a, the voltage 
for patch-pair one was an order of magnitude higher than that for patch two. So for example with the beam 
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Figure 2: Mode 1, Free Time Response, (experimental) 

\ il>rntiiig ill the first mode, and q adjusted so the voltage for patch-pair one was just saturating its amplifier 
at ±3-3 volts, the voltage to patch-pair two would be near ±3.5 volts. Therefore the second patch-pair would 
not he activated. This is acceptable for controlling the second mode. Vet the first mode is best controlled 
with the middle patch. So it was decided to implement the controller for each mode separately. Since mode 
one is best controlled from the second patch-pair, a was adjusted so that when the controller was first turned 
on, the amplifier for patch-pair two just barely saturated. This of course makes for a control voltage for 
the first patch-pair that is far to high for the available amplifier. Therefore the connection for the first 
patch-pair was disconnected. In this configuration, significant damping was achieved as seen in figure (4). 
I he simulation for controlling mode one is shown in figure (5) and correlates well with (4). As previously 
.-i'ated mode two is best controlled using the first patch-pair. Therefore mode two was controlled using the 
same a with the patch-pair at the root connected. The particular a used to control mode one made the 
amplifier that supplied the first patch-pair, just saturate. This is exactly what is needed to get the most 
control authority from the patches. The free time response plot for the beam with mode two controlled 
is shown in figure (8) and the simulation is shown in figure (9). Comparison of the two show a very nice 
correlation between the modeled and actual implementation. Future implementations will include optimal 
passive damping in conjuction with this optimal active control. 

7 Conclusion 

An optimal method for actively controlling a flexible slewing piezelectric laminate beam has been presented. 
The method uses the partial fraction expansion of exact open loop transfer functions of the beam and 
is implemented using Independant Modal Space Control (IMSC). Implementation of the controller was 
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Figure 3; Mode 1, Free Time Response, (simulated) 
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Figure -I: Mode 1. Controller On, Free Time Response, (experimental) 
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Figure 5; Mode 1, Controller On, Free Time Response, (simulated) 
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Figure 9: Mode 2, Controller On, Free Time Response, (simulated) 


conducted at the US Air Force Academy, Frank J. Seiler Research Laboratory. The results show significant 
vibration reduction for the very flexible beam presented here. This work leads directly to the implementation 
of an optimal active and passive control design which uses a constrained viscoelastic damping layer. 


8 Acknowledgement 

This work has been supported, in part, through the Air Force Office of Scientific Research Summer Research 
Program. I wish to thank the staff at the US Air Force Academy, F.J. Seiler Research Laboratory for their 
assistance with the experimental research. In particular, I wish to acknowledge my laboratory focal point 
C’apt. Daniel J. Stech, as well as Dr. Ralph Quan and Dr. Tae Lim. I also want to thank my advisor Dr. 
Thomas E. Alberts for the challenging and rewarding experiences he has offered me. 


A Partial Fraction Expansion 

The partial fraction expansion of a transfer function with a product of m second order terms s~ + 2(^,2iS + 2 ? 
as the numerator, and n second order terms of the form s* + 2C»PiS + p? as its denominator can be expressed 
as; 


n{®+ 

G{s) = i|i- 

n (*++dj 
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therefore 


]^(s^ +2C.r.s + z.*) 


i = l 


ji[(s* + 2Cip,s + p') 

i = » 


ai=0-i l>i = 2i \/l - C? 

Ci = C.Pi = p<\/l - C,^ 

Next we can express the function in the general partial fraction expansion form: 

C(s) = ^( -+= 

+ s + c,-jdiJ ^ \(s + Ci)-+ dr J 


where 


‘={i 

For any r'^ term the first order residue Ar is: 


= m 
m 


(57) 


Q(s + a,)- +6f) 


Ar = {s + Cr+jdr) 


1 = 1 


jj(s + ei)-+cff) 


i = l 


1 


SS~Cr-jdr 

JJ(s + a,)^ + bj 




i-2jdr)’-^ " 

J][(s +C,)-+ <f7 {s + Ci)^ + d- 

’ = * '='•+> \, = -Cr- 3 dr 

m 

]^((-Cr + a,) - jdr)'^ + br 

_l _^=J_ 

n((-Cr + Ci) - jdr)^ + d] n {{-Cr + C. ) - jdrf + d? 

i = l i=r+l 


J^(a, - Cr)* + i? - d* - 2jdriai - Cr) 


isl 


(2dr)^^ 

JJ(c,- - Cr)* + d* - d* - 2jdr{Ci - Cr) (c, - Cr)* + d? - d* - 2jdr(c, - Cr ) 

»=I isr+l 


(58) 


It is not necessary to get ,4" explicitly. From (58): 

Sr S A" Cr Ar 


Al 


(S + Cr )* + d? 5 + Cr + jdr S + Cr - jdr 
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_ /l,(a + Cr- jdr) + v4;(s -H Cr 4- jdr) 

{s + Cry + d^ 

_ (s + c,)(^r + Al) + jdr(A‘ - ^,) 

(s + Cry + dj 

_ 2(5 + Cr)Re{Ar) + 2drIm(Ar) 

(s + Cry+d^ 

Note that the Re() term means real part and Im() means imaginary part. From here we see; 

Br = 2Re{Ar) 

Cr = 2iCrRe{Ar) + drIm{Ar)) 

These two formulas seem a little clumsey with the Re and Im comments, but are very straight forward to 
implement in a computer code. 
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The Rational Design of Novel Semi-Organic Salts for Use as Nonlinear Optical Materials 
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Abstract 


The primary focus of this project was the rational design of nonlinear optical materials 
based on the concept of crystal engineering. It was proposed to use hydrogen bonding and 7t- 
stacking effects to design NLO optimized organic/inorganic salts. The salt series was based 
around the imidazolium cation with anions varying in size and hydrogen bonding ability. The salts 
were screened for second harmonic generation using the powder technique developed by Kurtz 
and Periy. In addition, the salts were characterized by IR, NMR, Raman, melting point, and mass 
spectroscopy data. 
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Introduction to Nonlinear Optics 


Joan Fuller 

Nonlinear optics is the study of the interaction of intense laser light with matter and how 
this interaction modifies the optical properties of the material * Nonlinear optical materials are 
divided into two classes; second-order and third-order materials. Second-order materials have 
changes associated with changes in the first nonlinear term of the polarization equation. They are 
characterized by the molecular hyperpolarizability term, p. Third-order NLO effects are 
characterized by the second hyperpolarizability, y. At this time, all commercial applications use 
second-order nonlinear optical effects; no third-order devices are being produced. In order to 
design materials that maximize the second-order NLO response it is necessary to understand the 
molecular and macroscopic polarization of the material. Equation 1 expresses the molecular 
polarization of a material in an applied electric field.2 

P = p + a E + PEE + y E E E + .. 

Equation I 

p = permanent dipole moment 
a = linear polarizability 
P & y = first and second hyperpolarizabilities 
E = applied electric field 
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The individual molecular responses interact to give the bulk response, also called the 
macroscopic polarization of the material, expressed by Equation 2.2 

P = Pq ^^ E + x^ E E + x^ E E E + 

Equation 2 

Po = spontaneous polarization of the material 

X ^ = linear susceptibility 

X 2 & X second and third order susceptibilities 

NLO Effects 

Several nonlinear optical effects are known but perhaps the best known is second- 
harmonic generation (SHG). Originally observed in single crystal quartz, SHG was the first NLO 
interaction to be described.^ In an SHG measurement a laser beam of frequency co is focused on a 
sample, and coherent light at twice the initial frequency 2 (d is detected. Because SHG 
measurements can be made on single crystals, poled polymers, films, or powders, it is a common 
screening technique used to determine NLO activity. The molecular hyperpolarizability term 
responsible for SHG effects is (3. In general, the higher the |3 value, the larger the NLO effect of 
the material .2 

Two other important second-order nonlinear effects are sum-frequency generation (SFG) 
and difference-frequency generation (DFG). In SFG, two input waves, © i and © 2 , are added to 
generate a new wave of frequency ©3 (©1 + ©2 = © 3 ). In DFG, two input waves, ©1 and © 2 , are 
subtracted to generate a new amplified wave of frequency ©3 (©1 - ©2 = © 3 ); therefore, DFG is 
also referred to as optical parametric amplification.** 
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Second Order NLO Materials 


Compounds that exhibit nonlinear optical properties are typically divided into three 
categories: 


* inorganic * organic * semi-organic 

Inorganic 

Inorganic compounds are the only commercially available NLO compounds, in part 
because they were the first class of materials to be studied, but also because they have the desired 
physical properties needed for useful materials. However, there are limitations to this class of 
materials. They have small y}- values, and there are few opportunities to alter the properties 
beyond what can be achieved through a restricted number of substitutions. Furthermore, many of 
the inorganic crystals must be grown from high temperature melts which leads to growth 
impurities and the introduction of crystal strain due to the cooling process. Problems such as 
these, result in sample to sample variations that make it difficult to design operative systems. ^ 
Table 1 lists examples of inorganic NLO systems that are currently in use and their relative merit 
as materials.^ 


Table 1: Inorganic NLO Systems ^ 


Compound 

FOMa 

b - BaB 204 

26 

KH 2 PO 4 (KDP) 

1 

KTi0P04 (KTP) 

215 

Urea 

6 


a. FOM: Figure of merit, a measure of the operative efficacy relative to KDP of the NLO material in a systems 
application. 
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Organic 

Organic compounds offer many advantages over inorganic compounds. The design 
flexibility of organic chemistry provides the opportunity to modify the response and spectral 
properties of the material In addition, there are examples of organic systems with frequency 
conversion responses orders of magnitude larger than the best inorganic systems. However, 
organic systems do have their problems. First, it is difficult to grow high quality crystals. 
Second, the interactions between molecules and between the crystal planes is usually relatively 
weak, which causes the majority of these materials to lack the physical properties needed for NLO 
materials ^ Finally, organic compounds crystallize in a centrosymmetric space group 75% of the 
time and, therefore, are useless as NLO materials unless poling experiments are performed. 

The largest single class of second-order organics consists of donor-acceptor substituted 
aromatics The donors can be amino or dimethyl amino groups, oxazole, dithiafulvene, and 
dihydropyrazole. The acceptors can be nitro, pyridinium ions, llilvenes, dicyano or tricyanovinyl 
groups, sulfones, enones, or imines. Aromatics are used to establish highly polarized charge 
distributions and can include, benzenes, stilbenes, diarylacetylenes or biaryls.‘^ Table 2 lists some 
promising donor-acceptor systems along with their molecular hyperpolarizability values. It is 
obvious that these donor-acceptor materials should have high bulk SHG values if they can be 
made to crystallize in acentric space groups. 
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Table 2: Donor-Acceptor Aromatic Systems 


Donor 

Aromatic 


3 X 10*^® esu 

NH2 

Benzene 

CN 

3.1 

NMe 2 

Benzene 

CN 

5.0 

NMe 2 

Benzene 

NO 

12 

NMe 2 

Benzene 

N02 

12 

NMe 2 

Stilbene 

N02 

73 

NMe2 

Stilbene 

CN 

36 

OMe 

Stilbene 

NO 2 

34 

NH2 

Stilbene 

NO 2 

40 


* Urea has a (3 value of -^.5 x lO'^^esu. 


Semi-organic 

Semi-organics, or mixed salts, are the newest class of materials studied for NLO 
applications They combine the best of the organic and inorganic systems. Perhaps their greatest 
advantage, is their straightforward synthesis that allows them to be rapidly screened for SHG 
activity. Two general categories of semi-organics have been defined: Type I: Salts containing an 
NLO active organic ion and an inorganic counterion, and Type II; Coordination complex salts 
containing a neutral NLO organic ligand coordinated to a metal ion. Tables 3 & 4 lists examples 
of Type I and Type II materials. 
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Table 3 ; Type I Semi-organic Systems 


Cation 

Anion 

SHG Value vs Urea 

4-CH30C6H4CH=CHC5H4NCH3 

CF 3 SO 3 

50 

4-CH30C6H4CH=CHC5H4NCH3 

BF 4 

0 

4-CH30C6H4CH=CHC5H4NCH3 

CH 3 C 6 H 4 SO 3 

120 

4-CH30C6H4CH=CHC5H4NCH3 

Cl 

60 


Table 4: Type II Semi-organic Systems 


Compound 

SHG Value vs KDP 

Cd(thiourea) 2 Cl 2 

2.75 

Zn(thiourea) 3 S 04 

3.6 

Pb(diethyldithiocarbamate )2 

12 


In these systems the organic molecule provides the source of nonlinear polarization, and offers the 
flexibility needed to design working NLO systems. The inorganic component provides the strong 
three dimensional binding network that gives the mechanical stability. Another advantage of these 
systems is the ability to crystallize them at ambient temperature, thus, creating less strain in the 
crystal and making commercial applications less expensive. 

Methodology 

A series of Type I semi-organic NLO materials listed in Table 5 were prepared and 
characterized. The salt was composed of the 1,3-disubstituted imidazolium cation, illustrated in 
Figure 1, and specific anionic groups that have been shown to alter the NLO properties of 
materials. These salts often stack in flat, one-dimensional layers which are promoted by hydrogen 
bond networks. Because of the inherent lack of symmetry that can be built into the cation, it is 
possible to promote crystallization into acentric space groups. The imidazolium salts prepared in 
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this study are soluble, stable, cheap, and are easily grown into large single crystals. The research 
was conducted in two phases. 



Rj=H, Methyl, Benzyl 
R 2 =H, Methyl, Ethyl 

R3=H 


Figure 1 

Phase I: Synthesis of stable imidazolium salts as prototypes for SHG active n^aterials 
Listed in Table 5 are the salts which have been prepared. The synthesis of 1- 
benzylimidazolium tetrafluoroborate is given as an example of the method used to prepare the 
salts. In 50ml of water 2.005 g (0.0127 moles) of l-benzylimidazole was dissolved with slight 
heating. While stirring this solution 2.891 g (0.0158 moles) of 48% HBF 4 was added dropwise. 
The solution was then rotoevaporated down to dryness. The salt was dried under vacuum at 60°- 
90® C for 24 hours. The reaction occurred in high yield and the product was confirmed with mass 
spectroscopy. Recrystallization was attempted from appropriate solvents. 


Phase II. SHG screening and structural characterization 
The apparatus for measuring the SHG activity of powders has been previously described 
by Kurtz and Perry. Preliminary powder SHG measurements were determined by using 1064- 
nm (SH frequency at 532-nm) fundamental radiation. All measurements were referenced to urea 
which was measured at the beginning of each experiment. Because of problems with particle size, 
absorption, and fluorescense it was not possible to obtain quantitative information on the value of 
these salts as nonlinear optical materials. However, by noting the presence of green light (532- 
nm) at the sample, it was possible to determine a qualitative SHG result. Table 6 lists the results 
from these initial qualitative measurements. It is hoped in the future that more quantitative 
measurements can be made. This is not presently possible with the equipment that is available at 
the host laboratory. 
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Table 5; Meltine Point (OC) and Physical State 


S 


1-Methyl- 

imidazolium 

69-73 
white solid 

60-62 
white solid 

1 -Benzyl- 
imidazoiium 

94-95 
white solid 

clear oil 

2-Methyl- 

imidazolium 

158-159 
white solid 

176-178 
white solid 

2-Ethyl- 

imidazolium 

93-95 
white solid 

clear oil 



hydroscopic clear oil 
white solid 


white solid 


60-62 120-122 
hydroscopic white solid 
white solid 


122-123 

clear oil white solid 


82-85 
white solid 




Table 6; Second Harmonic Generation Results 


S 


1 -Methyl- 
imidazolium 

NO 

1-Benzyl- 

imidazolium 

NO 

2-MethyI- 

imidazolium 

NO 

2-Ethyl- 

imidazolium 

NO 



* sample fluoresces at 514-nm 
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Melting point data was determined from a Perkin Elmer differential scanning calorimeter. 
In the cases where a visual observation was necessary to check for decomposition a standard 
capillary melting point apparatus was used Infrared data was obtained on a Nicolet 510P FT-IR 
spectrometer. Spectra for the crystalline salts were made by using KBr pellets. Salts which were 
oils were measured neat. Raman spectra was recorded at the University of Alabama on a Yobin 
UlOOO Laser Raman spectrometer. Measurements were made using the 514-nm wavelength. It 
was believed that some of the salts were undergoing fluorescence as they generated the SHG 
frequency, 532-nm. If this was the case then it is possible that they would also fluoresce when 
exposed to the 514-nm wavelength of the Raman measurement. 

Mass spectroscopy was used to evaluate the purity of the salts. These measurements were 
performed on a Hewlett Packard 5985B mass spectrometer by Mr. Jack Greathouse. The salts 
were analyzed via a direct insertion probe over a temperature range of 250-250° C. Nuclear 
magnetic resonance data was obtained on a Varian Gemini 300 spectrometer. Both and 
were recorded for the salts. The solvent used was dimethylsulfoxide with tetramethylsilane as the 
internal reference. 
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Table 7a: IR, Raman, and NMR Data 


S 


1-Methyl- IR:3150, 
imidazolium 3500, 1790, 
1400, 1040, 
819 cm' ^ 


Raman; 
1C39, 1322, 
716 cm" * 


! NMR: JhS;; 
8.9(s), 

7.6(d), 4. l(s) 


13c 5,; 136, 
123, 119.5, 
35 


IR: 3150, 
3450, 1590, 
1554, 1015, 
400 cm' 1 


Raman: 1105, 
611,983, 
450 cm'i 


NMR: 1H5; 

8.35(s), 
7.4(s), 7.3(s), 
3.8(t) 

13c 5; 13 7, 
124.2, 122.2, 
34.4 


IR; 3450, 
3120, 1120, 
1063, 
400 cm'l 


Raman; 1070, 
777, 533, 
360 cm'l 


NMR: 1H5; 
9.0(s). 7.7(d), 
3.8(m) 

13c 6; 136.2, 
123.5, 120.1, 
35.7 


Br 


IR: not 
available 


Raman; not 
available 


NMR: 1h 6; 
9.3(s), 7.8(5), 
7.7(s), 4.0(5) 


IR: 3157, 
1253, 1032, 
837, 758, 640, 
518 cm'l 


Raman; not 
available 


NMR; 1H8; 

9.03(5), 
7.7(m), 3.87(5) 


13c 5; 135.8, 13C5;136.2, 

123.5,119.7, 123.5,120.1, 

36.63 35.8 


Table 7b; IR, Raman, and NMR Data 


NOi' S 


1-Benzyl- IR: 3200, IR;3170, 

imidazolium 1049,1039, 1190,810, 

970, 820 cm' 420 cm'l 

1 


Raman: Raman; 1105, 



1354, 1209, 
1039, 
826 cm' 1 


NMR; 1H5; 


983, 611, 
450 cm'l 


NMR: 1H5: 


9.3(5), 7.8(5), 8.7(5), 7.5(5), 

7.7(5), 7.4(5) 7.4(5), 5.4(5) 

5.4(5) 


IR; 3410, 
3100, 1820, 
1060, 1015, 
900 cm'l 


Raman; 1070, 
777, 533, 
360 cm'l 


NMR; 1H6; 
9.2(5), 7.8(5), 
7.7(5), 7.4(5), 
5.4(5) 


Br 


IR: 2850, 
2600, 910, 
650 cm'l 


Raman: 1434, 
997, 803, 

608, 327 cm' 
1 


NMR; 1H6; 
9.4(5), 7.8(5), 
7.5(5), 7.4(m) 
5.5(5) 


13c 5; 135, 13c 6; 136.6, I3c 6; 129.0, 13c 5; 135.5, 

128.5, 122, 128.6, 124, 120.8, 52.0 129.0, 122.4, 

120,51.5 121.4,51.1 120.6,51.9 


IR: not 
available 


Raman; not 
available 


NMR: 1h5; 
9.3(5), 7.8(5), 
7.7(5), 7.4(5), 
5.4(5) 

13c 6; 135.8, 
129.4, 120.7, 
52.1 









































Table 7c: IR, Raman, and NM 

R Data 


NOr 

SO 4-2 

BF 4 - 

Br 

■ssssm 

2-Methyl- 

imidazolium 

IR: 3111, 
1576, 1388, 
1078, 

711 cm"^ 

IR: 3190, 
3005, 2650, 
1185, 850, 
610 cm‘l 

IR: 3450, 
3060, 1700, 

1050, 865 cm" 

1 

IR: 3420, 
2950, 2670, 
1300, 

840 cm‘1 

IR: not 
available 


Raman: 
1510, 1110, 
1042, 670, 
150 cm"^ 

Raman: 1105, 
983,611, 

450 cm*l 

Raman: 762, 
730, 668 cm-1 

Raman: not 
available 

Raman; not 
available 

NMR: 1H5; 
7.5(s) 

13c 6; 144.5, 
118.9, 114 

NMR: 1H5; 
9.3(s). 7.3(s), 
2.5(s) 

13c 5; 119.9, 
12.4 

NMR: 1h 5; 
7.5(s), 2.5(m) 

13c 5, 144 5, 
118.9 

NMR: 1H5; 
7.5(s), 3.4(s) 

13c 5; 144.0, 
118.9, 

NMR: 1H5; 
7.5(s), 3.9(s), 
2.6(s) 

13c 5; 144.5, 
118.9, 11.4 


Table 7d; IR, Raman, and NM 

■1 Data 


NOv 

SO 4-2 

BF 4 - 

Br- 

WESSSSSM 

2-Ethyl- 

imidazolium 

IR; 3400, 
3150, 1640, 
1410, 1050, 
930cm" 1 

IR; 2810, 
2000, 1680, 
1040, 985, 

510 cm"l 

IR:not 

available 

IR: 3010, 
2740, 

810 cm-1 

IR: not 
available 


Raman; 
1511, 1116, 
1043, 643, 
136 cm"l 

Raman; 1105, 
983,611, 

450 cm’l 

Raman: 1070, 
777, 533, 

360 cm'l 

Raman: 1510, 
1123,957, 
917 cm'l 

Raman: 
1516, 1127, 
1034 cm'l 

NMR: 1h 5, 
7.6(s), 2.9(m), 
1.3(m), 0.9(s) 

13c 5; 148.9, 
118.9, 19.2, 
11.6 

NMR; 1H6; 
9.4(s), 7.2(s), 
2.8(m), 1.3(m) 

13c 6; 149.1, 
119.9, 20.1, 
12.3 

NMR: 1h 6; 
7.5(s), 2.9(m), 
1.3(m) 

13c 5; 149.1, 
118.8, 19.1, 
11.4 

NMR: 1H6; 
3.5(s), 2.9(m), 
2.8(s), 1.3(s), 

13c 5; 116.4, 
37.6, 9.14 

NMR: 1H5; 
10.6(s), 
7.5(s), 
2.9(m), 
1.3(m) 

13c 5; 149.1, 
119.2, 19.4, 
11.6 
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Conclusion 


Materials research and development for nonlinear optical applications has progressed to 
the point where several systems are commercially available. So far, no material system has proven 
to possess all the desired properties for nonlinear optics. Thus, considerable research is still 
needed to develop materials that can fulfill the critical requirements of devices used in information 
processing, optical frequency conversion, integrated optics, and telecommunications. 

This project was limited to a specific salt series based around the 1, 3-disubstituted 
imidazolium cation. However, given the fact that there are virtually limitless numbers of 
permutations of cations and anions this project could lead to a massive growth in the number of 
known SHG active materials. Based on the positive results achieved in this preliminary study, 
work is continuing on extending the salt series to improve the NLO activity and the desired 
structural properties. In addition, improvements on the currently available SHG apparatus are 
planned which would allow quantitative measurements to be recorded. A separate project using 
X-ray crystallography to evaluate the structure of these salts is also underway. 
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A bstract 

Structural damage detection in a twenty-bay planar truss was accomplished using an 
artificial neural network. Instead of using natural frequencies and mode shapes, the frequency 
response functions (FRF’s) experimentally obtxiined from accelerometers at two locations on the 
truss were directly used to distinguish among damage cases and to train the network. Unlike 
conventional approaches based on system identification techniques, the neural network approach 
does not require an analytical model of the structure. The direct use of the FRF’s eliminates the 
need for modal parameter identification. Out of the 60 damage cases considered in this study, the 
neural network was able to identify uniquely 21 damage cases and narrowed 38 others down to 
two possible damaged struts. 
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AN ARTinCIAL NEURAL NETWORK APPROACH TO STRUCTURAL DAMAGE 
DETECTION USING FREQUENCY RESPONSE FUNCTIONS 


Clinton R. Povich 
Tae W. Lim 


Introduction 

To maintain the performance and safe operation of aerospace vehicles, structural integrity 
must be monitored periodically. Numerous studies have been conducted to detect structural 
damage using measured mode data and system identification techniques [1-4], Recently, attention 
was given to the application of the artificial neural network to structural damage detection [5,6]. 
The natural frequencies and mode shapes are used to train the network and to conduct damage 
detection. Unlike the conventional approaches based on system identification techniques, the 
neural network approach does not require a refined model of the undamaged structure. The fact 
that damage detection can be performed using the measured mode data alone without an analytical 
model relieves the burden of conducting analytical model refinement to match measured mode data, 
which is typically costly and time consuming. 

This paper investigates an effective use of artificial neural network to conduct structural 
damage detection. In a typical modal test, the modal parameters are identified using the frequency 
response functions (FRF’s). Instead of using mode shapes and frequencies, the neural network 
can be trained directly using the FRF’s. Some of the advantages using the FRF’s directly include; 
(1) the modal parameter identification process is not required and thus the error and cost associated 
with it are eliminated, (2) any nonlinear behavior exhibited in FRF’s can be incorporated in the 
damage detection process, and (3) the number of transducers required to perform damage detection 
can be significantly reduced as will be shown in the twenty-bay truss structure example. 
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Test Apparatus and Procedure 


The focus of this study is damage detection in a twenty-bay two-dimensional truss (see 
Fig. 1). The 7.07 m long truss fixed at one end and free at the other lies flat on a table. The truss 
consists of Meroform hollow aluminum tubes and steel node balls. Underneath each batten, a steel 
bar is bolted to the node balls. Each steel bar rests on two steel ball bearings that allow the truss to 
move in the horizontal plane. Detailed information on the truss is available in Ref. [7]. At the end 
of the truss, a MB Dynamics Modal 50 Exciter (shaker) is attached to provide force excitation. 
Accelerometers are located at the numbered positions as shown in Fig. 1. Since the minimum 
number of accelerometers needed was unknown, eleven accelerometers were installed to provide 
sufficient modal test data. Then, the minimum number of accelerometers and their locations were 
to be determined by inspecting the FRF’s and evaluating the damage detection performance of the 
trained neural network. Consequently, only the FRF’s from numbers 5 and 11 are selected to train 
the neural network and to conduct damaged detection. Accelerometer 11 is chosen because the 
FRF’s at this point displayed the greatest diversity. The other accelerometer was used to provide 
uniqueness to the few damage cases that had similar FRF’s from accelerometer 11. 

A different strut is removed for each of the 60 damage cases investigated to simulate 
complete damage in that strut. In each bay, only the diagonals and longerons are removed. It is 
difficult to remove the battens because of the steel bars connecting the nodes of each batten. Thus, 
they are eliminated from damage consideration in this study. The FRF’s are then found using the 
Tektronix 2630 Fourier Analyzer. The Fourier Analyzer controls the random force input to the 
shaker and records the FRF’s over a range of 0 to 50 Hz, which were obtained using 20 
averagings. The FRFs obtained for all damaged and undamaged cases contain at least four 
bending modes below 50 Hz. Figures 2 and 3 show the FRF’s measured from accelerometers 5 
and 11, respectively. The figures show the changes in FRF’s due to damage in upper longerons, 
diagonals, and lower longerons. Each figure has 21 FRF’s corresponding to the undamaged case 
and the damage cases in bay 1 through 20. Thus, the FRF corresponding to damage case number 
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10 for the upper longeron case in Fig. 2 indicates the FRF measured at accelerometer 5 with upper 
longeron removed in bay 9. Figure 4 shows in detail the FRF of a typical damage case plotted 
with the FRF of the undamaged case at accelerometer 11. 

The sample rate employed to measure the FRFs was 2.56 times the bandwidth (50 Hz) 
and 512 samples were taken for each test. Since only data up to 50 Hz is desired, each FRF 
consists of 200 frames (201 data points). Of these 201 points, the first four are not included in the 
neural network training due to the large ‘spike’ in the data at these points in nearly every FRF of 
every damage case (see Fig. 4). The resulting FRF’s then have a range of 1.25 to 50 Hz (197 data 
points). Finally, the FRF data is convened to a Matlab file to conduct neural network training. 

Artificial Neural Network 

The Matlab Neural Network Toolbox [8] was used to construct and train the neural 
network. The type of network used in this study was a tree-layer neural network using log- 
sigmoid neurons. The backpropagation method was used to train the network. 

A neural network is composed of one or more layers and each contains a series of 
processing units (neurons) that act in parallel within the layer. Data arranged in a vector (or a 
matrix of vectors) is transferred through the layers. Each data value within a vector is multiplied 
by the weight associated with each neuron in a layer and the results is added to the neuron’s bias. 
Then, the results are summed at each neuron creating a new vector with a size equal to the number 
of neurons in that layer. A transfer function converts this new data set and sends it to the next 
layer. The transfer function used in this study is a log-sigmoid which converts any real number to 
a value between 0 and 1. In the target output for the network, a ‘ 1 ’ signified no damage in a strut 
while a ‘0’ indicated complete damage (missing). 

Initially, the weights and biases for all of the neurons are randomly chosen. After one pass 
(epoch) through the network, the output matrix is compared to the target output matrix. The neural 
network computes the delta error vectors (derivatives of error) which are backpropagated through 
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the network until delta vectors are found for each layer. The learning rate is then used to change 
the weights and biases in the direction of steepest decent with respect to the delta vectors. For 
example, the change in a weight, delW, and a bias, delB, are found as 


delW(i,j) = lr*D(i)*P(j) 

(1) 

delB(i) = lr*D(i) 

(2) 


where 


Ir = learning rate 

D(i) = delta vector of a layer 

P(j) = input vector to a layer 

i = number of neurons in layer 

j = number of elements in input vector to a layer 

The inputs are sent through the network layers again and the process is repeated. This process is 
referred to as training. Training stops when the sum of the squared errors for each training vector 
is less than the error goal specified by the user. In this study, the error goal is selected as le-10 to 
confine individual errors within about +le-6. 

An adaptive learning rate allows the network to be trained faster by comparing the new 
error vector after each epoch with the previous (old) error vector. If the new error vector is larger 
than the old (by 1.04), the new weights and biases are discarded and the learning rate decreases by 
0.7. If the new error vector is less than the old, the learning rate is increased by 1.05 and the new 
weights and biases are kept. Also, momentum is used such that the network avoids getting stuck 
in a local minimum of the overall error surface. Basically, momentum allows the network to 
ignore small changes and respond to the general trend in the error surface. A momentum constant 
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is chosen before training with a value between 0 and 1. The quantity 0 signifies a weight change 
based solely on the gradient while 1 sets the new weight change equal to the last weight change and 
the gradient is ignored, i.e., 

delW(i,j) = mc*delW(i,j) + (1 - mc)*lr*D(i)*P(j) (3) 

where the quantity me is a momentum constant A typical value is 0.95. 

For the network in this study, 61 input vectors are used to train the network; 60 damage 
cases and 1 undamaged case. Each input vector contains the FT^F’s with the FRF from 
accelerometer 5 listed first than that of accelerometer 11. Thus, the size of the input matrix 
becomes 394 by 61. Two hidden layers of 125 and 40 neurons, respectively, are included in the 
network to achieve the desired error goal with a reasonable amount of training time. The target 
output vector corresponding to each input vector was 60 element vector consisting of 1 ’s (denoting 
struts with no damage) and a 0 (damaged strut). The undamaged case contains all I’s. Therefore, 
the size of the target matrix is 60 by 61. Each element in the target vector corresponds to a truss 
strut with following pattern 


Element Strut description 


1 

2 

3 


^ Bay 1 


upper longeron 
diagonal 
lower longeron 


5 ^Bay2 


upper longeron 
diagonal 
lower longeron 


58 V 

59 > Bay 20 

60 


upper longeron 
diagonal 
lower longeron 


14-7 




The planar truss structure that is used to perform damage detection exhibits interesting 
dynamic characteristics mainly due to the lacing pattern of the diagonals. Each longeron does not 
produce unique changes in the FRF’s. The two longerons located at the open end of each ‘V’ 
created by two diagonals in adjacent bays produced nearly identical FRF’s when either was 
removed (see Fig. 5). This resulted in 19 pairs of longeron damage cases with indistinguishable 
dynamic characteristics. Therefore, the best that can be achieved for longeron damage cases was 
too narrow down from 60 damage -andidates to two adjacent longerons. Also, when the lower 
longeron in the 20th bay is removed the resulting FRF is similar to that for the undamaged case; 
thus making it difficult to detect damage in that strut (see Fig. 6). To accommodate these 
characteristics of the truss structure, the target output vectors are modified as follows: for the 19 
longeron pairs, each vector contains two O’s representing the two possible damage cases and for 
the 20th lower longeron, no O’s are present. All the diagonals produced dynamic characteristics 
that are uniquely identifiable. The number of epochs taken for the network described above to be 
trained was between 1200 to 1500 when random initial weights and biases were used. Figure 7 
shows the error distribution between the target output and actual output computed using the final 
weights and biases of the trained network. Most of individual errors are confined within ±le-6 
thanks to the small summed squared goal of le-10. 

To evaluate the performance of the trained network, the FRF’s of several damage cases 
which were not used to train the network, were presented to the network. Figure 8 shows the plot 
of the resulting output vectors of several damage cases. To clarify the location and magnitude of 
damaged elements, each element of the output vector was subtracted by 1. Thus, the quantities 0 
and 1 represent no damage and complete damage, respectively, in the figure. The location and the 
magnitude of damage is clearly indicated in the plot In nearly all damage cases, the location and 
magnitude of damage was clearly identified. For longeron damage cases, two struts are identified 
because of the reasons discussed previously. 


14-8 



Conclusions 


Successful damage detection was accomplished using the artificial neural network trained 
with experimentally obtained FRF’s from just two accelerometers on a two-dimensional truss. The 
neural network successfully located a missing strut in 21 cases out of the 60 damage cases 
investigated. The 21 cases include all diagonal damage cases and the lower longeron damage in the 
first bay. In 38 of the 60 damage cases, the network narrowed the 60 candidates down to two 
adjacent longerons. For one case, which corresponds to the lower longeron in the 20th bay, it was 
not successful in locating the missing strut. The success of damage detection was directly tied to 
the uniqueness of the FRF that was produced by a certain damage case. For instance, due to the 
lacing pattern of the diagonals of the truss, 19 longeron pairs produced almost identical FRF’s. 
Thus, pinpointing the damaged member was impossible. 

A neural network trained using natural frequencies and/or mode shapes may have produced 
the same results but the number of accelerometers required would have been at least four to be able 
to distinguish four bending modes below 50 Hz. The direct use of FRF’s investigated in this 
report eliminates the need for the modal parameter identification process and helps avoid the cost 
and error associated with it. Also, a refined analytical model is not required unlike the 
conventional approaches based on system identification techniques. Thus, the burden of 
performing modal refinement to obtain high fidelity analytical model is removed. 

Further study should include the response of the trained network to partially damaged struts 
and/or combinations of damaged struts and the need to include these in training. Additional FRF’s 
from other accelerometers may also be needed for these cases. The FRF of an axially mounted 
accelerometer may help separate the indistinguishable pairs of damage cases. A similar study on a 
three-dimensional truss might prove more successful since the number of damage cases that 
produce almost identical FRF’s may be fewer. 
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Abstract 


Frequency doubling in germanium doped silica planar waveguides deposited on fused silica 
substrates is studied. It is demonstrated that an externally applied, periodic DC field can cause 
instantaneous frequency doubling in these waveguides. The periodicity which causes frequency 
doubling corresponds to the beat length betsveen fundamental and second harmonic light 
propagating in the waveguide. A current has been measured from these periodic electrodes which 
corresponds to the generation of a periodic DC field internal to the glass when illuminated with both 
fundamental and second harmonic light. 
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ELECTRIC FIELD INDUCED SECOND HARMONIC GENERATION IN GERMANIUM DOPED 


SILICA PLANAR WAVEGUIDES 


Peter S. Weitzman 


Introduction 

EfRcicnt frequency doubling of L064nni radiation was first observed in germanium doped 
silica fibers in 1986 (1). Recently this has also been achieved in planar dielectric waveguides made 
from the same materials |2|. Second harmonic generation (SHG) in amorphous materials is 
interesting due to the fact that the second order susceptibiity, is zero. There have been several 
theories proposed to explain this phenomenon (3,4,5,6|. All of these theories are based on Stolen and 
Tom's original hypothesis that an internal DC electric field is responsible for an effective >nd 
that the periodicity of this field provides the quasi phase matching necessary for efficient SHG |7]. 
Assuming that phase matching has been achieved, the conversion efficiency is dependent on the 
overlap integral between the fundamental and second harmonic optical fields and the DC electric 
field, w hich is given by the equation 


x'"\lE^E,^E.E,dA 


( 1 ) 


w here q is the conversion efficiency and the integral is over the transverse cross section of the 
waveguide. This effect has been studied extensively in fibers, however the planar waveguide 
geometry provides easy access to the core region of the waveguide where this internal electric field 
can be modified and probed. 

In 1975, Levine and Bethea (8| applied a periodic electric field to a liquid planar waveguide. In 

1989, Kashap |9| applied an external periodic DC field to the exposed core of an optical fiber and 
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wa$ able to induce SHG. In this paper we present results of using a similar technique on a planar 
dielectric waveguide. We have applied a periodic electric field to the planar waveguide and obsened 
SHG and have also applied both 1.064 micron (IR) and 0.S32 micron (green) radiation 
simultaneously to a planar waveguide and observed the buildup of the internal electric field by 
monitoring the current in the periodic electrode structure. 

In order to provide the correct periodicity for quasi phase matching, the nonlinearity must be 
strong only at points where the fundamental and second harmonic are in phase. This periodicity, 
A/ is given by the expression 


Ar = 


2;r 




( 2 ) 


where |)^ is the propagation constant for the IR light and p2cD propagation constant for the 
green light. The planar w aveguide geometry provides a simple means for determining these 
propagation constants by prism coupling. Figure 1 Illustrates the geometry of the planar waveguide 
and coupling prisms. The propagation constants are determined from the internal prism angles by 
the equation 


P=nk,^ cose. 


(3) 


w here np is the prism refractive index, 6p is the incidence angle of the light internal to the prism 
and kq is the optical wavenumber in free space. These propagation constants can also be 
determined from knowledge of the film thickness and refractive index by solving the appropriate 
eigenvalue equation |10]. The same eigenvalue equation can also be used to determine the film 






thickness and index of refraction front the observed coupling angles for a number of modes in the 
waveguide. 

Methodology 

The waveguides used to obtain the results in this paper were fabricated using argon ion 
beam sputtering of silica and germanium onto commercial, optical grade fused silica substrates. The 
sputtering rate was O.llnm/scc in a partial pressure of 10*^ Torr of oxygen. The two films used in 
this study were approximately 2.0 and 3.S microns thick. Auger analysis indicated that the Ge02 
content of the films was approximately 6m%. The refractive index of the substrate is 1.4607 at 
0.5J2^m and 1.4496 at 1.064pm. The refractive index of the film was measured to be I.S46 at 0.532 
pm and 1.532 at 1.064pm using the prism coupling method. 

Figure 2 illustrates the geometry of the planar waveguide used along with the periodic electrode 
structure mounted on the surface of the waveguide. This waveguidc/electrode combination is 
mounted in the setup illustrated in figure 3. Some of the output from a Q switched, mode locked 
Nd:YAG laser is frequency doubled in a KTP crystal, subsequently the Fundamental and SH beams 
are separated into two paths. Using this technique, the input coupling angle of the two wavelengths 
can be adjusted independently, as can their power levels, polarization and relative phase. The 
disadvantage of separating the beams is that the slight relative phase fluctuations between the paths 
can wash out the grating as it is forming using the seeding technique (11). Despite these phase 
fluctuations, we have still been able to prepare these films for SHG. The output SH signal is 
detected with a photomultplier and averaged with a boxcar integrator. Typically 3000 q-switched 
pulses are averaged and read by the computer. 

We have performed r>vo different types of experiments. In the first experiment, the voltage source is 
applied to the electrode structure, only IR light is coupled through the waveguide and the angle of 
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the electrodes with respect to the optical path is varied in order to optimize the detected SH signal. 
By adjusting the angle of the electrodes, the effective periodicity can be modifled. Once the SH 
signal is optimized, the SHG is measured as a function of applied voltage. 

In the second experiment, the optimized electrodes are kept in place but they are connected to an 
electrometer capable of measuring current as low as lOOaA. The waveguide/electrode assembly and 
associated wires are placed in a metal box to shield them from electromagnetic interference. At this 
point, both IR and green light can be coupled simultaneously into the waveguide. Care was taken to 
insure that they travel along the same region of the waveguide. The current induced in the 
electrodes is measured by the electrometer as a function of the propagation of various optical 
frequencies within the waveguide. 

Results 

Using equation (2) it > ;4S determined that the quasi phase matching period for the 3.5pm 
waveguide is 26.7 pm and for the 2pm waveguide it is 20.2 pm. This was for the lowest order 7 
mode in both the green and IR. The TM modes were used due to the fact that more efficient SHG 
has been obsers’ed in these modes. ( 2) We chose to use the lowest order modes because they have 
the largest overlap integrals with the applied DC field. We found the greatest SHG on the 3.5 
micron thick waveguide using an electrode structure with a 26pm periodicity rotated 7.9 degrees to 
the optical path. This provides an effective electrode periodicity of approximately 26.25pm, very 
close to what we calculated using equation (1). The maximum SHG with 50 volts applied across the 
electrodes was a visible green spot at the correct angle to indicate that it was in the TMp mode. The 
power produced by the SHG was strongly dependent on the periodicity of the waveguide. A rotation 
of less than 1 degree from the optimum caused a 90% decrease in the SHG produced. With the 
electrode periodicity fixed at 7.9 degrees the SHG was measured as a function of applied voltage 
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(Fig 4). The observed quadratic dependence of the electric field induced SHG is expected due to the 
dependence of the conversion efficiency on the square of the electric field in equation 1. 

We were also able to generate a signal in the TM 2 mode in this film by changing the electrode 
periodicity to 70pm, close to the calculated Az value for this mode. The data for SHG in the TM 2 
mode vs applied voltage is plotted in figure S. It is several orders of magnitude lower than what was 
observed for the TMq mode. We also tried to generate the TM} mode but were unable to see any 
signal. This is expected due to the cancellation in the overlap integral. We also tried to produce 
SHG with the input light in the TEq mode, and the electrodes aligned to produce quasi phase 
matching between theTE^ modes in IR and green. For this configuration, no SHG was observ'ed. 
This is most likely due to the fact that the applied electric field was almost perpendicular to the 
electric field of the optical modes. 

With the electrodes optimized for quasi phase matching the TMq mode, we disconnected the voltage 
source and connected the electrometer to the electrodes. We then measured the current induced 
into the electrodes with both the IR and green light coupled into the film, the IR only, the green 
only, and no light. The measured current changed with time at a rate which was dependent on the 
coupling of light through the waveguide. With no light present, an increasing current was 
detected, due to the discharging of the electrodes through the electrometer. When a combination of 
both IR and green light was coupled through the waveguide, the current was obsevred to decrease. 
When only green or IR alone was coupled in, the current increased. This data is plotted in figure 6. 

Discussion 

We have presented a device consisting of a germanium doped glass waveguide and a 
periodic electrode structure which is capable of instantaneously frequency doubling light We have 
demonstrated that a periodic DC electric field is capable of inducing this frequency doubling in glass 






if the periodicity is correct for quasi phase matching. We have been able to selectively excite 
different modes in the SH by varying the electrode periodicity. The efficiency of the electric field 
induced SHG is greatly reduced for higher order modes. This can be explained by examining the 
overlap integrals involved. The conversion efficiency is determined by the overlap integral in 
equation 1. 

The three lowest order modes and the electric field distribution in the waveguide are 
illustrated in figure 7. From this figure, it is observed that the positive and negative lobes in the 
higher order modes will cancel out when integrated across the applied electric field. The applied 
field does decrease as the distance from the electrodes is increased. This will result in a slightly 
nonzero overlap integral for the first order mode, however the separation between the electrodes is 
much larger than the film thickness so the change in electric field across the film is less than 10% 
based on the conformal mapping calculation presented by Ramer [12]. 

The obsen ation of a decreasing current when both the fundamental and SH light are propagating 
through the waveguide, and an increasing current when only fundamental or only SH light is sent 
through strongly supports the SHG models which rely on the buildup of a DC field in the glass. 
When no light is present, a current was observed to flow due to the charging of the electrodes by the 
electrometer. When a combination of green and IR is present, the buildup of an internal electric 
field in the glass causes a current to flow in the electrodes which opposes the existing current and the 
total current is observed to decrease. 

We would like to acknowledge the developmental optical facility at Phillips Laboratory, 
Kirkland Air Force Base, New Mexico for film production, the National Nanofabrication Facility at 
Cornell University where we fabricated the electrode structures and the Air Force office of 
Scientific research and National Science Foundation for sponsoring this research. 
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Partial Demineralization of Canine Mandibular Cortices 

Sanjay Jain 
Medical Student 

Northeastern Ohio Universities College of Medicine 

Abstract 

Partial demineralization is a technique used to process cortical bone 
alloimplants. Obtaining consistency in this procedure while maintaining a 
prescribed level of mechanical strength is the chief concern of this protocol. 
The canine mandible model was used in this study due to its cancellous 
property and application to mandibular reconstuction. Bone segments were 
demineralized with 0.6N HCl over a range of time to determine an appropriate 
time to attain a level of 23%-60% demineralization. A time of approximately 34 
hours was established. A new set of samples were later demineralized in 34 
hours to analyze for relative consistency. 92% of the segments achieved 
partial demineralization within a range of 23%-60%. Despite a linear relation 
with time and degree of demineralization, several other variables exist which 
may influence this process. 
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Partial Demineralization of Canine Mandibular Cortices 

Sanjay Jain 


Introduction 

Cortical bones usually undergo a particular type of processing to 
prepare it as an alloimplant. Demineralization of bone grafts is a technique 
involved with cortical preparation that has been shown to help promote new 
bone formation. It has been theorized that growth factors from the edge of a 
defect are released into the extracellular matrix causing nearby mesenchymal 
cells to differentiate into osteoblasts and chondroblasts.' This 
osteoinduction, as defined by Urist*, is further assisted by the 
osteoconducting nature of a graft which provides a scaffold into which the 
cells migrate. Demineralized bone enhances the osteoconduction by reducing the 
obstacles for migratory activity.^ Completely demineralized bone, although 
more easily incorporated, is not a viable replacement for defects in areas 
requiring a relatively high degree of mechanical strength. Thus it is 
important to obtain a level of demineralization which can provide such 
mechanical capabilities. 

Guo et al* address this issue by investigating partially demineralized 
bone and its affects on mechanical strength. Their protocol used cortical bone 
plates from the dog femora. Using HCl as their solvent, they subjected their 
specimens to various periods of demineralization. They were able to conclude 
that a partial demineralization of 23%- 40% of normal bone mineral content 
provided "reasonable mechanical strength and normal potential for 
osteoinduction". The average time it took to attain this level was between 6- 
10 hours. 

Their study, however, does not take into account several possible 
variables involved with demineralization. These variables include whole 
specimen demineralization, cortical thickness, and biological differences. In 
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our study, our main concern was to establish a protocol that would allow for a 
consistent partial demineralization of bone while maintaining the prescribed 
level mechanical strength for internal fixation and biological function. By 
doing so, we were able to identify many variables which can affect bone 
demineralization. 

In proceeding with this study, we chose to use the canine mandible model 
because of its cancellous property and its similarity to the healing and 
remodeling processes of the tibia.^ It is the ideal model for mandibular 
discontinuity defects, where it is one of the most frequent sites in the 
maxillofacial region requiring bone grafts for reconstruction.* 


Materials and Method 
Preparation 

11 mandibles were harvested, but not aseptically, from mature female 
mixed breed dogs that had been euthanized for use in surgical teaching 
laboratories. The collected mandibles were cleaned of soft tissue and had the 
teeth extracted. The cancellous bone was removed with a high speed bur cooled 
by a saline drip. Each hemimandible was cut into two segments of approximately 
2.5cm. These segments were identified as either the left or right side and 
anterior or posterior segment from a particular dog. These segments were then 
wrapped in saline soaked sponges and stored in a -70 "C freezer to prevent 
dehydration.' 

Bone Analysis 

Bone segments were removed from storage and thawed to room temperature. 

These specimens were later weighed. Mineral content of all segments were 

determined by dual energy x-ray absorbtiometry (DEXA)^. A standardized scan of 

an area of .4 cm' ' s performed in the center of each segment. A scan of the 

t The use of a Lunar DPX model was provided by The Southwest Foundation for 
Biomedical Research 
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entire bone segment was also taken for total bone mineral content. Using 
vernier calipers, the cortical thickness was measured at the anterior, 
posterior, and superior borders of the medial and lateral wall. A mean value 
was then determined for the 6 sites per segment. 

Demineralisation 

Process (for parts I £ II}: Demineralization was achieved by determining 
the amount of acid needed to provide for the total weight of bone. A 
proportion of 100 ml of acid solution / gm of bone was used. The acid 
concentration was 0.6N HCl.^ Each bone specimen used was held in a permeable 
plastic cartridge and placed in the acid solution. A constant circulation was 
maintained with a magnetic stirrer. After removal of the specimen at the 
appropriate time, demineralization was halted by rinsing with deionized water. 
Mineral content of the treated specimens were again determined by DEXA to 
determine the change in bone density. 

Part I: The purpose of this procedure was to determine an appropriate 
time in which to demineralize bone to attain the prescribed range of 23%-40%. 
For this step, 5 dogs were used providing a total of 20 samples. 
Demineralization was performed at one-half hour intervals to a maximum of 5 
hours. 


Part II: Using the data obtained in Part I, we proceeded to demineralize 
bone at a constant time. 6 dogs were used for this procedure giving a total of 
24 samples. One segment from each dog was randomly selected for partial 
demineralization. The corresponding segment from the opposite hemimandible was 
used as the control only. The remaining 2 segments from that particular dog 
(one from each hemimandible) underwent the opposite treatment from its 
counterpart on that side. Thus 12 specimens were demineralized while the 12 
corresponding segments from the opposite hemimandible remained intact so that 
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each demineralized segment had a control segment from the same anatomic 
location on the opposite side. 


Statistical Analysis 

Data collected was compared using a paired ”t" test. A linear regression 
line was determined for the whole bone scan in part I. Correlation 
coefficients were also determined for the data. 

Results 

A scan of the entire bone as well as a smaller portion of the bone 
(.4cm*) showed a large variance in terms of demineralization vs. time in part 
I (see Table 1). The entire bone scan after demineralization within certain 
intervals showed a linear relationship as seen in Fig. 1. The Pearson 
correlation coefficient was -0.8460257 and r‘= 0.7157594. The bone scan 
involving the smaller, select area of bone had a wide array of 
demineralization with time. The Pearson correlation coefficient was 0.11325 
and r*= 0.1282. The data indicated that 3Jj hours was a suitable time to 
demineralize bone to arrive at a BMD level close to 40%. 

Using 12 randomly selected bone specimens for part II, all the specimens 
were subjected to a time of 3^ hours to determine relative consistency. 91.67% 
(11 of 12) of these bones fell close to a range of approximately 23%-60% while 
41.67% (5 of 12) fell close to a range of approximately 23%-40% as determined 
by an entire bone scan by DEXA. Ths smaller scan of .4cm’ showed a greater 
degree of variation similar to that described in part I. The comparison of 
cortical thickness' to demineralization showed no definitive relationship. 
There also was no evident relation of the particular bone piece with the 
demineralization. 


t The mean cortical thickness of all bone samples was 2.24 ± .11 mm, SD = .40 

I 
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Discussion 

Using the canine mandible model, this study shows that partial 
demineralization to an expected percentage is dependent on many factors 
involving a high degree of variability. In our experiments, all samples of 
bone were relatively the same size at 2.5 cm despite studies of self repair 
at this length without a graft.' The hydrogen ions in the medium were always 
kept in excess of what was truly needed. 

Urist et al’ were able to quantitate demineralization at a specific time 
by varying acid concentrations. They calculated a ratio of Ig of compact bone 
to approximately 16.3 ml of 0.6N HCl by basing it on two assumptions: "(a) 
that the chemical reaction Ca|o(P 04 ) 4 (OH); + 14 HCl = 7CaCl, + 3 Ca(H 2 P 04 )} + H.O is 
applicable to bone mineral as a whole, and (b) that bone contains 70% mineral, 
all in the form of hydroxyapatite, Ca,o(PO,) (OH),." They demonstrated that 
during the early stages of demineralization, the bone to acid ratio was very 
important. The more hydrogen ions which were available in the demineralization 
medium, the more the decalcification which would occur. In general, a ratio of 
10 ml HCl/g of bone (6.0 x 10'^ H ions/g bone) partially demineralized bone to 
a level close to 50% at 6 hours. An extrapolation of their results indicate 
that to achieve a range of approximately 23% - 40% demineralization, 15 ml 
HCl/g bone at 6 hours would be sufficient. 

Unlike the previous study which varied acid concentrations, the variable 
that we measured in accordance with demineralization was time. We used a 
shorter time interval than the conventional times used for demineralization to 
account for the thin cortical wall of the mandible. The data indicated a 
linear relationship between time and amount of demineralization. We were 
looking for a suitable time for demineralizatpion of mandibles so that 23% - 
40% could be attained. From our study in part I, our preliminary conclusion 
was that 3J} hours was an appropriate time for part mineralization at a 

higher concentration of 100 ml HCl/g bone. 
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A problem arises though when you standardize time because of other 
variables which play into account. One such variable is the surface area 
exposure. Guo et al^ used bone plates for their demineralization. This would 
insure a consistent flow of HCl throughout the bone surfaces. However most 
implants are performed using larger, intact specimens. We decided to use the 
cortical cribs of the mandible as a whole specimen for true application. As a 
result, the data showed several inconsistencies that is due to the altered 
flow of the free ions to some of the areas of bone. 

The thickness of the bone also plays an important part in the 
demineralization process. Upon harvesting of bone, a bur was used to eliminate 
the cancellous bone. However this creates a rigid surface level and thus 
variances in cortical thickness throughout the wall. This variance can be 
evidenced from a different BMD measurement from one area of bone scanned with 
another area of bone scanned by DEXA. 

Biological variability plays a rather minor role in this study since we 
used similar dogs that were all young females. The calcium content of normal 
cortical bone in dogs is 231.2 mg/g bone.* However the rabbits used by Urist 
et al had a BMC of 270.7 mq one.’ Overall, though, the ratio of inorganic 

to organic content i to 30%. Certain bones in the body may have a 

slightly higher mineral content depending on its functional adaptation. 

Other factors which may have distorted the outcome was the resistance to 
diffusible ions due to the plastic cartridges, an incomplete ionization of 
acid in the medium, and an inaccurate bone scans due to remnants of teeth left 
on the mandible. 

Throughout our trials, we had been using duel x-ray absorbtiometry 
(DEXA)"* as a tool to measure bone mineral content. This instrument is usually 
used clinically for measuring skeletal density in relation to osteoporosis and 
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other disease processes. It is based on the measurement of the relative 
absorption of 2 different energy, x ray beams that pass through the specimen. 
As a result, high resolution images are provided with a discrete reading of 
the bone outline. After scanning an area of interest, DEXA can calculate both 
bone mineral content (BMC) (gm) and bone mineral density (BMD) (gm/cm^ ). An 
example of bone scans by OEXA is shown in Fig. 2. Accuracy is excellent (1%) 
with a correlation between bone mineral content and dry weight of r«0.965 (p< 
. 0001 )." 

Conclusion 

Due to the nature of the thin cortical bone of the mandible versus the 
femur, it is best to maintain a slightly higher level of BMC in order to 
provide enough strength for internal fixation. A broader range of partial 
demineralization involving 23%-60% is recommmended in comparison to the range 
of 23%-40% determined for the femur although further studies should be done 
involving strength tests to actually test this. The disadvantage of decreased 
demineralization is that there would be a slower rate of incorporation of the 
alloimplant with a minimal new bone yield. There are also many other variables 
involved that affect the process of demineralization that must be accounted 
for. A time range can be used as a standard for demineralization as we have 
been able to establish 3*? hours as an appropriate time for mandibles. A 
standard method can also be achieved by varying acid concentrations. In 
general, due to the differences in various bones, a preliminary analysis of 
partial demineralization should be performed on similar samples before 
demineralizing actual samples. 


Bibliography 


1. Hollinger, Jeffrey O., Kleinschmidt, James C.: The Critical Size Defect as 
an Experimental Model To Test Bone Repair Materials. J Craniofacial Surg 
1:1, Jan 1990. 


16-9 





2. Urist, M.R.: Bone Transplants and Implants. In Urist, 

M.R.(ed.):Fundamental and Clinical Physiology of Bone. Philadelphia, 

J.B. Lippincott Co., 1980, p.331. 

3. Nade, S.,Burwell, R.G.: Decalcified Bone as a Substrate For Osteogenesis. 

An appraisal of the interrelation of bone and marrow in combined grafts. 
J Bone Joint Surg. 59B 1977, 189-196. 

4. Guo, M.Z., Zhen, X.S., Lin, L.B.: The Mechanical and Biological Properties 

of Demineralised Cortical Bone Allografts in Animals. J Bone Joint 
Surg[BrJ 73B 1991, 791-794. 

5. Najjar, T.A., Kahn, D.: Comparative Study of Healing and Remodeling in 

Various Bones. J Oral Surg. 35:375, 1977. 

6. Kelly, J.K.: Maxillofacial Missle Wounds: Evaluation of Long-Term Results 

of Rehabilitation and Reconstruction. J Bone Joint Surg. 31:438, 1973. 

7. Friedlander G.E.: Current Concepts Review Bone Banking. J Bone Joint Surg 

64A:2 1982. 

8. Schmitz, J.P., Hollinger, J.O.: The Critical Size Defect as an 

Experimental Model for Craniomandibular Nonunions. Clin. Orthop. 205 
April 1986, 31-40. 

9. Urist M.R., Strates B.S.,: Bone Formation in Implants of Partially and 

Wholly Demineralized Bone Matrix. Clin Orthop 71:271-8, 1970. 

10. Cullum I.D., Ell P.J., Ryder J.P.: X-Ray Dual Photon Absorbtiometery: A 

New Method for the Measurement of Bone Density. Brit J Radiology 
62(739):5a7-592, 1989. 

11. Gluer C.C., et al: Comparative assessment of Dual-Photon Absorbtiometry 

and Dual-Energy Radiography. Radiology 174:223-228, 1990. 


16-10 





Tabla 1. % 


a in Bona Minaral Oanttty (BMOl at VariaUa Timaa 


5 

1 000 

3.060 

1 800 

5 

0 900 

1.100 

2.100 


10 

0.840 

1.020 

1.930 

10 

0.8SS 

1.020 

1 930 

15 

0 738 

0.551 

1 715 

1.5 

0.752 

0.519 

1.830 


> 0.324 

0.159 

2 350 

> 0.838 

0.919 

1.440 

> 0594 

0495 

1 710 


3 5 

0.436 

0.541 

1 950 

3 5 

0 502 

0.478 

2 120 


0.313 

0.127 

2.000 

0.424 

0421 

2.020 

0 187 

0.000 

1735 

0532 

0.478 

1 930 

0 306 

0.245 

1 830 

0 279 

0 171 

1 870 


TaWa 2. % Changa in Bona Minaiai Oa 


(BMDI at 3.5 Houra 


Maan Conical Thicknass immi' 

Partial Bona Scan 

(%) 

Puli Bona Scan 

1 27 

0 240 

0.339 

1 87 

0.096 

0 237 

1 95 

0.486 

0 587 

2 03 

0 402 

0.482 

2 08 

0 429 

0 357 

2 10 

0 222 

0.507 

2.27 

0 571 

0 611 

2.27 

0 160 

0.222 

2 32 

0.556 

0.781 

2 58 

0 218 

0 350 

2 67 

0 369 

0 433 


2 77 


0 367 


0 464 


















PARTIAL DEMINERALIZATION 

Time vs Percentage of BMD 



TIME (hrs) 


Fig. 1. Demineralization of bone analyzed in its entirety by OEXA. 

The linear regression line was defined by ysax-t-b where a=0.9976 and b=-0.1426909. 

(n = 20. wnere p < .0001) 
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